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The Goddard Modeling ar<v: Simulation Facility {6MSF) of the Laboratory for 
Atmospheric Sciences (GLAS) Is engaged in general circulation modeling studies 
related to global atmospheric and oceanographic research. The research activi- 
ties are organized into two disciplines: Global Weather/Observing Systems and 

Climate/Or.ean-Air Interactions. 

The Global Weather activities are grouped in four areas: 1) Analysis and 

Forecast Studies, 2) Satellite Temperature and Wind Retrievals, 3) Analysis and 
Model Development, 4) Atmo;;pheric Dynamics and Diagnostic Studies. During the 
period of this review. Global Weather Research was directed primarily toward 
the development of techniques for the utilization and analysis of the First 
GARP Global Experiment (FGGE) data sets. This involved acquisition and data 
file management efforcs, development of Improved techniques for the analysis 
of FGGE data, observing performance assessments, extended-range forecasts, and 
diagnostic studies. Observing system studies were concerned with the develop- 
ment of a GLAS TIROS-N £ funding retrieval system and preparation for the joint 
NOAA/NASA Advancod Moisture and Temperature Sounder simulation study. Brief 
description?: of several theoretical studies in the areas of analysis, forecast, 
and atmospheric dynamics are also included. 

Climate research activities are focused in three areas: 1) Data Analysis, 

2) Climate Modeling, and 3) Sensitivity Experiments. Ocean/air interaction 
studies are concentrated on the development of models for the prediction of 
upper ocean currents, temperatures, sea state, mixed-layer depths, and upwelling 
zones, and on studies of the interactions of the atmospheric and oceanic circu- 
lation systems on time scales of a month or more. 

This research review contains a synopsis of extended abstracts in the Glo- 
bal Weather/Observinq System studies, and Climate/Ocean-Air Interaction studies 
conducted at the Modeling and Simulation Facility for 1980-81 (Sections 2 and 3). 
Abstracts of the Summer Lecture Series, jointly sponsored by GMSF and the Univer- 
sity of Maryland, are present in Section IV. Section V lists the recent publi- 
cations by the GMSF and visiting scientists. 
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THE EFFECT OF VARYING AMOUNTS OF SATELLITE DATA, OROGRAPHY AND DIABATIC 
PROCESSES ON THE NUMERICAL PREDICTION OF AN INTENSE CYCLONE 

R. Atlas 


Studies of the growth of forecast error for numerical models have resulted 
In quantitative estimates of the amount of error contributed by each of the 
major features of the dynamical prediction system (Robert, 1976). In recent 
years, a number of experiments have been conducted to assess the usefulness uf 
satellite temperature sounding data In weather prediction (Ohring, 1979). 
Objective measures cf forecast accuracy have served as the primary tool for 
evaluating the predictive Impact of the satellite sounding data. In addition, 
a limited number of subjective evaluations and case studies have been per- 
formed. Recently, a detailed subjective evaluation of the effect of horizontal 
resolution on satellite sounding data Impact for the January-March 1976 Data 
Systems Test (DST-6) period was performed (Atlas et al . , 1979). Assimilation 
cycles were carried out from 29 January to 21 February 1976 with coarse (4® 
latitude by 5® longitude) and fine (2.5® latitude by 3® longitude) resolution 
versions of the GLAS second-order general circulation model. For each model 
resolution, an assimilation cycle (referred to as SAT) was run using both con- 
ventional and special DST data, wh'ich included teaiperature soundings from the 
NOAA-4 and Nimbus-6 satellites. A second cycle (referred to as NOSAT) was run 
usin^ the sare data but excluding the satellite-derived temperature soundings. 
Eleven 72 h forecasts utilizing the same resolution In the forecast model as 
In t^e data assimilation were performed at 48 h intervals for each assimilation. 
This evaluation revealed a larger beneficial influence of the satellite sounding 
data with the fine model than had previously been found with the coarse mod'-l 
(Ghil et al . . 1979). In this paper, descriptive results from a study of cyc ione 
evolution In the DST-6 forecast case from 0000 GMT 19 February 1976 will be 
presented. This study was aimed at assessing the effects of satellite data, 
orography and diabatic processes on the numerical prediction of cyclone develop- 
ment and displacement. 


2. SYNOPTIC SITUATION 

At 0000 GMT 19 February 1976, a moderately intense low-pressure system, 
associated with an upper-level short-wave trough, was located off the north- 
west coast of the United States. As the system moved inland, a new low devel- 
oped along an existing stationary front, and became the dominant feature by 
1200 GMT on February 19. During the next 24 h, the low moved southeastward 
and intensified, after which time it recurved and then accelerated toward the 
northeast. This cyclone was accompanied by heavy snow, blizzard conditions, 
and severe thunderstorms and tornadoes. 


3. THE EVOLUTION OF SATELLITE SOUNDING DATA IMPACT 

Fig. 1 shows the sea-level pressure analysis at 0000 GMT 22 February and 
the observed cyclone track. Figs. 2 and 3 depict the 72 h NOSAT and SAT sea 
level pressure forecasts and predicted cyclone tracks that were obtained using 
fine (2.5® lat. by 3® long.) resolution In the assimilation and In the forecast 
model. Comparisons of these charts reveal that a very significant Improvement 
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In the predicted displacement of the surface low resulted from the Inclusion of 
satellite data. The fine model forecast from NOSAT Initial conditions predic- 
ted the low to move southward during the last half of the forecast, while the 
fine model SAT prediction correctly Indicated recurvature of the low to the 
northeast. As a result, there was a more than 1000 km reduction In the vector 
error of the low pressure center's predicted position. 

A study was conducted to determine (1) why the surface lows In the two 
forecasts suddenly diverge and more differently during the latter half of the 
forecast, and (2) If differences In the forcing for the cyclone, at the time the 
two forecasts diverge, can be traced to Initial state differences between the 
SAT and NOSAT systems. To answer these questions, maps and cross-sections of 
a number of primary variables and derived quantities related to the dynamics 
of the cyclone were examined. The synoptic results of this study (presented 
by Atlas et al. , 1979) showed how the prognostic Impact evolved from an enhance- 
ment of the thermal and vorticity fields In the SAT Initial state. As shown 
by Atlas et al. . amplification of the Initial state differences resulted In 
differing phase relationships between the upper-level vorticity maximum and the 
surface low and differing thickness advectlon patterns by 0000 GNT 21 February, 
which produce the diverging paths of the surface low throughout the remainder 
of the forecasts. 

The Initial state differences for this forecast were the result of a time- 
continuous assimilation of satellite temperature sounding data which began on 
29 January 1976. To further understand the role of the sounding data, addi- 
tional fine model forecasts were generated In which varying amounts of satel- 
lite data were excluded from the assimilation. Exclusion of satellite soundings 
after 1800 GNT 18 February produced no noticeable change In the Initial state 
at 0000 GMT 19 February or the 72 h forecast from this time. However, major 
Initial state and prognostic differences were observed when satellite data was 
excluded during the 12 h or 24 h of the assimilation preceding 0000 GMT 19 
February. Figs. 4 and 5 show the 72 h sea-level pressure forecasts from 0000 
GMT 19 February produced by the fine model when satellite data Is excluded from 
the assimilation after 0000 GMT and 1200 GMT 18 February respectively. Compar- 
ison of these figures with (Figs. 2 and 3) shows the progressive Improvement 
of the prediction which results from the assimilation of additional satellite 
data. Examination of the Initial states for these forecasts (charts not shown) 
revealed a gradual enhancement of the thermal vorticity field with the assimi- 
lation of additional satellite data. 


4. EFFECT OF PHYSICAL PROCESSES 

A fine model forecast from SAT Initial conditions was performed In which 
all moist and radiative processes were excluded from the model. Comparison of 
this prediction with the corresponding forecast which utilized full model 
physics, showed the effect of diabatic processes to be very small during the 
first 48 h of the prediction. Following the recurvature of the low, diabatic 
processes became more Important, as warm moist air was advected northward In 
advance of the cyclone and severe convective activity developed. From 48-60 h, 
the rate of Intensification of tlie cyclone In the full physics forecast was 
twice as large as that for the forecast which did not Include diabatic processes. 
During the last 12 h of the forecast this effect did not amplify. Fig. 6 shows 
the 72 h sea-level pressure prognosis for the forecast which excluded diabatic 
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processes. Comparison of this figure with Fig. 3 shows the small diabatic 
effect on the Intensity of the cyclone. We can conclude that moist and radia- 
tive processes In the model did not significantly affect the prediction In this 
case. 


5. EFFECT OF OROGRAPHY 

A fine model forecast from SAT Initial conditions was performed with the 
height of all land above sea level set equal to one meter. Although this 
Introduces an artificial representation of the atmosphere below high terrain. 

It still provides an Indication of what the evolution of the cyclone would be 
If orographic effects were not present. During the first 12 h of the forecast, 
as the cyclone approaches the Rocky Mountains, there are no noticeable differ- 
ences between the prediction with and without orography. After this time, as 
the cyclone crosses the mountains, major differences develop between the predic 
tions, and the role of orography In modifying the evolution of the cyclone 
becomes evident. In particular, the maximum generation of low level cyclonic 
vorticity Is shifted southward, to the lee of the largest slope of the terrain. 
In the run with orography. Fig. 7 shows the 72 h sea level pressure prognosis 
and cyclone track for the fine model forecast without mountains. In addition 
to the more northerly track, the cyclone In this forecast Is larger In scale 
than the cyclone In Fig. 3, 
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Figures 1-7. Sea-level pressure analysis and 72h forecasts for 0000 GMT 22 
February 1976. Cyclone positions at 12h Intervals denoted by 
dots. See text for additional details. 
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THE RELATIVE CONTRIBUTIONS OF INCREASED RESOLUTION IN THE DATA ASSIMILATION 
AND IN THE FORECAST MODEL TO SATELLITE DATA IMPACT 

R. Atlas, M. Ghil , and M. Halem 


The results presented by Atlas et a1, (1979) Indicate that the GLAS model 
forecast accuracy during the January-MarcK 1976 Data Systems Test (DST>6) was 
significantly Improved by the assimilation of satellite sounding data as well 
as by Increased horizontal resolution, and that the practical utility of the 
sounding data was actually enhanced by the Increased resolution. The question 
remains as to whether Increased resolution had a larger effect on the forecast- 
ing skill of the numerical prediction model, or on the assimilation cycle. 
Assimilation cycles were carried out with data from 0000 GMT 29 January to 0300 
GMT 21 February 1976 with two versions of the GLAS second-order 6CM: a coarse 
version with 4® latitude by 5® longitude resolution, called the C -model . and a 
fine version with 2.5® latitude by 3® longitude resolution called the F-model . 
For the two OST-6 cases- where the combined Influence of satellite data and 
model resolution are at a maximum at sea level, I.e., the forecasts from 11 
February and 19 February 1976, the relative contributions of Increased resolu- 
tion In the data assimilation and In the forecast models have been evaluated. 
F-model forecasts were generated from the C-model SAT assimilation Interpolated 
to the F-grId, and C-model forecasts were generated from the F-model SAT assi- 
milation Interpolated to the C-qrld. We shall denote the former by CFS and the 
latter by FCS. These forecasts were then compared with the corresponding fore- 
casts which had utilized the same grid resolution In the data assimilation and 
forecast models, CS and FS. 

Fig. 1 presents the S] scores for the sea level pressure forecasts from 
19 February 1976 as verified over North America at 12 h Intervals; Fig. 2 shows 
the S] scores for the 11 February case. The combined effect of Increased resol- 
ution In the data assimilation and In the forecast model Is represented by the 
difference in S] between the F-model forecasts from the F-assImilatlon (FS) and 
the C-model forecasts from the C-assIrallatlon (CS). The Impact of Increased 
resolution In the data assimilation alone Is obtained by subtracting the S] 
scores for the C-model forecasts from the F-assImllation (FCS) from the corres- 
ponding S] scores for the C-model forecasts from the C-assImllation (CS). The 
Impact of Increased resolution in the forecast model, but not In the assimila- 
tion, Is obtained In a similar manner by subtracting the Si scores for the 
F-model forecasts from the C-assImllation (CFS) from the S] scores for the 
C-model forecast from the C-assImllation (CS). 

Comparison of the forecasts from 19 February (Fig. 1) indicates that both 
Increased resolution In the assimilation and increased resolution in the fore- 
cast model contributed to the improvement of forecast accuracy at most verifi- 
cation tli.’es; the effects, however, are not strictly additive. The Improvement 
of forecast accuracy resulting from Increased resolution In both the assimila- 
tion and forecast model Is larger than the sum of the Impacts of Increased 
resolution In the assimilation and increased resolution In the forecast model 
at all verification time periods. During the first 24 h. Increased assimilation 
resolution results in a small improvement In accuracy, while Increased forecast 
model resolution results In a slight degradation of forecast skill. At 36 h, 
both effects are beneficial and of comparable magnitude. From 48 h to 60 h. 



the beneficial effect of forecast model resolution is much larger than the 
effect of assimilation resolution. At 72 h, they are once again of similar 
magnitude. 

The verification for the forecasts from 11 February (Fig. 2) shows 
what difference results. In this case, the effect of increased assimilation 
resolution (CS-FCS) and that of increased forecast model resolution (CS-CFS) 
taken separately are both slightly negative and of nearly identical magnitude 
at 12 h. In contrast, the forecast which utilized increased resolution in 
both the assimilation cycle and in the forecast model shows a modest improve- 
ment at this time. At 24 h, there is a negligible impact of increased resolu- 
tion in the assimilation, but a modest beneficial impact of increased forecast 
model resolution. The latter impact by itself is slightly larger than the 
combined effect of increased resolution in both the assimilation and the fore- 
cast model. From 36 h to 60 h, the beneficial influence of increased forecast 
model resolution is much larger than the improvement resulting from increased 
resolution in the assimilation and nearly as large as the combined effect of 
increased resolution in both the assimilation and forecast model. At 72 h, 
there is a negligible influence of increased assimilation resolution, while 
the improvement of increased forecast model is significant and larger than the 
combined Impact. 

Our discussion of objective impact measures indicates that both increased 
resolution in the assimilation and increased resolution in the forecast model 
can act to improve forecast accuracy. In the 19 February case, each effect was 
important although the influence of forecast model resolution dominated at most 
time periods. In the 11 February case, the effect of forecast model resolution 
clearly was responsible for the large increase in accuracy between the C- and 
F-model results. 

4 subjective evaluation of the forecasts for these cases was conducted in 
order to determine if similar results would be found for the improvement of 
practical forecast utility. Comparison of the sea level pressure prognostic 
charts for each of the forecasts from 19 February showed the forecast which 
utilized increased resolution in both the assimilation and forecast model to 
be superior to all other forecasts at each time period, in agreement with the 
S] verification. However, the C-model forecast from the F-assimilation (FCS) 
was substantially less useful than the F-model forecast from the C-assimilation 
(CFS) at all time periods after 12 h. In particular, the track and Intensity 
of the cyclone over the United States was better predicted by CFS. The subjec- 
tive evaluation of the forecasts from 11 February was in close agreement with 
the S] verification. Comparison of these forecasts confirmed a greater benefi- 
cial influence of increased resolution in the forecast model. 
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Fig. 1 S] score verification for the sea-level pressure forecasts for North 
America from 0000 GMT 19 February 1976. Results are presented for 
the coarse-model forecasts from coarse SAT initial conditions (CS) 
and fine SAT initial conditions (FCS), and for the fine model forecasts 
from fine SAT initial conditions (FS) and coarse SAT initial conditions 
(CFS). 



Fig. 2 Same as Fig. 1 for the forecasts from 0000 GMT 11 February 1976. 
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THE EFFECTS OF TROPICAL WIND DATA ON THE PREDICTION OF ULTRALONG WAVES 

U. E. Baker 


Summary: This paper is a preliminary study of the influence of tropical wind 

data on the prediction of planetary waves in which two assimilation experiments 
are performed, one with and one without FG6E tropical winds. The planetary 
wave error is then analyzed in 72 h forecasts from the initial conditions pro- 
vided by the two assimilations. 


1. INTRODUCTION 

Inaccurate specification of the initial conditions for numerical models is 
recognized as an important source of forecast error. Of particular concern is 
the contribution of the planetary wave error variance to the total error variance, 
estimated to be as much as 50% in the first two planetary waves by Baumhefner and 
Bettqe (1981). The augmented observational database provided by FGGE offers a 
unique opportunity to further examine the role of initial data in numerical 
weather prediction, especially the effects of the data on planetary wave fore- 
casts. 


Recent model studies (Paeqle, 1978) indicate that tropical disturbances 
may influence the subtropical jet streams, while analyses of a DST-6 assimila- 
tion experiment (Paeqle et al. , 1979) suggest that low-latitude convective 
activity may propagate energy outside the tropics on short time scales. Kung 
and Burqdorf (1978) also found the A-B scale disturbances over the GATE region 
to be a source of eddy potential energy outside the region. In a simulation 
studv (Gordon et a 1 . , 1972), the insertion of tropical wind data resulted in a 
reduction of the wind error in the extratropics. Recently, Somerville (1980) 
has suggested that tropical data may affect the forecast of ultralong waves. 

In the present study we examine the effects of tropical wind data from FGGE 
on the numerical prediction of ultralong waves by performing the real -data 
counterpart to the simulation study of Kasahara and Williamson (1972), which 
indicated a preferred global data coverage of winds in the tropics and tempera- 
tures outside the tropics, in agreement with geostrophic adjustment theory 
(Blumen, 1972). Using the Goddard Laboratory for Atmospheric Sciences (GLAS) 

4th order general circulation model (Kalnay-Ri vas and Hoitsma, 1979), two 
assimilation experiments are conducted. Experiment 1 utilizes wind data in the 
tropics (20®S to 20®N) and temperature data outside the tropics. Experiment 2 
is the same as Experiment 1 except that no wind data are assimilated. Both 
experiments begin from the initial conditions of 0000 GMT 9 January 1979 provided 
by the 6 h first guess of the GLAS FGGE assimilation cycle. The data are then 
assimilated for a period of several days followed by a series of 72 h forecasts 
generated from initial conditions provided by each experiment. In the following 
sections, we describe the assimilation/forecast model, the FGGE data utilized 
in these experiments, the objective analysis and assimilation procedure, and 
the differences in the ultra-long wave predictions. 
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2. OBJECTIVE ANALYSIS AND ASSIMILATION OF THE F6GE DATA 

In this section the objective analysis and assimilation procedure utilized 
with the FGGE data Is described. 

2.1 The Asslmllatlon/Forecast Model 

The model utilized In this study Is the global 4th order 6LAS general circul- 
ation model described in detail In Kalnay-Rivas et al. (1977) and Kalnay-Rivas and 
Holtsma (1979). The model Is based on an energy conserving scheme with horizontal 
differences computed with 4th order accuracy. A 16th order Shapiro (1970) filter 
Is applied every 2 h to the sea level pressure, wind, and potential temperature 
fields. There are 9 vertical layers equal In sigma with a uniform horizontal 
grid (4® In latitude by 5* in longitude). A non-staggered horizontal grid, the 
Arakawa (1972) A grid. Is used In the current model. The Arakawa B grid was 
employed In the earlier 2nd order version of the model (Somerville etaK, 1974). 
With the exception of the computation of the longwave radiation (Wu, 1980), the 
parameterization of the physical processes In the 4th order model Is substantially 
the same as In the 2nd order model. 

2.2 The Data 

For this study we utilize the data collected during the first FGGE Special 
Observing Period (SOP-1) from 9-17 January 1979. This period was selected because 
of the presence of strong high-level divergence In the tropical Pacific during 
'Ms time (Paegle, 1981). Table 1 summarizes the data assimilated In the two 
experiments. 

Table 1. Tropical (20®S to 20®N) and extratropical FGGE data assimilated from 9-17 
January 1979. Experiment 1 is denoted by a "1" and Experiment 2 by a "2". 
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As may be seen in Table 1, the tropical wind database Is dominated by the 
cloud-track winds (denoted by CTW) from five geosynchronous satellites GOES 
East* West* and Indian Ocean* METEOSAT* and the Japanese satellite). The cloud- 
track winds are* therefore* likely to be the data primarily responsible for the 
differences In the numerical predictions discussed In Section 3. 

2.3 The Objective Analysis Procedure 

In the 6LAS objective analysis scheme (Baker et a1. * 1981)* eastward and 
northward wind components* geopotential height and relative humidity are analyzed 
on mandatory pressure surfaces. The 6 h model forecast provides a first guess 
for these fields at 300 mb and sea level* where sea level pressure and sea level 
temperature are also analyzed. The first guess for the other levels Is obtained 
from the model first guess* modified by a vertical Interpolation between the two 
closest completed analyses. Vertical consistency Is maintained through static 
stability constraints. The analysis at each level Is performed with a succes- 
sive correction method (Cressman* 1959) modified to account for differences In 
the data density and the statistical estimates of the error structure of the 
observations. The average distance between data points Is found In a circle 
with a radius of 800 km centered at each grid point. Three scans are performed 
with a radius of Influence R-| ° C'jd* where the coefficients c^ (1.6* 1.4* 1.2) 
were chosen to minimize the analysis error (Stephens and Stitt* 1970). However* 
the radius of Influence Is not allowed to become smaller than 300 km. During 
this process* all data are checked for horizontal consistency. The completed 
analyses are smoothed and then Interpolated to the model sigma levels. 

The assimilation procedure provides for the Intermittent analysis of batches 
of data grouped In a ^ 3 h window about each synoptic time. In these experiments, 
the wind and height fTelds are analyzed independently with no explicit coupling 
or balancing. However* a geostrophic correction Is applied to the first guess 
wind field every 6 h in the extratropics using a technique similar to that pro- 
posed by Kistler and McPherson (1975). The correction Is computed from the 
change to the mass field (surface pressure) only. 


3. RESULTS 

Only a limited number of forecasts have been completed, and only the results 
for the 72 h 500 mb forecast from 0000 GMT 15 January 1979 are presented here. 

In Fig. 1* the 72 h planetary wave (zonal wavenumbers 1-3) error (analysis-forecast 
difference), as verified against the NMC global O/I analysis on 0000 GMT 18 
January 1979* is shown for the western half of the northern hemisphere. Fig. la 
Illustrates the planetary wave forecast error from the initial conditions pro- 
vided by Experiment 2 (no wind data assimilated)* and Fig. lb Illustrates that 
for Experiment 1 (tropical wind data assimilated). The reduction in the error 
In Fig. lb is striking, particularly between 40*’N to 50®N, along the west coast 
of North America. There Is also an eastward shift in the maximum positive error. 
The wavenumber 2 dominance of the error field at 40**N Is noteworthy. A smaller 
reduction in the error over the Atlantic may also be seen. The ultralong wave 
error is also reduced over Eurasia (not shown) with the assimilation of tropi- 
cal winds by an amount comparable to that seen over the Atlantic. 

The strong divergence in the tropical Pacific during the time period of 
this study (Paegle* 1981) may be responsible In part for the differences In the 
ultralong wave predictions. This possibility Is being explored. Additi lal 
cases are being analyzed and will be reported c i at the Symposium. 
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Fig. 1. The 500 mb planetary wave error for the 72 h forecast 
from 0000 GMT 15 January 1979. The verification is the 
NMC global O/I analysis for 0000 GMT 18 January 1979. 
The forecast error for Experiment 2 (no wind data) is 
shown in la and that for Experiment 1 (tropical wind 
data) in lb. The contour interval is 30 m. 
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THE STATE OF THE ATMOSPHERE AS INFERRED FROM THE FGGE 
SATELLITE OBSERVING SYSTEMS DURING SOP-l^ 

M. Halera, E. Kalnay, W. E. Baker, and R. Atlas 


1. INTRODUCTION 

At the Goddard Laboratory for Atmospheric Sciences (GLAS), we performed a 
series of data assimilation experiments to test the influence of different ele- 
ments of the satellite observinq systems that were collected during the Special 
' Observing Period (SOP-1). In this paper we present results from some of the 
experiments. These findings show that the FGGE satellite systems are able to 
infer the three-dimensional motion field and improve the representation of the 
large-scale state of the atmosphere. We also present some preliminary results 
of the forecast impact of the FGGE data sets. 

The GLAS analysis/forecast system for producing a global gridded analysis 
consists of an objective analysis scheme making use of the continuity provided 
by a model first guess integrated from a previous analysis. The first guess 
is then corrected by data collected in a > 3 hour window about each analysis 
time. The results are obtained with the (>LAS analysis cycle (Baker et al. , 
1931). 

The assimilation/forecast model is the 4th order global atmospheric model 
described in Kalnay-Rivas et al . (1977) and Kalnay-Rivas and Hoitsma (1979). 

It is based on an energy conserving scheme with all horizontal differences com- 
puted with 4th order accuracy combined with the periodic use of a 16th order 
Shapiro filter. 


2. COVERAGE AND ACCURACY OF THE FGGE SATELLITE DATA 

There are approximately 6000 temperature soundings per day produced opera- 
tionally from the TIROS-N HIRS/MSU sounders, 1000 temperature soundings from 
the NOAA-4 VTPR sounders, and roughly 6000 cloud-tracked winds per day from 
the five FGGE geostationary satellites. The influence of this large volume of 
satellite data on a global analysis system is one of the main concerns of this 
study. In this paper we briefly describe only the statistical properties and 
coverage of the temperature sounding data leaving the other data sets such as 
geostationary cloud-track wind motions, drifting buoys, etc., for a later 
detailed paper. 

2.1 TIROS-N Sounding Data 

■ Fig. la shows the geographical frequency distribution of TIROS-N temperature 

\ soundings observations per grid area per day for the period January 5 through 

I January 21, 1979. If at least one sounding occurs in an area within ^ 3 hours 

I 1 Published in Condensed Papers and Meeting Report: International Conference 

t on Early Results of FGGE and Large-Scale Aspects of its Monsoon Experiments, 

I Jan. 12-17, 1981, Tallahassee, Florida. Published by the World Meteorological 

j Organization, Geneva, Switzerland. 
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of an analysis time* then that area is considered to have been observed. Thus, 
the maxiinuin possible frequency of observations per day for any area is four. 

The frequency distribution shown in Fiq, la indicates that the tropics are beinq 
observed on the averaqe of once every other day by TIROS-N while the extratropics 
are qenerally observed more than once per day but less than twice per day. It 
is also seen that operational TIROS-N soundings are not produced over mountainous 
regions, such as the Rockies and Himalayas, although they are generated over 
other land areas. Fig, lb shows the combined frequency per day per grid area 
for the two satellite systems during FGGE, namely NOAA-4 and TIROS-N, Over 
oceans, the frequency is usually two to three observations per day. Although 
this frequency of observations still does not match the 6 h frequency at which 
analyses are produced operationally at NMC, i.e. four times daily, the asynoptic 
coverage is comparable to the nearly twice daily coverage of the United States 
rawinsonde network. 

The accuracy of the satellite observing systems is a more difficult pro- 
blem to assess because of the lack of “ground truth" for verification over data- 
sparse regions. Two methods that have been used to estimate temperature sounding 
errors appear to present differing accuracy estimates. The methods are: 1) rms 

differences of layer mean temperatures obtained from collocating rawinsonde 
observations and satellite temperature profiles in space and time, (e.g. Phillips, 
et al., 1979), and 2) rms differences of layer mean satellite temperature sound- 
ings and the NMC analysis (e.g. Schlatter, 1981). Phillips' statistics show the 
rms differences at mandatory pressure levels for the TIROS-N type A soundings 
at selected ocean stations to be ~ 2.2®C in the extratropical troposphere and 
l.S'^C in the tropics during SOP-1. The errors above 300 mb are larger, up to 
3.0®C at 200 mb. Schlatter's statistics indicate a surface to 50 mb rms differ- 
ence in thickness temperature to be ~ 0.8°C over the United States. He ascribes 
this higher accuracy to the fact that small biases of opposite sign in the 
lower and upper tropospheric sounding temperatures cancel when integrated over 
large thicknesses. 

The GLAS analysis system attempts to take advantage of the accuracy in the 
sounding thicknesses by performing a geopotential height analysis. Thus, we 
show in Fig. 2a and Fiq, 2b the vertical rms differences between the final 
height analysis and the rawinsonde profiles, TIROS-N clear-colMinn (type A) 
soundings, and the microwave soundings for land and ocean, respectively. The 
TIROS-N type A soundings over land are seen to be in slightly less agreement 
with the analysis than the rawinsondes up to 300 mb and in better agreement 
with the analysis above 300 mb. The rms differences in "microwave only" (type 
C) soundings over land consistently increase from 850 mb to 300 mb and then 
decrease above 300 mb. The surprising agreement for the type A soundings is 
not only consistent with Schlatter's findings but also shows the geopotential 
height to be accurate at the lower levels as well. Another factor contributing 
to these statistics may be that the rawinsonde captures more of the sub-grid 
scale variance and this variance cannot be represented on a coarse 4° x 5° 
analysis. Additionally, the rawinsondes utilized in this study were not cor- 
rected for the effects of solar radiation. Over the oceans, the rms differences 
in Fiq, 2b show that both infrared and microwave TIROS-N retrievals fit the 
analyzed height better than the rawinsondes. This may be explained by the fact 
that the satellite soundings are the dominating data influencing the analyses 
over the oceans. Nevertheless, the larger differences in the microwave soundings 
over land and oceans, which constitute 50% of the retrievals, must be spuriously 
affecting the analysis. 
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3. EXPERIMENTAL ANALYSIS/FORECAbT SYSTEM 

Several experiments to assess the Influence of the F6GE observinq systems 
on the qlobal analysis were performed for the period January 5 to March 5, 1979. 
Each experiment is characterized by a 6 h analysis/forecast cycle using a speci- 
fied subset of the FGGE Il-b data at each synoptic time. We confine ourselves 
in this section mainly to the study of the effect of satellite observations 
on the analysis. Table 1 describes the acronyms used to denote the different 
experiments and the combinations of FGGE data that were utilized in the assimi- 
lation cycles. The initial state for all of the assimilations was the January 
5, 1979 final analysis provided by the National Meteorological Center (NMC). 

3.1 Rms Differences Between First Guess and Analysis 

A measure of the influence of the observational data on an analysis cycle 
can be obtained by calculating the rms correction made between the forecast 
first quess and the completed analysis. In Figs. 3-5, we present the rms dif- 
ference fields between the 300 mb 6 h forecast (i.e. first guess field) and the 
final analysis for the NOSAT, FGGE, and NORAOB experiments, respectively, for 
the period January 5 to January 21, 1979. The shaded areas in these figures 
indicate regions where the average difference during the 16 d period is less 
than 20 m; diagonal hatching indicates large errors greater than 60 m; blank 
areas correspond to intermediate values between 20 and 60 m. 

Fig. 3a shows that the 6 h rms forecast error over the vast oceanic regions 
in both hemispheres is less than 20 m, while the differences over regions with 
dense rawinsonde observations are greater than 60 m. The observed structure 
in the 6 h forecast error can be explained by the fact that the 6 h forecast 
itself defines the atmosphere in data-sparse regions leading to lower "errors." 
Those small corrections are just an indication of the serious data gaps of the 
conventional network. We note that large 6 h errors are present in the Southern 
Hemisphere and in the tropics. By comparing with Fig. 4a which shows the dis- 
tribution of rawinsondes, it is clear that the large corrections occur in data- 
sparse regions at isolated ocean and island stations and can be more than 120 m. 

The rms differences are smaller over Europe than over North America, probably 
due to more frequent rawinsonde observations. It is interesting to note the 
maximum 6 forecast error of 60 to 80 m which occurs along the full length of 
the western boundary of the North American continent, as well as in other regions 
downstream of data-sparse regions. Along the northwest coast of Europe this 
sharp error gradient does not appear because of the presence of rawinsonde ship 
stations in the Atlantic. 

Fig, 3b depicts the FGGE 6 h forecast rms error in the 300 mb height field. 
Comparison with Fig. 3a reveals that large errors have been significantly 
reduced, especially along coastal boundaries and at the isolated oceanic sta- 
tions. For example, the large errors along the west coast of North America 
have been reduced to about 40 m. On the other hand, we now observe increases 
in areas of large spatial extent such as over the northern Pacific ranging from 
40 m to 60 m. This indicates that satellites and other FGGE observing systems 
have eliminated some of the data gaps in the NOSAT system and are introducing 
substantial modifications to the 6 h forecast fields used in the analysis. 
Moreover, the closer agreement of the rawinsonde heights to the 6 h forecast 
in the FGGE experiment suggests that the satellite data over the ocean have 
produced a more accurate initial state. 
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The rms difference between the final F6GE and NOSAT analyses is shown in 
Fiq. 4b. In the data-sparse regions, the differences range from 80 m over the 
Pacific to more than 200 m in the Southern Hemisphere. The differences in the 
tropics and oceanic regions of the Southern Hemisphere are as large as those 
between randomly chosen synoptic situations. Only over the North American and 
Eurasian continents are the differences between the final analyses smaller 
than 20 m. 

The sharp contrast in the difference fields between the land and sea in 
the two analyses raises a very important question. Is it possible that the 
inability of present models to accurately forecast the planetary scale waves 
may be partly due to the land-sea contrast in data coverage? It might be 
that large errors over the oceans give rise to spurious planetary scale wave 
components in the initial conditions which incorrectly travel during the 
forecast, adversely af fleeting the forecast skill. 

To test the hypothesis that the satellite data alone can produce analyses 
which are consistent with the g^"X>sphere . we performed an additional exper- 
iment without rawinsonde data (HORAOB). Figs. 5a and 5b illustrate the 300 mb 
rms height differences between the 6 h forecast and N0RA06 analysis and the 
difference between the FGGE and NORAOB analyses, respectively. The low rms 
error in Fiq. 5b over the oceans is due to the predominance of the same satel- 
lite data in both experiments. This is analogous to the small differences we 
noted in Fig. 4b over land, where both experiments are dominated by rawinsondes. 
The largest differences occur downstream of the major mountain chains . For 
example, the rms differences are smaller over the western half of the United 
States and Canada than over the eastern regions of North America. This results 
from the fact that neither temperature sounding data nor cloud-tracked wind 
data were available operationally over the Rocky Mountains as indicated in the 
data coverage maps of Fig. 1. Around isolated island stations, there is good 
agreement, even without rawinsonde data. 


4. SYNOPTIC SCALE COMPARISONS 

In this section we compare the NORAOB inferred height and wind fields with 
the FGGE analyses over data dense regions. The purpose of such a comparison 
is to assess the quality of the analysis over the oceans where little or no 
data is available for verification. 

Figs. 6a and 6b present typical comparisons of the vertical structure of 
inferred wind profiles at two low-latitude isolated ocean stations with coinci- 
dent rawinsondes. The panels on the left of the figures depict the wind direc- 
tion as a function of pressure, and the panels on the right illustrate the wind 
speed. Also shown are the aircraft and cloud-tracked wind data available within 
400 km of the station at 3 h of analysis time. In both cases the NORAOB 
analysis defines the vertical wind speeds accurately at low and upper levels 
but underestimates the wind speed at 500 mb by about 10 to 20%. This is gener- 
ally the case with many other profile comparisons we have examined. We also 
see in Fig. 6b that the NORAOB analysis can infer strong wind shear, change 
abruptly within a 100 mb interval from easterly to westerly, or exhibit multiple 
directional changes with height in a light wind situation as seen in Fig. 6a. 


Flqs. 7a and 7b show the 300 mb analysis of the qeopotentlal height field 
for the F66E and NORAOB experiments. The position of the trough In the NORAOB 
experiment Is In excellent agreement with the FGGE analysis, although Its 
Intensity Is slightly weaker. Part of the weakness In the trough Is due to the 
lack of satellite data over the Rockies (see Figs, la and lb). The difference 
between the FGGE and NORAOB analyses Is small over the central United States, 
with the NORAOB analysis underestimating the ridge In the northeast. Figs. Ba 
and 8b show the corresponding 300 mb wind vector and Isotach maps. The position 
of the Isotach minimum and maximum Is nearly Identical as Is the flow direction. 
The Intensity of the system Is, however, systematically underestimated over most 
of the United States by about 10% exceot In the Intense jet streak located over 
the Gulf of Mexico, where It Is nearly 40T. weaker. This difference Is traced 
to a weaker ageostrophic wind component In the NORAOB analysis. The absolute 
vorticity, vorticity advectlon, and thermal advectlon fields have a similar 
structure and position but weaker Intensity. Whether this deficiency can be 
corrected over the western part of the United States through the use of cloud- 
tracked wind data or the extraction of temperature sounding data over high 
terrain Is still a conjecture at this time. Nevertheless, the ability to 
capture much of the day-to-day detailed changes In the temperature and wind 
fields is highly encouraging. 


5. FORECAST IMPACT RESULTS 

Clearly the most Important application of an objective analysis and. at 
the same time, the most sensitive test of Its accuracy Is the Impact It has on 
forecasting skill. We have performed a series of short-term forecast Impact 
tests with our coarse resolution model and a few selected extended range fore- 
casts with a higher (2.5° latitude by 3° longitude) resolution model. The 
statistical sample of fourteen forecasts Is small, but examining a number of 
geographically distinct regions Increases the Independent sample size, making 
the results statistically significant. Each forecast was of five days duration 
from Initial conditions taken every fourth day from January 9 through March 2 
from the analysis cycles described earlier. Figs. 9a and 9b present the SI 
skill score for the three day forecasts evaluated relative to the NMC operational 
analysis over North America, Europe, and Australia. A consistently significant 
positive impact of the FGGE data occurs over Australia In the sea level pressure 
and 500 mb SI skill scores. A smaller and consistent Impact occurs over Europe. 
There Is a small positive impact over North America at 500 mb, but a less con- 
sistent Impact In the sea level pressure field. Also shown In these figures 
are the SI scores for the forecasts from the NORAOB experiment. Over Australia, 
the accuracy of the three day forecasts Is comparable to that of the FGGE 
system at 500 mb but slightly less accurate at sea level. Over Europe and 
North America for the limited number of cases, the NORAOB forecasts are, in 
general, poorer by about 24 h in forecast skill. 


6. CONCLUSIONS 

The major findings of this study are: 

(1) The rms geopotential height difference over land between the dear- 
column infrared TIROS-N retrievals and the completed analysis Is only slightly 
larger than the rms difference between the rawinsonde height and the analysis 
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up to 300 mb. Above 300 mb, the completed analysis is substantially more con* 
si stent with the qeopotential height of the infrared retrievals than with the 
rawinsondes. Microwave sounding differences over land and oceans are substan- 
tially poorer than infrared retrievals. 

{?.) The 6 h forecast error over land downstream of data-sparse regions 
is significantly reduced in the F6GE analysis as compared with the NOSAT 
analysis. 

(3) The large-scale synoptic structure of the global geopotential height 
field produced by the NORAOB and FGGE systems is remarkably similar. The only 
exceptions occur in regions over and downstream of major mountain chains where 
soundings are not available. Similarly, the instantaneous wind vector field 
inferred from the NORAOB experiments over land is in good agreement with the 
FGGE analysis but has weaker intensity. 

(4) The skill of the forecasts from initial conditions obtained with the 
NORAOB experiment are somewhat degraded when compared with forecasts from the 
full FGGE system. 

(5) The forecast impact from the FGGE system over Australia is consis- 
tently significant and positive. The impact over North America and Europe is 
also consistent but slightly weaker at 500 mb and less consistent in sea 
level pressure over North America. 

A serious limitation of the NORAOB system is the lack of observational 
soundings over mountains and cloud-tracked wind vectors over land. 
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DATA UTILIZED FROM THE FG6E OBSERVING SYSTEM 


EXPERIMENT 

FGGE 

NOSAT 

NORAOB 

NOGEOS 

NO TIROS 

NO BUOYS 

Rawinsondes • 2A 

X 

X 


X 

X 

X 

NAVAIDS 

X 


X 

X 

X 

X 

Dropwi ndsondes 

X 


X 

X 

X 

X 

Aircraft - 2A 

X 

X 

X 

X 

X 

X 

ASDAR 

X 


X 


X 

X 

TIROS-N 

X 


X 

X 


X 

VTPR 

X 


X 

X 


X 

NESS CTW 

X 


X 


X 

X 

European CTW 

X 


X 


X 

X 

Japanese CTW 

X 


X 


X 

X 

Wisconsin CTW 

X 


X 


X 

X 

(East/West) 







Wisconsin CTW 

X 


X 


X 

X 

(Indian Ocean) 







TWERLE 

X 


X 


X 

X 

(const, lev. balloon 







Surface stations 

X 

X 

X 

X 

X 

X 

(land) - 2A 







Ships - 2A 

X 

X 

X 

X 

X 

X 

Drifting buoys 

X 


X 

X 

X 





26 











TABLE 1 





i 


ORIGINAL PAGE 18 
OF POOR QUALIIY 






ORIGINAL PAGE IS 
OF POOR QUALITY 


Fig. lb. COMBINED TIROS-N AND NOAA-4 TEMPERATURE SOUNDINGS (FREQ. /DAY) 
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Fig. 3a. RMS DIFFERENCE BETWEEN 6 hr FORECAST AND NOSAT ANALYSIS 
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Fig. 3b. RMS DIFFERENCE BETWEEN 6 hr FORECAST AND FGGE ANALYSIS 
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Fig. 4b. RMS DIFFERENCE BETWEEN NOSAT AND FGGE ANALYSES 
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Fig. 5b. RMS DIFFERENCE BETWEEN NORAOB AND FGGE ANALYSES 
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Fig. 5. WIND PROFILE' FROM NORAOB ANALYSIS COMPARED WITH RAWINSONDES 
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NUMERICAL PREniCTION OF THE LAROE SCALE TROPICAL FLOU 
E. Kalnay, W. E. Baker and J. Shukla 


1. INTRODUCTION 

The problem of numerical forecasting is more difficult in the tropics than 
in mid- and hiqh-latitudes for the following reasons: 1) There has always been 

a lack of reliable large scale observations in the tropics. This is partly due 
to the low land coverage in the tropics, and more importantly, to the need to 
measure the wind field rather than the mass field, which results in a challenging 
problem of instrumental accuracy. 2) The fact that the dominant instability in 
the tropics is convective rather than baroclinic as in the extratropics. This 
implies very fast growth rates, and requires an accurate representation of sub- 
grid scale forcing, especially convective heating, which is one of the most 
difficult problems of atmospheric modeling. 3) The "no skill competition" 
persistence, against which numerical forecasts are compared, is more accurate 
in the tropics than in the extratropics, and therefore more difficult to be 
improved upon (Shukla, 1R81). 

During the Global Weather Experiment (GWE or FG6E) year of 1979, a special 
effort was made to gather data in the tropics including unconventional observ- 
ing systems such as satellite temperatures and cloud-track winds, dropwindsondes , 
research aircraft, etc. 

In this study we examine the mean error characteristics of two series of 
forecasts of the tropical flow. By using as initial conditions analyses made 
with and without the FGGE special observing system, we estimate the impact that 
initial data has on the accuracy of the forecasts. 


2. DESCRIPTION OF THE EXPERIMENTS 

The results reported here are obtained from the global assimilation and 
forecast experiments performed by Halem et al . (1982) for the purpose of assess- 
ing the impact of satellite data upon extratropical analysis and forecasting. 

The GLAS analysis/forecast system for producing a global gridded analysis 
consists of an objective analysis scheme which makes use of the continuity 
provided by a first guess which is a 6 hour forecast from the previous analy- 
sis. The first guess is then corrected by all the data collected within a _+ 3 
hour window about each analysis ’ime. The analysis scheme (Baker et al . , 1981) 
is a modified successive correction method (Cressman, 1959) which takes into 
account the density and the quality of the observations. The model used in 
both the analysis cycle and the forecast is the GLAS fourth-order global atmo- 
spheric model described in Kalnay-Rivas et al . (1977) and Kalnay-Rivas and 
Hoitsma (1979). It is based on an energy conserving scheme with all horizontal 
differences computed with fourth-order accuracy. A 16th order Shapiro filter 
is applied periodically to remove unresolved scales. The parameterization of 
subqrid physical processes is identical to that of the GLAS climate model 
(Shukla et al . , 1981). It includes long and short wave radiation with a diurnal 
cycle, which allows a convective cloud parameterization, conditional instability 
supersaturation clouds, a bulk formula parameterization of surface fluxes and 
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a realistic orography. The resolution used in these experiments, 4° latitude 
X S'' longitude and 9 vertical levels, is somewhat coarse, but this is partly 
compensated by the improved accuracy of the finite differences used in the 
model . 


Two analysis cycles were performed for the first FG6E Special Observing 
Period (SOP-I), from January 5 to March 5, 1979. In one of them, denoted FGGE, 
all available FGGE Il-b data were assimilated. In the second experiment, 
denoted NOSAT, only conventional data (rawinsonde, pilot balloon, aircraft and 
surface land and sh’p reports) were utilized (Table 1, from Halem et al . , 1982). 

Fourteen 5-day numerical forecasts were then generated every four days 
from the initial conditions of both the FGGE and the NOSAT assimilation experi- 
ments. 


3. RESULTS 

We present here preliminary comparisons of the mean and standard deviations 
of the forecast error. The mean errors represent the systematic forecast errors 
which may be due either to the parameterization of forcing, or to systematic 
observational errors or lack of data. The standard deviation of the forecast 
error is a measure of the skill in predicting the evolution in time of the 
atmospheric systems. 

The error has been computed by subtracting the GLAS FGGE analysis frr.ii the 
forecast. Even though the choice of analysis clearly influences the "error," 
over most of the globe, the forecast error after one day is larger than the 
uncertainty in the analysis. A comparison made with the NMC operational analy- 
sis, based on a very different analysis/forecast scheme, and which used only 
NESS operational winds and no satellite temperatures in the Northern Hemisphere, 
agrees well with the results presented here in most regions. 

3.1 FGGE forecasts 

Figs. 1 to 3 correspond to the forecast error in the meridional velocity v 
as computed from FGGE initial conditions. Fig. 1 presents the average error in 
V at 850 mb after 1, 3, and 5 days. It may be seen that at low levels the 
systematic errors are dominated by large scales, both in the tropics and in the 
extratropics. This, combined by the fact that their phase is rather constant, 
indicates that they are associated with forcing, both thermal and orographic. 

For example, the fact that the forecasts overpredict the equatorward flow over 
the Andes even after one day, indicates that the mountains are generating more 
drag in the model than in the real atmosphere. At 300 mb, the average error in 
the tropics is still of planetary scale, but in the extratropics the error is 
of cyclonic scale. This, and the change in phase in the error after 1, 3 and 
5 days, indicates that the error in the extratropics is dominated by the 
systematic component in the forecast of moving cyclones. 

The systematic error grows in amplitude from day 1 to day 3. There is 
further growth from day 3 to day 5 in the extratropics, indicating further 
forecast skill. In the tropics the systematic errors seem to have attained 
their maximum amplitude by day 3. 
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Fiq. 3 presents the ratio between the average of systematic errors and the 
standard deviation of the error at 850 mb and 300 mb. Values smaller than one 
indicate that the error is dominated by transient features, and values larger 
than one indicate that the systematic error is more important. It may be 
observed that the extratropical error is dominated by the transients, whereas 
in the tropics, the systematic error is very important, especially at low levels. 
Fiq. 4 presents the heatinq rate at 500 mb as computed by the model during the 
January 1979 assimilation cycle. A comparison of Figs. 3 and 4 confirms that 
the large systematic errors are associated with regions of strong heating, as 
well as with orographic forcing. 

2,Z Comparisons of FGGE, NOSAT and PERSIS Forecasts 

Here we compare the systematic and transient errors of the forecasts 
obtained from the FGGE and NOSAT analysis cycles. Persistence forecasts, in 
which the forecasts coincide with FGGE initial conditions, are also presented 
and denoted PERSIS. Figs. 5 and 6 present the 3-day mean and standard devia- 
tion of the error in the zonal velocity u at 300 mb. 

It may be seen that the systematic errors from NOSAT initial conditions 
(Fiq. 5b) are only slightly larger than those of FGGE initial conditions. This 
indicates that systematic errors are due more to model parameterization defi- 
iencies than to initial data. It is interesting that both forecasts show 
characteristics similar to those of a "warm episode" of the Walker circulation, 
with enhanced easterlies and stronger subtropical jets in the Pacific (Horel 
and Wallace, 1981; Julian and Chervin, 1978). At low levels, not presented 
here, the error is reversed, completing an eastwest circulation. The system- 
atic error in the PERSIS forecast (Fiq. 5c) is much smaller than either NOSAT 
or FGGE in the tropics. This is not surprising because, for a large enough 
sample, there should be no svstematic errors in persistence forecasts. In the 
extratropics, the average PERSIS errors are dominated once again by the 
cyclonic scales that have the largest changes after 3 days. It should be 
remembered that the 14 forecasts are spaced by intervals of 4 days. 

Fiq. 6 presents the standard deviation of the 3-day forecast errors of u. 

The regions with errors smaller than 6 m sec"^ have been hatched, and those 
with errors larger than 12 m sec"^ are indicated by bold contours. A comparison 
of Figs. 6a and 6b indicates that the use of the FGGE special observing system 
has improved both the tropical and extratropical forecast of the transient 
features. We see in Fiq. 6c that the PERSIS forecast of transient features 
is, not surprisingly, much worse than either FGGE or NOSAT in the extratropics. 

The FGGE forecoSt errors are better than those of PERSIS in the subtropics 
and slightly better in the tropics, indicating some skill in predicting transient 
features. 

The improvenient in the forecast of the transient features using the FGGE 
analysis compared to either the NOSAT or PERSIS forecast is also clear after 5 
days, both in the tropics and in the extratropics. 
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4. CnNCLllSIONS 


From this prelininary study several conclusions may be drawn. We have 
found that the systematic error dominates the tropical forecast error. This 
error seems to be more dependent on model deficiencies than on the initial data 
and becomes large amplitude in a few days. The model forecast error then becomes 
comparable to the persistence forecast error in 3 to 5 days. On the other hand, 
the model retains some skill in the prediction of transient features in the 
tropics after 3 days and there is a significant impact of the FGGE data in the 
prediction of these transients. 

This study suggests that a major obstacle in accurate low latitude fore- 
casting is the prediction of the large scale quasi -stationary tropical circu- 
lation. 
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Figure 1. Mean forecast errors in the meridional velocity v at 850 mb 
Interval : 3 m sec'^. 
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Figure 4. January assimilation heating rates at 600 mb (deg/day) 
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Mean errors in the 3-day forecast of the zonal velocity at 300 mh 
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INTERACTIVE PROCESSING OF SEASAT SCATTEROMETER DATA 
R. Atlas, G. Cole, A. Pursch, and C. Long 


The objectives of the interactive processing of Seasat data are to (1) 
develop and apply techniques for the removal of alias in the SASS winds using 
all available data, (2) study the properties of SASS wind data in a variety of 
synoptic situations, and (3) provide improved analyses for meteorological diag- 
nosis and prediction. 

The Man-computer Interactive Data Access System (Me IDAS), developed by the 
Space Science Engineering Center (SSEC) of the University of Wisconsin and 
adapted for the interactive processing of satellite-derived temperature sound- 
ings and cloud-track winds for the FGGE Special Effort, has now been adapted 
for the processing and evaluation of Seasat data. The implementation of the 
McIDAS Seasat processing system required (1) extensive modifications to the 
data base structure to store and display SASS winds, as well as all corrobora- 
tive level II data, model first guess fields and level III analyses, and (2) the 
development of software to dealias and analyze SASS wind vectors interactively. 

The procedure for the subjective dealiasing of SASS winds that has been 
developed for use on the McIDAS consists of the following steps: 


1. Display SASS vectors. 

2. Superimpose auxiliary data. 

3. Remove alias where only two alias directions or ship wind data are 
present. 

4. Use pattern recognition, kinematic consistency, and meteorological 
experience to streamline entire SASS wind field. 

5. Remove alias of all remaining SASS vectors. 

6. Insure consistency with cloud imagery. 


This procedure has been applied to selected cases during the period 8-15 
September 1978. For situations in which nearly uniform wind directions exist 
over a large area, the synoptic analyst utilizing the McIDAS may dealias the 
entire wind field for that area by positioning a cursor over a single SASS wind 
vector which indicates the correct direction. For cases in which there is 
large variability of the wind directions, each SASS report is dealiased separ- 
ately using the same technique. 

Fig. 1 presents an example of the dealiasing of SASS winds on the McIDAS. 
In Fig. la, alias SASS wind directions and wind speeds are shown for a small 
section of the North Atlantic on 8 September 1978. The dealiased wind field 
for this section is shown in Fig. lb, where streamlines and ship report wind 
directions have been superimposed. Fig. 2 shows the dealiased wind field for 
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the entire North Atlantic from 0000-0036 GMT 8 September 1978, which incorpor- 
ated the section shown in Fiq. 1. The SASS winds are in good general agreement 
with most conventional ship reports and serve to define the synoptic and sub- 
synoptic scale flow with greater accuracy than is possible from the conventional 
data alone. Comparisons of the convergence and divergence fields associated 
with the SASS winds with geostationary cloud Imagery shows excellent agreement. 


a. 



SER SRT. .l^a_ND S 


i*c 
















330 



Fig. 1 Alias SASS vectors (a), and deallsed vectors and ship wind directions 
(b) for 0024-0030 GMT 8 September 1978. 
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DERLIRSED WINDS 



Fig. 2 Dealiased SASS winds, ship reported wind directions and surface 
streamlines for 0000-0036 GMT 8 September 1978. 
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PRELIMINARY EVALL' JION OF THE FGGE SPECIAL EFFORT FOR DATA ENHANCEMENT 
R. Atlas, A. Pursch, C. Long, G. Cole, and R. Rosenberg 


The FGGE Special Effort is being conducted for the two FGGE special observ- 
ing periods (January to March and May to June 1978) as a joint project between 
NASA, NOAA, and the University of Wisconsin. Man-computer Interactive Data 
Access System (McIDAS) terminals developed by the Space Science Engineering 
Center (SSEC) of the University of Wisconsin are being utilized by experienced 
meteorologists at the National Meteorological Center (NMC) for data evaluation 
and quality assessment and at the Goddard Laboratory for Atmospheric Sciences 
(GLAS) for satellite data enhancement. The interactive procedures used for 
sounding and cloud-motion wind enhancement were originally developed by the 
National Environmental Satellite Service and SSEC. (Smith et al . , 1978, and 
Menzel et al . , 1978). (See Greaves et al . , 1979 for a detailed description of 
the devel opment and implementation of the FGGE Special Effort.) 

For each synoptic period there is an initial editing of FGGE data by NMC. 

At this stage, quality indicators are assigned to temperature soundings from 
TIROS-N and to cloud-tracked wind vectors from geostationary satellites. The 
determination of data quality is made on the basis of synoptic considerations 
including horizontal, vertical and temporal consistency. Regions where data 
deficiencies exist or where higher resolution data is needed to adequately 
represent the atmospheric thermal structure are then selected for enhancement. 
The eastern North Pacific is routinely enhanced. In addition, cases have been 
preselected on the basis of meteorological interest, by an ad-hoc international 
committee of participating scientists. These cases include situations of 
blocking, cut-off low development, cycloqenesis, and tropical circulation. 

The sounding data enhancement is aimed at supplementing the operational 
satellite sounding data set with higher resolution soundings in meteorological 
active regions and with new soundings where data voids or soundings of question- 
able quality exist. The algorithms for retrieving temperature profiles from 
the TIROS-N observations of radiance displayed on McIDAS are essentially the 
same as those used by NESS for the objective generation of operational tempera- 
ture profiles. Three ty.es of temperature retrievals are possible: (1) dear- 

column, and (2) partly cloudy retrievals utilizing infrared observations from 
the HIRS instrument aboard TIROS-N, and (3) cloudy retrievals utilizing only 
microwave observations from the Microwave Sounding Unit (MSU) on TIROS-N. The 
operational temperature profiles have a horizontal resolution of 250 km, whereas 
special effort temperature profiles can be retrieved at the resolution of the 
fueasurements (30 km for HIRS and 150 km for MSU). 

The sounding enhancement process at GLAS is performed as follows: once 

an area has been selected, the McIDAS operator begins by displaying all avail- 
able conventional and special FGGE data for that area and visible, infrared, and 
microwave images from TIROS-N. The operator notes where data deficiences exist 
and, from the TIROS-N images, determines the extent of cloudiness and where 
the most intense atmospheric thermal gradients are located. High-resolution 
infrared temperature retrievals are then generated for the area. This is fol- 
lowed by a comprehensive manual editing of the retrievals to remove small-scale 
discontinuities due to cloud-induced noise while retaining significant meteoro- 
logical structures. Enhanced microwave retrievals, consistent with neighboring 
infrared retrievals, are then generated in cloudy areas. 
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After the entire enhancement process Is completed, a final editing and 
quality assessment of the enhancement data is performed at NMC. The data is 
then archived as supplementary Level Ilb and will be utilized by N0AA‘s Geo- 
physical Fluid Dynamics Laboratory (6F0L) in generating Level Illb analyses. 

Enhanced sounding data can be displayed on the McIDAS as either mandatory 
level temperatures or as the thickness between any two levels. Analyses of 
these quantities as well as thermal winds, vertical temperature profiles, and 
vertical cross-sections of potential temperature are routinely produced. If 
conventional surface data is available, then the 1000 mb geopotential height 
is analyzed, and the satellite-derived thickness and thermal winds are added 
to the 1000 mb height and 1000 mb qeostrophic wind to yield geopotential height 
and qeostrophic winds on constant pressure surfaces. In addition to the high 
resolution analyses performed directly on the McIOAS, enhanced soundings may 
be assimilated into larger scale objective analyses which then serve as Initial 
conditions for numerical model predictions. 

Subjective comparisons of enhanced soundings with all available operational 
soundings and radiosonde reports are performed at the end of the enhancement 
process. These comparisons have shown that the enhanced and operational retrie- 
vals tend to be similar in cloud-free areas. However, large differences occa- 
sionally occur. The differences are accentuated in and around cloudy areas 
where the interactive processing at high resolution allows for the generation 
of substantially more Infrared retrievals and where the microwave retrievals 
are enhanced. 

Examples illustrating some of the effects of the sounding data enhancement 
are presented in Figures 1 and 2. In Figure 1, 1000 mb to 300 mb thickness 
analyses of operational (solid lines) and Special Effort retrievals (dashed 
lines) and plotted radiosonde reports are shown for a section of the North 
Pacific at 0100 GMT 7 January 1979. Large differences in the orientation, 
gradient, and absolute value of the enhanced and operational thickness contours 
are evident. In addition, the errors relative to colocated radiosondes have 
been reduced by 30 to 80 geopotential meters. 

In Figure 2, vertical temperature profiles of unenhanced and enhanced 
microwave retrievals are plotted relative to the colocated radiosonde observa- 
tion at Yakutat, Alaska at 0000 GMT 11 January 1979. Differences between the 
retrievals are small throughout most of the troposphere in this case. However, 
the enhancement process has resulted in a better definition of the tropopause 
and a reduction of error of 4*C near the 100 mb level. 

As an example of the application of Special Effort sounding to large-scale 
meteorological analysis and prediction, the enhancement for 0000 GMT 9 January 
1979 is presented here. This was a case in which a large gap in the operational 
TIROS-N soundings occurred in the eastern North Pacific. Figure 3 shows the 
orbit of TIROS-N and the location of operational soundings from 2100 GMT 8 
January to 0300 GMT 9 January. Reprocessing of observed radiances on the McIDAS 
allowed for the generation of enhanced soundings over a large portion of the 
data gap from 30* to 50*N. Additional high resolution soundings were also 
generated over regions adjacent to the gap to better define the atmospheric 
thermal structure and to insure consistency. 
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Figure 2 . Vertical temperature profiles derived from unenhanced 
microwave (M), enhanced microwave (E), and radiosonde (R) for 
Yakutat, Alaska, at 0000 GMT 11 Jan 79. 



Figure 3. Location of TIROS-N retrievals from 2100 GMT 8 Jan 79 - 0300 
9 Jan 79. 


Figures 4 and 5 display the 1000 to 300-mb thickness analyses and ther- 
mal winds for the enhanced regions. These analyses are based completely on 
the high-resolution soundings, with each thermal wind barb plotted at the 
location of the retrieved vertical temperature profiles. 
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Figures 6 and 7 illustrate the effect of the enhanced soundings on the 
GLAS ohiective analysis for 0000 GMT 9 January. In Figure 6, the 300 mb height 
and wind analysis for all conventional and operational satellite data is pre- 
sented. The corresponding analysis, which includes enhanced soundings, is 
presented in Figure 7. Comparison of the figures reveals that a major modifi- 
cation to 300 mb heights and winds has occurred just off the west coast of the 
United States. The enhanced soundings have increased the amplitude of the 
ridge at 135“W and shortened the half-wavelength between this ridge and the 
trough along the west coast. Elsewhere the effects are smaller due to the 
consistency of the operational and enhanced soundings over most of the eastern 
North Pacific and the averaging of the effects of enhanced soundings by the 
objective analysis. However, changes to the intensity of the geopotential 
gradient and orientation of the winds and contours are evident throughout this 
region. 

Figures 8 and 9 show the 24-hour GLAS model forecasts of 300 mb geopoten- 
tial heights and winds, which were generated from initial conditions corres- 
ponding to Figures 6 and 7, respectively. The verifying analysis is presented 
in Figure 10. Comparison of these figures reveals that the inclusion of 
enhanced soundings resulted in a large improvement in the predicted amplitude 
of the trough over the southwestern United States. Additional improvements 
are evident over the North Pacific. However, the verifying analysis is less 
certain in this area. 



Figure 4. 1000-300 mb thickness and thermal wind at 2300 GMT 8 Jan 79 
derived from high resolution TIROS-N soundings for the region 30-55“N 
and 12n-145“W. 
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Fiqure 6. 300 mb qeopotential height and wind analysis without en 

hanced retrievals for 0000 GMT 9 Jan 79. 








Figure 8. 24h 300 mb forecast generated without enhanced retrievals 


Figure 9. 24h 300 mb forecast generated with enhanced retrievals 


Figure 10. Verifying 300 mb geopotential height and wind analysis 
for 0000 GMT 10 Jan 79. 


Subjective evaluation of the effects of sounding data enhancement Indicate that: 


1. Enhanced and operational IR retrievals are similar In most situations. 

2. Large Improvements In 1000 to 300 mb thickness and mandatory level 
temperature and Intensification of atmospheric thermal gradients 
occasionally occur. 

3. Interactive processing can fill In major gaps In the TIROS-N coverage. 

4. Assimilation of Special Effort soundings can result In major modifica- 
tions to large-scale analyses and prognoses. 
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COMPARISON OF OBSERVED AND COMPUTED BRIGHTNESS TEMPERATURES 
FOR THE HIRS2 AND MSU SOUNDERS ON TIROS-N 

J. Susskind, J. Rosenfleld, and D. Reuter 


Global temperature soundings from HIRS and MSU data on the TIROS-N satel- 
lite are produced at the Goddard Space Flight Center Laboratory for Atmospheric 
Sciences Modeling and Simulation Facility using a direct physical retrieval 
method. Intrinsic in the success of such a method is accurate computation of 
expected radiance observations » given a set of atmospheric and surface conditions. 
McClatchey (1976) indicated very poor agreement between radiances observed by 
the DMSP sounder under clear conditions and those calculated from colocated 
temperature humidity profiles, raising questions about the accuracy in which 
radiative transfer calculations can be performed in the infrared. Here, we 
briefly summarize our findings, showing excellent agreement between observed and 
computed brightness temperatures for HIRS2 and MSU under both clear and partly 
cloudy conditiors. More details are given in Susskind et al . , (1982). Table 1 
shows the channels, centers, and peaks of the weighting functions, or other 
relevant information, for the channels on MSU and HIRS2. The HIRS2 and MSU 
channels used in the GLAS physical inversion method are designated by 

Table 1 

HIRS2 and MSU Channels 


Channel 

v(cm“l ) 

Peak of 
dt/dlnp (mb) 

Peak of 
Bdr/dlnp (mb) 

HI 

668.40 

30 

20 

*H2 

679.20 

60 

50 

H3 

691.10 

100 

100 

*H4 

703.60 

280 

360 

H5 

716.10 

475 

57-. 

H6 

732.40 

725 

875 

H7 

748.30 

Surface 

Surface 

H8 

897.70 

Window, sensitive to 

water vapor 

H9 

1027.90 

Window, sensitive to 


HID 

1217.10 

Lower tropospheric water vapor 

Hll 

1363.70 

Middle tropospheric water vapor 

H12 

1484.40 

Upper tropospheric water vapor 

*H13 

2190.40 

Surface 

Surface 

*H14 

2212.60 

650 

Surface 

*H15 

2240.10 

340 

675 

H16 

2276.30 

170 

425 

H17 

2310.70 

15 

2 

*H18 

2512.00 

Window, sensitive to 

solar radiation 

*H19 

2671.80 

Window, sensitive to 

solar radiation 

♦Ml 

50.30a 

Window, sensitive to 

surface emissivity 

♦M2 

53.74a 

500 


♦M3 

54.96a 

300 


♦M4 

57.95a 

70 



a values in 6^2 
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All HIRS2 observations utilized In global analysis of the data are composed 
of two averages of spots within a 125 x 125 km area. The two fields of view 
represent relatively clear and relatively cloudy areas. Both sets of radiances 
are used to determine a reconstructed dear-column radiance* that Is, the radi- 
ance which would be observed If the areas were cloud-free (Susskind and Rosenfleld* 
1980). For verification purposes* HIRS2 observations during the period Jan 5 
1979-Feb 25 1979 were chosen which were colocated with a select set of oceanic 
radiosondes (see Phillips et al . * 1979) In space by 110 km* and* In time, by three 
hours. Only those radiosondes reporting temperatures at all mandatory levels 
from 1000 mb - 30 mb and humidities from 1000 mb - 500 mb were used In the 
comparison. 


Olven the radiosonde temperature and humidity profile* brightness tempera- 
tures for the temperature sounding IR channels on HIRS2 were calculated accord- 
ing to Susskind et al . * (1982). The temperature profile was interpolated to 
53 levels between 1000 and 30 mb assuming temperature to be linear In the logar- 
ithm of the pressure between the mandatory levels. The 11 levels above 30 mb 
were extrapolated according to climatology. The surface temperature was taken 
to be the climatological sea surface temperature. The specific humidity was 
interpolated between mandatory levels assuming a dependence. A specific 
humidity corresponding to a climatological water vapor mixing ratio of 2 ppmv 
was assumed at and above 100 mb. Zonal ly averaged climatological ozone profiles 
were used to compute the ozone component of the transmittance. 


Infrared radiances are potentially cloud contaminated. An estimate of the 
dear-column radiance* R^, for channel 1 can be reconstructed from the observa- 
tions In adjacent fields of view according to 


Ri = Ri,i + n (Ri,i - Ri,2) 


( 1 ) 


where R^^j is the observation for channel i In the field of view j and n Is a 
parameter given by al/(a2 - al)* where al and a2 are the cloud fractions in the 
two fields of view (Susskind and Rosenfield, 1980). For the purposes of this 
study* the parameter n* which Is independent of channel, is best determined 
according to 


n 



13*1 


R13*1 - R13*2 (2) 

cc 

where Rj3 is the computed dear-column radiance for channel 13, using the radio- 
sonde temperature profile. 

A value of n = 0 corresponds to field of view 1 (always taken as the field 
of view containing the larger radiance In the 11pm window channel) being clear. 
If n was less than 0, n was set equal to -0,5. In this case, the dear-column 
radiance Is the average of the, observations in both fields of view. 

Two sets of statistics were calculated comparing observed and computed 
radiances. The first sample represents comparlslons made only in those cases 
where at least one field of view was found to be clear* that Is n<0. The second 
sample contains only partially cloudy cases* n>0* in which case compa* isons were 
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made between radiances reconstructed from equations 1 and 2 and computed clear- 
coluiwi radiances. The first set of statistics refers primarily to the accuracy 
of the forward problem calculation. The second set Includes the accuracy of 
the ability to account for cloud effects. 

Table II shows the comparison of observed and computed brightness tempera- 
tures for the air and surface temperature sounding channels of HIRS2 and MSU in 
the 7? cases where the cloud filtering algorithm showed at least one field of 
view to be cloud-free. Those channels used in the retrieval of temperature 
profiles are marked by *. Columns 2 and 3 show the mean and standard deviation 
of the differences between computed and observed brightness temperatures. The 
standard deviation between observed and computed brightness temperatures is 
generally better th^n 1°K for the HIRS2 channels used in the calculation and 
somewhat worse than 1®K for the MSU channels. Significant bias errors are 
found in some channels, however. To a large extent, these biases are due to 
errors in the line oarameters used in the line-by-line calculations. The 
largest standard deviations found in the infrared channels are in channels 1-3 
and 16, whose radiances are quite sensitive to the temperature profile above 
30 mb, and channels 8, 18, and 19, the surface temperature sounding channels, 
where agreement is limited by the accuracy of the climatological sea surface 
temperature. Channel 16 is also affected by non-local thermodynamic equilibrum, 
which was not accounted for in the calculation. Channel 17, which is affected 
even more greatly by non-local thermodynamic equilibrium, was not included in 
the calculation. Observations in channels 18 and 19 are affected by contribu- 
tions of reflected solar radiation to the radiances. In addition, agreement in 
all channels is degraded somewhat by sampling differences between radiosonde 
and satellite observations. 

The systematic errors can be removed by empirically tuning the transmit- 
tance functions in the form 

MC) * (3) 

which is equivalent to multiplying the effective channel optical depths by 1+e^. 

Tuning coefficients were found for channels, 4-7, 13-15, and M2-M4 which 
minimized the bias and standard deviations in the brightness temperature errors 
r-'-^puted in the sample of clear cases. The tuning coefficients obtained for 
inese channels are shown in column 6 of the table and the resulting mean and 
standard deviation of the errors of the brightness temperature using the tuned 
transmittances are shown in columns 4 and 5. No attempt was made to tune the 
remaining channels because of the considerations discussed above. The empirical 
tuning coefficients found in Table II are similar to those found by Weinreb 
(1979). 

After tuning, the mean errors for all channels except 1, 3, and 16, which 
are not used in analysis of the HIRS data, are insignificant. With the exception 
of the stratospheric sounding channel 2, the standard deviation of the 15pm and 
4.3um brightness temperatures in channels used for retrievals of temperature 
profiles are of the order of 0.7®. These errors are sufficiently small so as 
not to seriously effect the accuracy of retrievals. Agreement in channel 7, 
used only for cloud height determination, is somewhat worse, possibly because 
of the effects of incorrect surface temperature. Likewise, agreement in the 
microwave channels is about 50X worse than the infrared channels, possibly 
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because of higher noise in the microwave measurements. The sixth column of 
Table II shows the standard deviations of the measurements for each channel. 

The ratio of the standard deviation of the observations to the standard devia- 
tion of the errors may be thought of as a signal to noise for each channel. 

Agreement in channels 18 and 19 appears to be very poor, with a large 
negative bias in the observations. This is due to the effects of solar radia- 
tion, which were not included in the calculations of brightness temperatures 
for these channels. In the physical inversion of sounding data, the two obser- 
vations for channels 18 and 19 are analyzed simultaneously to compute a surface 
temperature and a correction for solar radiation (Susskind and Rosenfield, 

1980; Chahine, 1980). That procedure was not done here because it would force 
agreement between observed and computed brightness temperatures. To assess the 
accuracy of the forward calculation for these channels, statistics were gathered 
for observations at night (21 cases). These are shown in the last two columns 
of Table I for the window channels 18, 19, and 8. The general agreement is 
good and again limited by the accuracies of the climatological sea surface 
temperature. The two shortwave windows are shown to be somewhat more accurate 
than the 11 wm window channel. This is most likely due to the lesser effect 
of absorption due to water vapor continuum on the radiances in the 3.7 ym 
region. Absorption due to water vapor continuum is not thoroughly understood 
and less well parameterized than the other sources of atmospheric absorption. 

Agreement of reconstructed and computed brightness temperatures in an 
independent sample of 149 partially cloudy cases is shown in Table III. The 
tuning coefficients shown in Table II, which were determined from the sample 
of clear cases, are used in the calculations. The third column shows the 
root mean square difference of the reconstructed clear column brightness 
temperatures and the observations in the warmer field of view. In channels 
1-3, sounding the stratosphere, the reconstructed clear column radiances are 
taken as the average of the radiances in the two fields of view, both of which 
are assumed to be unaffected by clouds. In these channels, the difference of 
observed and reconstructed brightness temperature is a measure of the instru- 
mental noise. In the other channels, it represents the correction made to the 
brightness temperature in accounting for cloud effects. 

The mean and standard deviations of the errors of computed compared to 
reconstructed dear-column brightness temperatures are shown in columns 2 and 3. 
The bias errors are small, indicating the applicability of tuning coefficients 
determined from a set of clear observations for use on a set containing differ- 
ent, partially cloudy, observations. The general agreement in those channels 
used in the analysis is comparable to, though slightly worse, than that in the 
clear cases. The lower tropospheric 15 ym channels 5-7 are degraded somewhat 
under partially cloudy conditions, but are not used in the analysis. Agreement 
in the window channels is shown at right in the last two columns and show 
comparable accuracy to that in clear cases. The "errors" in channel 13 are 
zero because this channel is used for the cloud correction. 

The results of this study indicate that infrared radiative transfer cal- 
culations can be done with the same, if not better, accuracy than microwave 
radiative transfer calculations in cloud-free and partially cloudy conditions. 

The agreement of observed brightness temperature with those computed using 
colocated radiosonde information is sufficient to allow for accurate physical 
retrievals from the satellite data in both clear and partially cloudy conditions. 
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TABLE II 


COMPARISON OF OBSERVED AND COMPUTED BRIGHTNESS TEMPERATURES C*K), CLEAR CONDITIONS 















TABLE III 


COMPARISON OF RECONSTRUCTED AND COMPUTED BRIGHTNESS TEMPERATURES. CLOUDY CONDITIONS 


DAY AND NIGHT (140) NIGHT (64) 


TUNED 

MEAN 

(COMP-RECON) 

STD DEV 

RMS 

(RECON. 

-OBS.) 

VAR. OF 
RECON. 

MEAN 

(COMP-RECON) 

1 

2.72 

3.98 

1.34 

3.97 


* 2 

-.06 

1.65 

.38 

5.14 


3 

1.25 

1.47 

.29 

5.02 


* 4 

.04 

1.13 

1.10 

3.19 


5 

.08 

1.34 

2.52 

4.83 


6 

-.00 

1.61 

3.82 

5.91 


7 

.27 

1.82 

5.06 

5.84 


ft 

.39 

1.44 

7.31 

7.21 

.35 

*13 

-.00 

.00 

3.87 

8.11 


*14 

-.09 

.56 

3.14 

8.23 


*15 

.03 

.96 

2.26 

6.86 


16 

.64 

2.19 

.64 

3.32 


*18 

-1.05 

2.74 

4.14 

9.64 

.37 

*19 

-4.38 

5.19 

3.71 

9.54 

-.54 

Ml 

0 

0 


10.69 


M2 

.03 

.98 


7.61 


M3 

-.15 

1.06 


4.52 


M4 

-.52 

1.51 


8.55 
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SOLUTION OF LINEAR SYSTEM OF EQUATIONS 
WHEN THE MATRIX IS CIRCULANT 


\ R. C. Balqovind 
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The objective of this note Is to make "could be users" aware of the exis- 
tence of such a methodoloqy. This type of tool has been In use since the Intro- 
duction of the Fast Fourier transform (see all references except the following). 
Here we Introduce the general circulant matrix and show how to "Invert" It. 

Evans (19B0) has shown that the Quin-diagonal system can also be solved by L-U 
factorization. 


1. PROBLEM AND THE METHOD OF SOLUTION 
Our Interest Is to solve the system 

Ax = b 

Where A Is the non-singular matrix of the form 


Co 

Cn-1 

Cn-2 

Cn-3 Cn-4**«C3 C2 

Cl 

Cl 

Co 

Cn-1 

Cn-2 Cn-3«“C4 C3 

C2 

C2 

Cl 

Co 

Cn-1 Co-2*»»C5 C4 

C3 


^n-3 ^n-4 ^n-5 ^n-6 ^n-7‘» 
^n-2 ^n-3 ^n-4 ^n-5 ^-6»* 
^n-1 ^n-2 ^n-3 ^n-4 ^n-5** 

One can express the above matrix in the form; 


^0 ^n-1 ^n-2 
^1 ^0 ^n-1 
Cz Cl Co 


A = f(P) = Cq 

where 


P Is the permutation matrix 


I + P + C2 P^ + 

r 

j 0000... 
1 000... 
0100... 
001 0 ... 
0001 ... 


0000... 


0 1 
0 0 
0 0 
0 0 
0 0 

1 0 
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I would like to recall here that 

Pe-| “ 61+1 » i“l »2,3»***»n-l > 

and Pep “ e] where ei are the usual unit vectors. 

Let w * exp(at1/n) and » l.//n (l.w^, w^'', Then A 

has eigenvalues x^« f(w”^) for r » 0,1 ,2,... ,n-l corresponding to eigenvectors 
Vp, since P'^p “ w‘‘'‘vp. 

Our algorithm relies on the fact that the eigenvalues Xp and eigenvectors 
Vp of P are known. Let Q be the matrix whose columns are^Vg, , V2»**»»'^n-1 
and let A = diag(XA, Xi, Xof.*! ^n-O* Q Q * I. Q AQ “ A and it follows 

from equation (1) that 


X * Q( A-l(Q*b)) . 

Since these correspond to a direct and an inverse Fourier transform, they 
can be performed efficiently by Fast Fourier transform algorithm. 


2. EXAMPLE 


Let us look at an example that arises when 

{ c + ii_ ) u = RHS 
3xSx 

(for u periodic and c constant) is approximated with a 9-point fourth order 
scheme. This gives us a matrix of the form: 


A 


Co Cl Cz C3 C4 0 
Cl Co Cl Cz C3 C4 
C2 Cl Co Cl C2 C3 

Cl C2 Cl Co Cl C2 
C4 C3 C2 Cl Cq Cl 
0 C4 C3 C2 Cl Co 

o*‘o**o’*o**o**o 

0 0 0 0 0 0 

C4 0 0 0 0 0 

C3 C4 0 0 0 0 

C2 C3 C4 0 0 0 

Cl C2 C3 C4 0 0 


0 

0 

C4 

C3 

C2 

Cl 

0 

0 

0 

C4 

C3 

C2 

0 

0 

0 

0 

C4 

C3 

0 

0 

0 

0 

0 

C4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Co 

cl 

C2 

C3 

ci' 

o' 

Cl 

Co 

Cl 

C2 

C3 

C4 

C2 

Cl 

Co 

Cl 

C2 

C3 

C3 

C2 

Cl 

Co 

Cl 

C2 

C4 

C3 

C2 

Cl 

Co 

Cl 

0 

C4 

C3 

Cz 

Cl 

Co 
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In this special case (where matrix A Is symmetric) we have that 

f(x) » Cq + X + C^x^ + C3X^ + C4X^ 

+ x"”^ + + C^x''’^ + C4x”‘^ . 

So the characteristics corresponding to the eigenvectors Vp are 
Xp = f(w“^) 

= Co + 2Ci cos(r (2if/n)) + 2 C 2 COs( 2 r( 2 ir/n)) 

+ 2 C 3 cos(3r(2x/n)) + 2 C 4 COs( 4 r( 2 ir/n)) 
for r = 0 , 1 , 2 ,...,n-l . 

This scheme, combined with block reduction method, may be used to solve 2«d 
problems. 
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AN ANALYSIS OF INTERMITTENT UPDATING 
FOR THE SHALLOW WATER EQUATIONS 

K. P. Bube 


1. INTRODUCTION 

Numerical weather prediction (NWP) is the attempt to solve an initial- 
value problem for a system of nonlinear partial differential equations. One 
of the major difficulties in NWP is the specification of initial conditions; 
that Is, the initial data are known only incompletely and inaccurately. For- 
tunately, o! ervations of the system distributed over a time interval preceding 
initial time are also available. The problem of four-dimensional data assimi- 
lation is to construct a complete, accurate set of initial data for a NWP model 
from all of the measured data. 

The most common practical procedure that used data available at times pre- 
ceding the initial time is intermittent updating. The model is provided with 
the best available data at some preceding time, and is integrated forward in 
time. Additional data replace the model values as they bicome available. When 
the integration reaches the initial time, an inital state has been constructed. 
Many variations of this procedure have been proposed and studied (Bengtsson 
(1975); Ghil, !!alem, and Atlas (1979)). Recently, the application of the esti- 
mation theory to improve the procedure has been studied (Ghil, et al. (1981)). 

The purpose of this work is to analyze the method of intermittent updating 
for a model system. We want to understand analytically how the updating of one 
variable of the system affects the other variables, what the factors are which 
infljence the rate of convergence of this procedure, and how measurement error 
affects the procedure. Preliminary results appear in Bube and Oliger (1977), 
Bube (1978), and Bube and Ghil (1981). Talagrand (1977, 1981) has also studied 
this procedure analytically. 


2. UPDATING WITHOUT ERRORS 

We present here results for the model system 


a b Uj^, 
Ut = b + a Ujf, 


( 1 ) 


with periodic boundary conditions in x. The linearized shallow water equations 
for purely zonal, one-dimenisonal flow around a circle of latitude can be 
written in this form, where «j> is the geopotential, u is the zonal velocity, and 
-a = U and -b = ♦ are the zonal velocity and oeopotential of the basic state. 
We assume that we have measurements of «>(x,t) ”or all x and for t = jr, j = 

0,1, ^....,m, where t > 0 is the time interval between observations of <|>. 

We start with the initial measurement of ^ and an initial approximation to 
u at t = C. As we integrate forward in time, we replace the computed values of 
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♦ by the measured values at times t = .It. The error in the computed u, which 
which we call e(x,t), will decrease because of the linkage between u and ♦ in 
the system. 


For the analysis, we expand u, and e in Fourier series in x; 


1 

♦ (C.t) = e-2niCx 

♦(x,t) dx, 

0 

♦ (x,t) = I <KC.t) e2”i^x. 


( 2 ) 

( 3 ) 


where ^ is the wave number (an Integer). We can show that 

K(C,(.1+l)x)| = o(^,T) le(C,jx)|. (4) 

where 

p(C,x) = lcos(2"b5x)l. (5) 

The effect of each update interval is to multiply each Fourier coefficient of 
the error in u by a decrease factor p(C,x) i. If p(C,x) « 1, then «^(^,jx) 
approaches zero slowly as j increases. 

For t = jx , 




( 6 ) 


where 


o(f;,x) = P(C,x)V\ 


( 7 ) 


For ^ 0, o(C,x) increases as x ♦ 0, and 


lim p(5,x) = 1 (8) 

x+0 . 


Updating ♦ more frequently (i.e. decreasing x) does not necessarily nidke 
e(x,t) 0 faster as t incr*»ases. We must allow enough time for the energy of 
the error to pass from u to and then out of the system as ^ is updated. 

Note that p and c depend on both the wave number ^ and the updating inter- 
val X. For C = 0, P = 0 = 1 ; no improvement is made in the mean velocity u{0,t). 
For other valuer of C, the best updating interval is 

1 . (9) 

X - 


If the entire error in u is in only one wave number then x can be 
chosen by (9) to decrease the error quickly. In practice, e(^,0) / 0 for many 


75 



values of If t » l/(4|hc|). then pU.t) is just as likely to be close ';o 
1 as it is to be close to 0. The best t for a small wave number may give a 
p -1 for some large wave numbers. If t is too large, waves with high wave num- 
bers may converge slowly; if t is too small, waves with low wave numbers will 
converge slowly. 

One possible solution is to use several different updating intervals 
Ti.....Tp, and to repeat them in a cycle. If there is a xj satisfying 
(9) for each wave number 5 of interest, then convergence should be reasonably 
rapid for all corresponding wave numbers. Numerical experiments with system 
(1) confirm these theoretical results. 

A similar analysis can be performed on a 3 x 3 version of the linearized 
shallow water equations, with Coriolis parameter f, and northerly wind component v 


^t u (>x + 4(ux+Vy) = 0. 
ut + U Ux + ♦x “ = 0, 
vt + U Vx + fu = 0. 


In this case, the decrease factor (depending on c and x). when <)> alone is meas- 
ured, is a 2 X 2 matrix. When a single updating interval x is used, the magnitude 
of the largest eigenvalue of this matrix determines the rate of convergence for 
this c and X (see Talagrand (1981)). If different x's are used in succession, 
the situation becomes more complicated and deserves further investigation. 


3. UPDATING WITH ERRORS 

Consider the situation discussed in the previous section using system (1), 
but now suppose that the measurements of ♦(x,jx) for j = 0,1..., m have errors. 
Suppose that the expected value of the square of the error in (C.jt) is 
r(^), for j = 0,1, ...,m. Let q(C.t) be the expected value of |e(^,t)|2. Then 
we can show 


q(C.nx) = r(^) + p(^) + p(C.x)2n (q(f,,0) - r(^) . 


Thus, the expected value of the error in the computed u decreases down to the 
expected value of the error in the measurements of <>, in a manner which is 
analogous to that of the error-free case. A corresponding result holds when 
a cycle of different updating intervals is used. 

Further analysis incorporating systematic error into these procedures is 
planned; the first model to be considered is the discrepancy between the solu- 
tions of a differential equation and sol u tens of corresponding difference 
schemes. 
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ASSIMILATION OF ASYNOPTIC DATA AND THE INITIALIZATION PROBLEMI 

K. P. Bube and N. Ghll^ 


We discuss a mathematical framework for the use of asvnoptlc data In deter- 
mininq initial states for numerical weather prediction (NWP) models (bengtsson 
et al . , 1981). A set of measured data, synoptic and asynoptic, is termed com- 
plete if it determines the solution of an NWP model uniquely. We derive theor- 
etical criteria for the completeness of data sets. The practical construction 
of the solution from a complete data set by intermittent updating Is analyzed, 
and the rate of convergence of some updating procedures Is given (Bube and Ghil, 
1981). 


It is shown that the time history of the mass field constitutes a complete 
data set for the shallow-water equations. Given that the time derivatives of 
the mass field are small at Initial time, we prove that the velocity field 
obtained by the diagnostic equations we derive will also have small time deriva- 
tives. Hence our diagnostic equations also solve the Initialization problem 
for this system; namely, they provide an initial state which leads to a slowly 
evolving solution to the system (Bube and Ghil, 1981). 

Finally, we review the bounded derivative principle of Krelss. It states 
that in systems with a fast and a slow time 'scale,' initial data can be chosen so 
that the solution starts out slowly. For such Initial data, the solution will 
actually stay slow for a length of time comparable to the slow time scale. The 
application of the principle to the Initialization problem of NWP Is discussed. 
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APPLICATIONS OF ESTIMATION THEORY TO NUMERICAL WEATHER PREDICTION 

Numerical weather prediction (NWP) is an initial^alue problem for a system 
of nonlinear partial differential equations in which the initial values are 
known only incompletely and inaccurately* Data at initial time can be supple- 
mented, however, by observations of the system distributed over a time interval 
precedinq it. Estimation theory bas^S^S^ successful in approaching such problems 
for models governed by systems of ordinary differential equations and of linear 
POEs. Wa diwe4np ^sti mat ion-theoretic methods for NWP^t^e. 

A model exhibiting many features of large/icale atmospheric flow Important 
in NWP is the one governed by the shal low^fluid equations. Wffsyaad y jhe ^ i.- / y, 

mation problem for a linearized formulation of these equations^''"TTftrfT ^^ sTudtf 
difference version of the equations is used as a forecast model to simulate the 
numerical models used in NWP. 

At each observation time, observed-minus-foreca't residuals are multiplied 
by a matrix of weighting coefficients, then added to the forecast itself to 
produce the "analyzed" field from which the forecast then proceeds. The standard 
estimation-theoretic approach is to use the so-called Kalman-Bucy (K-B) filter 
as the matrix of weighting coefficients. 

The specific character of the equations' dynamics gives rise to the neces- 
sity of modifying the usual K-B filter. The goal of the modification is to 
eliminate the high-frequency inertia-gravity waves which are otherwise generated 
by the insertion of observational data (Bengtsson et a1 . . 1981). Our numeri- 
cal experiments show that the modified filter accomplishes this, at the expense 
of errors which are only slightly larger than those produced by the optimal, 

K-B filter (Ghil et al . , 1981), The modified filter, therefore, combines in an 
optimal way dynamic initialization (i.e., elimination of fast waves) and four- 
dimensional (space-time) assimilation of observational data, two procedures 
which traditionally have been carried out separately in NWP, 

The modified filter converges rapidly to an asymptotic, constant matrix, 
which can then be used as a time-invariant filter. This filter is much more 
efficiently computable than its time-varying counterpart, and our numerical 
results show that these two filters perform equally well. 

Finally, we compare uur methods with the existing method in widest opera- 
tional use, namely, optimal interpolation (01) (Cohn et al. , 1981). Our results 
indicate that the asymptotic tilter yields a substantial improvement over 01. 

We expect that a suitable nonlinear version of our asymptotic filter, or a 
compromise between it and the one used in 01, can be implemented in an opera- 
tional NWP model. 


S. Cohn, M. Ghil^ and E. Isaacson 


^ Also Laboratory for Atmospheric Sciences, NASA Goddard Space Flight Center, 
Greenbelt, Maryland 20771. 

Courant Institute of Mathematical Sciences 
New York University, New York, New York 10012 


79 



REFERENCES 


Benqtsson, L,» M. Ghil, and E. Kallen, 1981: Dynamic Meteorol ogy : Data Assimi - 

lation Methods* Springer-Verlag, New York/Hei del berg/Berl in , to be published. 

Cohn, S., M. Ghil, and E. Isaacson, 1981: Optimal interpolation and the Kalman 

filter. Preprint , 5th Conf. Numer. Weather Prediction , Monterey, California. 

Ghil, M. , S. Cohn, J. Tavantis, K. Bube, and E. Isaacson, 1981: Applications 

of estimation theory to numerical weather prediction. In Bengtsson et al. 
(1981). 


80 



’‘^ 82 ' 25568 


A STOCHASTIC-DYNAMIC MODEL FOR GLOBAL ATMOSPHERIC 
MASS-FIELD STATISTICS 

M. Gh1ll, R. Bdlgovlnd^ and E. Kalnay^ 


Global atmospheric mass-field error correlations based on satellite obser- 
vations and on numerical forecasts show strong and systematic latitude dependence 
A model for the latitude-dependent spatial correlation structure of mass-field 
forecast errors is derived from dynamical considerations. The model is governed 
by the potential vorticity equation forced by random noise; 

(v^-CgSin^O) = F(x,(j»;w) (1) 


where v2 is the Laplacian operator on the unit sphere, x and 0 are longitude 
and latitude, ()> is the geopotential error field at 500 mb and F is white noise 
corresponding to a random realization u. 

The spatial covariance function r is defined by 

r(Xi ,O;X 2 , 02 ) * E [(J» (Xi,0i;u) (X2f02>)l» ,o\ 


where E is the expected value. 

Three methods of solution have been tested. In the fir^t method. Eg. (1) 
was solved by expansion in spherical harmonics, and the cor* elation function was 
computed analytically using the expansion coefficients. In the second method, 
the finite-difference equivalent of Eq. (1) was solved using a Fast Poisson 
Solver. The correlation function was computed using stratified sampling of the 
individual realizations of F((d) and hence of (u). In the third method, a 
higher-order equation for r was derived from Eq. (1) and solved directly in 
finite differences by two successive applications of the Fast Poisson Solver. 

The three methods were compared for accuracy and efficiency, and the third 
method was chosen as clearly superior. 

The results agree well with the global latitude dependence and local aniso- 
tropy of observed atmospheric correlation data. The value of the parameter Cq, 
which gives the best fit to the data, is close to the value expected from dynami- 
cal considerations for a baroclinic atmosphere. Differences between model and 
data are typically smaller than those between data for the Northern Hemisphere 
and those for the Southern Hemisphere. These results provide the basis for an 
optimal choice of coefficients for statistical analysis of atmospheric data. 


1 Courant Institute of Mathematical Sciences, N.Y.U., New York, New York 10012 
and Laboratory for Atmospheric Sciences, NASA/GSFC. 

2 Sigma Data Services Corporation, Greenbelt, Maryland 20771 

3 Laboratory for Atmospheric Sciences, NASA Goddard Space Flight Center, 
Greenbelt, Maryland 20771 
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A variable vertical resolution weather model with an 

EXPLICITLY RESOLVED PLANETARY BOUNDARY LAYER 
H. M, Helfand 


Sftasat A scatterometer surface wind stress observations over the world's 
oceans and seas can he very useful in the production of qlobal weather analysis, 
especially since many of the reoions that they cover are very sparse in conven- 
tional weather data. The surface wind stresses relate directly to atmospheric 
winds in the lowest few meters of the atmosphere and they are also a key to 
the determination of surface fluxes of heat and moisture into the atmosphere. 

There was a problem, however, in assimilatinq surface wind stress data 
into the standard 6LAS Weather Model. The model, with its lowest level at 
about SB mb above the ocean's surface, did not offer sufficient vertical resol- 
ution to incorporate surface data. The structure of the wind profile in the 
lowest 55 ml) of the atmosphere cannot be adequately predicted in order to 
relate the winds at the lowest level to surface wind stresses. 

However, there is enough known about the structure of the surface layer 
(which is roughly ID mb deep) to relate winds at its top to surface wind stress. 
One can appeal to the Monin-Obukhov similarity theory to establish such a rela- 
tionship. 

A new Variable Vertical Resolution version of the Fourth Order Weather 
Model has been developed so that the Monin-Obukhov formulation can be used. 

The R o-levels of this model can be arbitrarily specified, but we. are currently 
using the following sine-squared distrubution. The boundary between the 

nth and the (n+l)th laye.*s of the model is specified as 

= sin^ (nw) 

The center o'’ of the nth layer is 


0-1 lOg * Oq ) 

7 

This distribution is analytically smooth in n and it has its hiohest 
resolution at the upper and lower boundaries of the atmosphere (og*^^^ is 
quadratic in nir/lB and in (9-n)n/18 as n approaches 0 and 9 respectively) where 
such resolution is most needed. 

The packing of o-levels in the stratosphere will, in part, make up for the 
inherent lack of resolution in tne upper parts of the at.riosphere whan pressure 
of o-pressure is used as a vertical coordinate. The packing of o-levels near 
the ground should allow an explicit resolution of the PBL, for the lowest o- 
level is about 15 mb or 150 rr, above the ground, and the layer above is located 
at about 75 mb or 750 m above the ground. Even the layer above that, centered 
at about 185 mb or 1850 m above the ground, is sometimes included within the PBL. 
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The location of the lowest model layer is close enough to the ground to 
expect the Monin-Obukhov surface layer similarity theory to give a reasonable 
approximation to surface fluxes of heat, moisture and momentum (see Panofsky, 
197:1), A scheme has been devised to use this theory together with the KEYPS 
enuatlon (Panofsky, 1973) for the wind profile structure, and a similar Inter- 
polatlve formulation for potential temperature and moisture profiles In the 
unstable case and with a formula for the stable case which Interpolates between 
the Rusinqer et al . (1971) functions for a nearly neutral PBL and Clarke's 
(1970) observdtions for a moderately stable atmosphere. 

An additional mechanism for turbulence Is needed In this model to describe 
the mean fluxes of heat, moisture and momentum between those model layers that 
are contained within PRL, The model currently uses Richardson number adjustment 
scheme that Is based on the Ideas of Chang (1979) to calculate these fluxes only 
when they are appropriate. Thus the eighth or next-to-lowest layer of the model 
is considered as part of the PBL only when turbulent fluxes of heat, moisture 
and momentum enter It from below. The seventh layer Is considered as part of 
the PBL only when there are fluxes entering both the eighth layer and the 
seventh layer from below. 

The Richardson number adjustment scheme can predict turbulent fluxes 
between any two adjacent layers of the model, and so it replaces the dry con- 
vective adjustment scheme previously used. 

Currently we are developing a Mellor-Yamada (1974) type of scheme to pre- 
dict turbulent fluxes of beat, moisture and momentum between layers Instead of 
the Richardson number adjustment scheme. 


REFERENCES 

Businger, J. A., C. Wvnqaard, Y. Irumi, and E. F. Bradley, 1971: Flux pro- 

file relationships In the atmospheric surface layer. J. Atmos. Scl., 28, 
181-189. 

Chang, S. W. , 1979: An efficient parameterization of convective and nonconvec- 

tlve planetary boundary layers for use In numerical models. 

Metenro. , 18 , 1205-1215. 

Clarke, R. H., 1970: Observational studies in the atmospheric boundary layer. 

Q uart. J. Roy. Meteoro. Soc. , 96 , 91-114. 

Mellor, G. L., and T, Yamada, 1974: A hierarchy of turbulence closure models 

for planetary boundary layers. J. Atmos. Scl. , 31 , 1791-1806, 

Panofsky, H, A,, 1973: Tower Micrometeorology. In Workshop on Hicrometeoroloqy , 

n, A. Haugen, ed., Boston, Araer, Meteoro. Soc., 392 pp. 


83 


r ^SZ 25565 

A FOUR DIMENSIONAL VARIATIONAL ANALYSIS EXPERIMENT 

R, Hoffman 1 


Since we are unable tc observe the atmosphere in a detailed and accurate 
enough fashion at a single Instant of time, we must use past Information. It Is 
highly desirable to use a truly four-d1mens1ona1 assimilation method (Bengtsson, 
1975) to accomplish this reconstruction. Sasaki's variational analysis method 
(1970) for four dimensions with exact constraints. In the special case when the 
constraints are the governing equations (hereafter referred to as VAM4D), may be 
restated as a direct minimization (or nonlinear least squares) problem. Con- 
ceptually, all we need to do Is find the solution of the governing equations 
which Is closest to the observations. The works by Tadjbakhsh (1969) and by 
Lewis and Bloom (1978) are Important precursors of this study. 

As further motivation for using the governing equations as constraints, 
consider that all successful meteorological analysis procedures, subjective or 
objective, use some a priori knowledge of the atmosphere. At one extreme, 
radiosonde reports wWich violate the hydrostatic relationship, ship reports at 
continental coordinates, or reports which are grossly miscoded must be rejected. 
At a more sophisticated level, the smoothness criteria which allow the rejection 
of data causing bulls eyes In the analyzed fields Is ultimately based on a 
priori knowledge of the spatial spectra of the meteorological variables, ^he 
best analyses are subjective because a human analyst can Incorporate a vast 
amount of a priori knowledge. Therefore, we conclude that we should use all our 
knowledge of the atmosphere In analyzing data. A large part of our knowledge 
of the atmosphere Is summarized by the equations governing the atmospheric 
dynamics; a convenient form of these equations is a numerical weather prediction 
(NWP) model. It is thoroughly consistent to require the analysis to satisfy 
the best formulation we have of the governing equations; that Is, the analysis 
should be a solution of the NUP model. (All this with the following provisoes: 
the time Interval analyzed must be smaller than the time for which the model Is 
valid, and It must be understood that only scales which can be represented by 
the model are Included in the analysis.) The more complete and accurate the 
model, the better will be the analysis; therefore, VAM4D should. In iqy opinion, 
eventually prove to be the most satisfactory objective analysis method. 

To test VAM4D, we have performed an Idealized continuous assimilation exper- 
iment using a very simple model which has been described In detail by Hoffman 
(1981). The small alterations to the imposed heating function discussed by 
Hoffman and Kalnay (1981) are also Included here. The experiment is an identi- 
cal twin experiment; nature and the dynamical model are the same quasi geostrophic 
model. The observing system contaminates nature with errors which are Indepen- 
dent for each spectral component and which have white noise temporal statistics. 
The standard deviation of the error added to each spectral component Is 7.5% 
of the climate variance of that variable. Stream-function Is not observed, only 
temperature Is observed. The assimilation period Is 4 days; observations are 
available every 12 h, and the entire experiment is 128 days long (I.e., there 
are 32 cases). 


^ National Research Council /Resident Research Associate 
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The VAM40 finds the model solution which qives the best fit to the obser ed 
data at t s 0 h, -12 h, -24 h«... -84 h and to the analysis at -96 h which is 
the result of the previous assimilation. (To avoid startup problems, the first 
case uses a perfect analysis at -96 h.) The results of the experiment are 
shown in Fiq. 1, which is the average over the 32 cases of the analysis error, 
as a function of time. We have normalized the analysis error by the observa- 
tional error which would have been present if stream-function had been observed 
(Fiq. la). The error is D, defined by 

* (6D^ + 2"^ (4(0^0^)^ a 2*^(4^)^ 


where e, g and ^ are the 500 mb nondimensional potential temperature, static 
stability and stream-function, the overbar is a horizontal average, 6X is the 
error in x and a * The analysis error is smaller than 25% of the observa- 
tional error; the squared error is reduced by ~95%. Note that the analysis 
error is smallest in the center of the assimilation period close to the center 
of mass of the data. Also, the analysis error is smaller at t = -96 h than at 
t = 0 h reflecting the information content of the t = -96 h analysis. Fig. lb 
is identical to Fiq. la except that the term involving ^ has been neglected 
(i.e., a = 0), so that D is now the error projected on the observable manifold. 
Apparently, the stream-function field is reconstructed as well as is the temper- 
ature field. 

In summary, we have shown for a perfect model and a simple (but nonlinear) 
system that the four-dimensional variational analysis method using the governing 
equations as exact constraints works successfully. First, the method is stable 
in an assimilation cycle. Second, it reconstructs the unobservable variables; 
in this case, no velocity data was observed. Third, the analysis errors are much 
smaller than the observing system errors. 
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Fig. 1. The analysis error as a function of time. The quantity plotted, 

R, is the ratio of the analysis error Danal observational 

error Dahs for (a) temperature and streamfunction (a = 1, Dobs ? 
556x10"®; and for (b) temperature alone, (a » 0. OqJjj » 261x10"°). 
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SASS WIND AMBIGUITY REMOVAL BY DIRECT MINIMIZATION^ 


Of) 


Ross N. Hoffman 2 


A qridded surface wind analysis Is obtained by minimizing an objective 
function, the magnitude of which measures the error made by the gridded analysis 
In fitting the SASS wind data, the conventional surface wind observations, and 
the forecast surface wind field. The ambiguity of the SASS winds Is then 
removed by choosing the alias closest to the analyzed wind. Because minimizing 
the objective function Is a problem of nonlinear least squares, the minimizing 
gridded analysis Is not unique, and a good first quess Is necessary to assure 
convergence to a reasonable solution. A good first guess may be generated by 
performing the analysis In stages. 

Illustrative results are shown for a limited region In the North Atlantic 
containing the QE II storm and for a limited amount ( ~ 12 minutes) of data 
observed near the synoptic time 1200 GMT 10 September 1978 (Fig. 1). Within 
the SASS data swath the resulting gridded analysis (Fig. 2a) Is a reasonable 
representation of the surface wind and Is not sensitive to the forecast surface 
wind field. The analysis Is sensitive to wind directions reported by ships 
near the center of the storm. The wind circulation center present In the fore- 
cast Is moved by the analysis. The resulting deallased SASS wind directions 
are noisy (Fig. 2b). 


' Presented at the A6U 1981 Spring Meeting. Baltimore. Maryland. 25-29 May 1981. 
This research will be reported In full In Mon. Wea. Rev. 
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Fig. 1. The SASS data used in the analysis. All observations are near 1200 
GMT 10 September 1978. Arrows originating at the data locations 
are proportional to wind magnitude and plotted away from the wind 
direction. Octagons are plotted at SASS nadir locations. (Detail shown). 
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Fig. 2. Minimizing analysis (a) and corresponding dealiased SASS winds (b) 

obtained by the three stage analysis procedure using the NMC analysis 
as the forecast wind field. The NMC conventional data base was used. 
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A SIMPLE ATMOSPHERIC MOOEL ON THE SPHERE WITH 100% PARALLELISMI 

E. Kalnay and L« L. Takacs 


1. INTROniJCTION 

One of the nost important applications of the forthcoming massively para- 
llel computers is numerical weather prediction and clir^ite simulations. The 
need for hioher spatial resolution in models for weather prediction, and of 
very long integrations for climate studies, requires computing capabilities 
several orders of magnitude larger than those presently available. Only mass- 
ively parallel computers will he able to fill this need. However, it will be 
also required to design models with a structure that can take advantage of the 
parallel configuration. 

We present here a very simple shallow water equations model, in which 
spherical geometry is transformed into a formally doubly periodic configuration, 
which can he easily implemented on a parallel computer at 100% parallelism. 
Shallow water equation models have been traditionally used to test numerical 
weather prediction schemes because their structure is analogous to 3-dimensional 
models of the atmosphere. 

This model has been tested on the serial Amdahl computer at the NASA/GSFC 
Goddard Laboratory for Atmospheric Sciences (GLAS). Its implementation on a 
simulator of the Massively Parallel Processor (MPP) being developed at NASA/GSFC 
bv nr. navid H. Shaeffer is part of a joint project with Drs. David Randall and 
Milton Halem. 


?.. OESCRIPTION OF THE SCHEME 

The shallow water equations on spherical coordinates are (Figure 1) 

ah = - 1 T ~a(hu) + a(hv cos4> )~ T 

at a cos<j> J__ax a<^ j[ 


a(hu) 

at 


at 


1 T~ a(hu)u + a(hv cos» }u ~T 

a coS(j. _L_ax a4» J[ 

qh ahT + ( f + u tanj> j hv 

a coS(|» ax ^ a ^ 


a cos<^ 


1 T” a(hu)v + a(hv cos<i> Ivn * 

os(^ 2_ax a(|) I 


gh 


ahT - ( f + u tan^ i hu 
a a<> ^ a ^ 


where f * 2sin^. 


^ Published in Advances in Computer Methods for Partial Equations - IV. 
Published by IMACS, 1981. 


( 1 ) 

( 2 ) 

( 3 ) 
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Lftts consider a spherical arid confiquration, with no qrid point at the pole 
(Fiq. ?a). We take "complete" meridians as great circles formed by a meridian 
on the "liqht side" of the sphere with a longitude \ and its continuation on 
"dark side", with a longitude X+ir. Me let "latitude" ^ be periodic along a com- 
olete meridian: - v/?. < ♦ < n/2, cos ^ > 0 on the "light side" and n/2 < 4 * ^ 3w/2, 

C 0 S 4 . < 0 on the "dard side." Similarly, sin^ > 0 in the Northern Hemisphere and 
< n in the Southern Hemisphere. The geometry becomes thus doubly periodic-like, 
with a configuration that is ideally suited for a massively parallel processor 
(Figure 2h). »^e assume that each processor is connected to its four closest 

neighbors, with connections "over the poles" as indicated in Fig. 2b. 

We use a nonstaggered grid, with all variables defined at the grid points 
(circles). Let 


«x^i.i 


^i+l/2.i-^i-l/2 

AX 




ij 


^i+l/2+f'i-l/2.i; so, 



eg, &x 




^i+ij“^i-i.i 




and similarly for *. Note that ,i=l,..., JMAX/2 on the liqht side, j=JMAX/2+l, 
..., JMAX on the dark side, AX = 2it/IMAX; A<t> = 2ir/JMAX, and (j»^=(.i-l/2)A4»-ir/2. 

2 

Then a very simple parallel scheme is the following: 

STEP 1 : Longitudinal cycle, performed for j=l,..., JMAX/2 with i cyclic 

from 1 to IMAX: 


^ 

3t ^ ^°S^ 

3hii = - ^ 

3t (x) a cos<|) 



3hv 

^ M 


+ r f + u tan<> ^ hv 
a 


1 


a COS 41 




(4a) 


(5a) 


( 6 a) 


Note that this step, as all the following ones, is performed simultaneously for 
all processors. 
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STEP 2: Latitudinal cyd •. performed for 1=1,..*, IMAX/2 with j cyclic from 1 

to jHAx. 

Chanqe siqn of wind components on “dark side" (l.e., multiply u, v, 
hu, hv by d, where 


( 1 on the light side 

d= ) 

/ -1 on the dark side 
2b: 


3h = “ 1 

66 hv COSA i 

(4b) 

3t (♦) e coS(t 

L i 


3hu 

♦ ♦ 


3t (®) = - 1 

6(|> (hv coS(|) u ) 


a cos^ 

L J 

(bb) 

3hv 

♦ ♦ 

r J 

3t (i) = - 1 

6(|> (hv coS(^ V ) 

- qh 6^ hr 

a coS(^ 

L J 

a 

■ ( f iL 

tan^ ] hu 
a 

(6b) 


Notes I 

.♦ 

i) hv coS(^ at the poles is not computed but defined as zero. 

11) 11 , V, cos ^ and 3/3^ chanqe signs on the dark side. As a result, the 
equations remain unmodified. 

2c : Same as 2a 

STEP T : Update the fields h, hu, hv using a suitable time step. For example, 

if leap frog is used. 
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3. IMPLEMENTATION OF THE MODEL 

Since no MPP or MPP simulator is yet available, we have implemented the 
proposed model on the Amdahl serial machine. The parallel structure was simu- 
lated by performinq STEP 1 in a FORTRAN array as io Fiqure 3b, chanqing the 
siqns of the winds rn the dark side. 

The model worked very well, qivinq results similar to those obtained with 
a standard shallow water equation model with a grid point at the poles. 


4. DISCUSSION 

The simple scheme proposed here as a prototype of a 3-dimensional model 
has the following characteristics: 

a) It is 100% parallel, so that on the MPP it should be IMAX*JMAX faster 
than on a serial machine of the same speed. In other words, with the exception 
of (hvcos^)^ at the poles, all operations are identical for all grid points and 
can he performed simultaneously. With the current values IMAX=72, JMAX»46 used 
in the GLAS models, this represents a factor of 3000 in speed. 

h) It conserves mass and energy. 

c) It is almost identical to a second order version of the Fourth-Order 
GLAS global atmospheric model (Kalnay-Rivas, et al . , 1977), the only difference 
being the lack of a grid point at the poles. 

d) It can he immediately generalized to Fourth-Order if four more connec- 
tors are included (Fiqure 4). 

e) The ISth order Shapiro filter currently used in the Fourth-Order GLAS 
Model can be trivially implemented at 100% parallelism. The Fourier filtering 
at high latitudes that eliminates short unstable waves near the poles requires 
the development of an appropriate FFT algorithm. 

f) If each processor has IK memory, as the MPP will, the implementation 

of 3-dimensionality can also be done at nearly 100% parallelism. This is because 
all the variables in one atmospheric vertical column (a few hundred with reason- 
able vertical resolutions) can be contained in one processor. Since the compu- 
tations within each vertical column are independent of other columns but similar; 
they can also be performed simultaneously. 

The impact that the utilization of massively parallel computers will have 
on weather and climate research may be as large as that produced when the 
development of computers during the late 1940's and SO's made numerical weather 
prediction feasible for the first time (Charney et al. , 1950). For example, 
current GLAS models require about one hour of computer time to simulate one day 
of atmospheric flow. At 100% parallelism, the same models would simulate eight 
years in one hour. 

It is very remarkable that in his monumental work, Richardson (1922) 
developed not only many of the basic physical and numerical ideas presently 
used in numerical weather prediction, but also the concept of a massively 
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parallel processor (Figure 5). Note that in his "fantasy", there are 2000 para- 
llel processors (columns), and 32 (human) computers for each processor. Note 
that he also suggested the use of simpler models for the development of improved 
procedures* In the preface of his book Richardson wrote: "The scheme [for 

weather prediction] is complicated because the atmosphere is complicated . . . 
Perhaps some day in the dim future it will be possible to advance the computa- 
tions faster than the weather advances and at a cost less than the saving to 
mankind due to the Information gained. But that is a dream." 

The application of massively parallel processor tu atmospheric modeling 
will truly fulfill Richardson's dream. 
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Figur e 4: Owinectlivis for a fourth orilor «ichifae. 


A ronrcAKT rACToRY 

C» ll/« THE iPEtn AND ORGANIZATION OF CXJUFLTINO 

It ttiok «• th« bwt part of «i wmIu to drmv up tko computing formt and to work 
out tko nov diotributioa in two oortieoJ oolurano for th» 6rit time. My ofioi woi n 
boop of hoy in n eoid runt billot With pmetioo tbc work of on ovora^ eomputtr 
might go porbopo too timoo laotor. If tbo time-otop wtt« 1 boort, tbon 37 individnolo 
oould juot eomputo two poioto to ao to kocp pace with tbo vaatbei, if wc stiow ootbii^ 
for tbo rtry graat gain in tpaod wbieb io invariably notioad whan a oomplioatad 
oporation ia dividod up into timplor parta, upon wbieb individua' apoeialiaa. If tbo 
oo-ordioata cboquor wart 200 km aquarc in plan, tbaro would be 3200 oolumna oo 
tba eomplata map of tbo globo. In tha tropiea tbo waatbar ia oftao foreknown, ao 
that wa may toy 2000 aetiva oolumna. So tbat 32 k 2000 > 64,000 eomputara would 
bo needad to raeo tba weatber for tba wbola globe. Tbat ia a ataggering Sgum. Par- 
bapa ia aoma yean' time it may ba poasibla to report a aimplifieatioo of the p r o c ma 
But in any eaaa, tba otganiaation iudioatad ia a ecntral foraeaat- factory for tba wbola 
globe, or for portiona aitending to bouodariaa wbaro tbc weatber ia ateady, with indi- 
vidu'J eomputara apoeialiting on tba aeparate aquationa Lat ui bopa for their ankaa 
tbat they are moved on from tima to time to new oparatiooa 

After to mock bard reaming, may ona play with a faolaay t ImagiiM a large ball 
like a tbaatre. escapt tbat tba eirclea and gnlJariat go right round through tha apaea 
uaually ocenpiod by tba atage. The walla of thia chamber arc painted to form a map 
of the globa Tba ceiling rapraaanta the north polar rogiona, England ia in tba gallery, 
the tropica in the upper eirele. Australia on the dreaa circle and the antarctic in the 
pit A myriad eomputara are at work upon tha weather of the part of tbc map where 
aacb aita, but aacb oomputar attaoda only to ooa aquation or part of an equation. The 
work of aacb region ia eoardinatad by an ofioal of higher rwak. N umermia little night 
signs * display the instantaneous values ao tbat neighbouring eomputeri can read them. 

number is thus displayed ia throe adjaeeot aooas so as to maintain eommunioa- 
tioo to tha North and South oo tha map. From the loor of the pit a tail pillar riaas 
to half tha haight of tbo halL It earriaa a large pulpit oo its top. In this aita tbc 
man in charge of tha whole thaatio ; ba is surrounded by aevatml sseistanta and ssas- 
aaegerm. Ona of bis dutiaa ia to maintain a uniform speed of prograas in all porta of 
tha globo In this raapoet ba is lika tba conductor of ao orchestra io which tba inatru- 
manta are alida-rulas and edculating maebinao But instead of waving a baton ha turns 
a beam of may light upoo any region that is running ahead of tha rest, aod a beam of 
blue light upon tboas wbo are behindhand. 

Four senior clerks in tha osntrul pulpit are eollacting the hitura waatbar as East as 
it ia being computed, and daspatchiog it by pneumatic carrier to a quiet room There 
it will be coded and telaphonad to the mdio transmitting station. 

Meaamgars carry pilas of used computing forms down to a sWebousa in tbe cellar. 

In a neighbouring building tbers ia a resaorch departmant. wbarc tbay invaiit 
improvamanto But there ia much etpanmenting oo a small acala bafora any changv 
ia made ia tha eoroplas rootina of tha computing thaotre In a basement an enthusiast 

is observing eddies in tha liquid lining of a huge spinning bowl, but ao far tba antb- 
matie proves tha better way. In another building arc all the usual financial, 
correspondeoea and adminiatrativa oAera Outside arc playing fields, iiouaes. mountains 
and lakes, for K was thought that those who compute the weather should braatha of 
it freely. 


F igure 5: 


"A Banscast Fbetory", fron Richariaon (1022). 
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ON THE EFFECT OF HIGH LATITUDE FILTERING 
IN GLOBAL GRID POINT MODELS^ 


L. L. Takacs^, R, C. Balgovind^, M. Iredell ^ and E. Kalnay 


1. INTRODUCTION 

The process of filtering at high latitudes In global grid point models Is 
fairly conmon In meteorological studies. This Is done in order to avoid the 
use of a prohibitively short time step required by fast moving waves near the 
poles. In the Arakawa-MIntz general circulation model (Gates et al., 1971), 
smoothing was performed on the zonal mass flux in the divergence terms and 
on the zonal pressure gradient in the prediction equations, thus eliminating 
the fast moving gravity waves In the zonal direction. Vanderman (1972) applied 
a weighted average in the east>west direction to the tendencies of the dependent 
variables near the poles to eliminate short wave disturbances. Umscheld and 
Sankar-Rao (1971), Arakawa (1972), and Kalnay>R1vas et al. (1977) used Fourier 
filtering to eliminate the unwanted waves. We have found that when these 
filtering techniques are applied In a coarse resolution grid point model, the 
non-linear effects responsible for change In the mean zonal flow become 
poorly simulated, resulting In spurious energy and momentum transfers. 

In this paper we will examine the effects of this phenomenon and its 
relationship to various filtering techniques. 


2. MODEL DESCRIPTION 

The model used to examine the effect of high latitude filtering is based 
on the shallow water equations on a sphere given by 
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1 Published In Preprint Volume: Fifth Conference on Numerical Weather Predic- 
tion, Nov. 2-6, 1981; Monterey, California. Published by the American Meteor 
ological Society, Boston, Massachusetts. 

2 Sigma Data Services Corporation 
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( 3 ) 


3(h) + 1__ 3(hu) + 3(hvcosfl) 

at acos6 ax ae 


0 


where h Is the height of the atmosphere, u and v are the zonal and meridional 
wind components, and f Is the Coriolis parameter. The latitudinal and longitu- 
dinal directions arc given by e and x respectively, while a Is the radius 
of the earth. For simplicity, bottom topography has been set to zero for 
these experiments. The equations at the poles follow those of the GLAS 4th 
Order General Circulation Model (Kalnay-Rivas et al. , 1977). 

All variables are defined on a non-staggered 'A' grid. Spatial derivatives 
are approximated by 4th order differences, while the time derivatives are 
approximated by the leap-frog scheme. A simple forward step Is applied with 
At equal to 1/4 of Its original value every 20 time steps to avoid the separa- 
tion of the solution produced by the leap-frog scheme. 

In these experiments the model was run under two general conditions, the 
first using a very small time step allowing all resolvable waves to remain com- 
putationally stable, and the second using a time step defined by the most un- 
stable wave located at 60 degrees latitude. Poleward of 60 degrees various 
types of filters and smoothing operators were used to control the amplification 
of the shorter waves. 


As the Initial condition for these experiments, we followed Phillips (1959) 
by taking an Initially non-dl vergent velocity field described by the stream 
function 



-Wsine + Kcos^^e sine cos(Rx) 


( 4 ) 


together with a height field derived by solving the Inverse balance equation. 

The constants K and VI In the stream function equation were both equal to 7.848xl0‘° 
sec" , while the wavenumber In the zonal direction, R, was equal to 4. The mean 
mean height of the atmosphere was equal to 8000m. 

As discussed by Hoskins (1973) and Ooron et al. (1974), this Initial con- 
dition produces a super-rotation of the atmosphere In addition to a Rossby- 
Haurwltz wave which moves westward relative to the fluid. For a non-divergent 
barotropic atmosphere, the total solution will move from west to east without 
change of shape. Since we are using a free upper surface, however, some changes 
In the solution can be expected due to the divergent nature of our model. In 
these experiments we are Interested In how closely the time evolution of the 
filtered run compares with the non-flltered run. 

3. RESULTS 

The following results were obtained by Integrating the above Initial con- 
ditions on a coarse grid having a resolution of 10 degrees In longitude and 
8.18 degrees In latitude. Figure (1) shows, using a rectangular latitude-longi- 
tude projection, the Initial height field and subsequent heights at 100, 200, 
and 300 hours of Integration for the nonflltered case, while Figure (2) shows 
the same for the filtered case. In this run, a Fourier filter with a response 
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similar to that used in the GLAS 4th Order Model was applied to the fields u, 

V and h at every time step. This filter will be discussed In detail in 
section (4). 

From Figure (1) we see that the initial disturbance moves eastward at a 
rate of about 90 degrees per 200 hrs. with relatively minor change. The u 
field (not shown) shows some slight tilting of the trough-ridge lines which 
results in u-momentum being transported meridionally. However, the overall 
solution appears close to that given by the non-divergent equations. 

Figure (2), however, shows drastic changes occuring between 100 and 200 
hours. The tilting of the trough-ridge lines becomes very pronounced in the 
northwest-southeast direction, and can be readily identified at 200 and 300 
hours. The appearance of a closed contour low is evident at 200 hours, 
while the amplitude of the wave steadily decreases from 200 to 300 hours. The 
zonally averaged u field, which was initially positive at all latitudes, 
became negative at mid-latitudes in conjunction with the cut-off low. 

Figures (3) and (4) show the mass weighted u-momentum transfer due to 
the zonally averaged eddy momentum flux at 0, 100, 200 and 300 hours. We 
can see from figure (3) that there is a transfer of u-momentum equatorward in 
the sub-tropics and a vacillating transport at mid-to high latitudes. In the 
filtered run, figure (4), there is no poleward transport at all but a very 
strong equatorward transport at mid-latitudes. As a result, the energetics 
of the two runs became vastly different since more energy was being fed into 
the zonal flow for the filtered case. Figures (5) and (6) show the zonal and 
eddy kinetic energy. From figure (5) we see that there is very little change 
in the energetics of the non-filtered solution once the initial state becomes 
balanced. In this run the zonal kinetic energy increased by 13% while the 
eddy kinetic energy decreased by 9% after 300 hours. The filtered case, however, 
shows significant changes at every stage. After 300 hours the kinetic energy 
due to the zonal flow has increased by 75% while that due to the eddy flow 
has decreased by 75%. 

These figures show the effect of high latitude filtering in our model. 

Since filtering can be thought of as a latitudinal dependent dissipation 
function, two further experiments were performed in which a non-fiUered 
model was run with the addition of linear damping in the u-and v-momentum 
equations. The time constant t associated with the damping was equal to 7 
days. In the first run j was constant in latitude while in the second run 
T was multiplied by sin^ (e), thereby simulating the effect of the 
filter. 

The results from the first run showed the Rossby-Haurwitz wave slowly 
moving to the east with steadily decreasing amplitude at all latitudes. No 
tilting of the trough-ridge lines was evident nor was there any formation of 
the cutoff low. The second run, however, showed results identical to those 
of the filtered runs. Strong tilting of the trough-ridge lines was present 
as well as the cut-off low in mid-to high latitudes. The results showed 
that the deformation of the Rossby-Haurwitz wave was a realistic response 
of the model to a latitudinal dependent dissipation function introduced by 
the filter. 
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These results can be understood in terms of Eliassen and Palm (1960) 
theory. The high latitude filtering Introduces a sink of wave energy, which 
In turn Induces a poleward flux of eddy energy. In accordance to the Eliassen* 

Palm theorem, this must be accompanied by an equatorward flux of westerly 
momentum, resulting In the strong subtropical Jets observed In the filtered 
numerical results. 

4. FILTERING TECHNIQUES 

In order to try to eliminate the errors produced by the filter, variations 
in the filtering technique and In the filter's frequency response were Investigated. 
As pointed out earlier, some atmospheric modellers filter other quantities such 
as the tendencies or the pressure gradient terms rather than the field variables 
themselves to eliminate the fast waves near the poles. We compared the results 
obtained from using these filtering procedures and found that all results showed 
the same behavior. It appeared that the deformation caused by the filter was 
independent of the filter's source. 

As noted earlier, the filter response used in this experiment was similar 
to that used by the GLAS 4th Order Band Model (hereafter referred to as the 
operational filter). In the filtering technique, the field variables u, v and 
h were Fourier decomposed into their wave constituents which were then altered 
by a wave dependent damping function. The frequency response of the damping 
function was empirically selected after performing experiments with the GLAS 
4th Order Model. The latitudes used by the filter were twice the distance 
from the poles than in the GLAS 4th Order Model due to the coarse resolution 
of the shallow water model. 


Figure (7a) shows the frequency response of the damping function used in 
the operational filter as a function of wavenumber. Primarily, this response 
uses a linear damping curve over six wavelengths until complete truncation is 
achieved. In addition to this response, three other damping functions were 
tested. As shown by Arakawa and Lamb (1977), a damping function based on the 
linearized gravity wave equations on a sphere was derived using a staggered 
grid and centered 2nd-order spatial differences. This function is given by 

F(k)=^ . cos e (5) 

A0 sin( k AX ) 

where k is wavenumber. The actual value used as the damping coefficient in 
the filtering procedure would be the minimum of 1 and equation (5). The 
response of this filter is shown in figure (7b). 

Using a similar type of analysis as Arakawa and Lamb, a damping function 
using a non-staggered grid and 4th-order spatial differences was also derived 
and tested. The response of this function is given in figure (7c), while 
Its functional form is given by 


F(k) = ^ 

A6 


cosfl 

4 sin(kAX)-l sin(ZkAX) 

■3 "5 


( 6 ) 


Note that In this type of filter the shortest waves remain unchanged. 


Finally, as an approximation to the Implicit bi-harmonic diffusion used 
in the ECMWF model, a damping function derived from the longitudinal component 
of the bi-harmonic diffusion was also tested, given by 

-1 

V (2cos(2kAX)-8cos(kAX) ♦ 6) + (7) 

_ (acos 9 AX)^ _]_ 

where the constant v was taken as 4.0x10^® mV sec. The frequency response 
of this function is shown in figure (7d). 

By using the convolution property of the Fourier transform, the damping 
coefficients in wavenumber space may also be viewed as the weighting co- 
efficients in grid space. These are shown in figure (8a-d) for the grid point 
nearest to the pole. 

All of the above filter types have been tested using the initial conditions 
described earlier, and all of these experiments produced results similar to 
those of the operational filter. The results obtained using the damping func- 
tion associated with figures (7c) and (8c) were slightly worse due to the in- 
clusion of the poorly simulated short waves. It again appears, however, that 
the filtering procedure itself causes the phenomena to exist and has little 
dependency on minor variations in the frequency response of the filter or on 
the method of application. 



5. CONCLUDING REMARKS 

We have seen from the preceding sections that high latitude filtering in 
a coarse resolution grid point model produces spurious energy and momentum 
transfers which affect the solution at all latitudes. This phenomenon seems 
to be independent of the details of the filtering procedure and occurs even 
if only the time derivatives or the divergence terms are filtered. In coarse 
resolution global models maximum integrations seem to be limited to one week 
or so unless sufficiently small time steps are taken to eliminate the need 
for filtering. 

It should be pointed out that the same experiment was also carried out 
with double resolution (i.e. 5 degrees in longitude X 4 degrees in latitude). 
Using a time step of 10 minutes and a filter identical to that of the GLAS 4th 
Order Model, the results were very similar to those of the nonfiltered coarse 
grid model and showed no evidence of spurious transfers of energy and momentum. 
One possible reason for this may be the fact that in the Rossby-Haurwitz wave 
of zonal wavenumber 4, there are strong nonlinear interactions of wavenumber 
8. In the coarse resolution model, this is close to the 4 AX wave (zonal wave- 
number 9) which is strongly filtered. In the high resolution model, the 4 
AX wave has zonal wavenumber 18. Thus the coarse resolution model may have a 
poorer simulation of the nonlinear effects when compared to the fine resolution 
model. 

These results point out a serious deficiency of coarse resolution grid 
point models used for climate simulations. We are presently stu((ying the 
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effect of high latitude filtering in higher resolution models used for weather 
prediction as well as in models using implicit time steps. 
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Figure 1. Initial height field (top) end eub- Figure 2. As In Figure 1 except for the filtered 

sequent height fields at 100, 200 run. 

and )00 hours for the non-flltered 
run. 
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Figure 3. Hess weighted u-aoaencua transfer Figure 4. As In Figure 3 except for the filtered 
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Figure 5. Hass weighted kinetic energy 

(X 10^ B^/sec^) for the non-flltered 
run. 
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run* 


103 
















Figure 


ORIGINAL PASS 13 
OF POOR QUALITY 


Di^W CWIkiotn n Mfcne Span 


CMileaaM m Ctcrf Sean 




7. Deaplng functione uaed In Fourier fil- 
tering for the 3 latitudes nearest to 
the poles for a) operational filter, 
b) 2nd order filter on a staggered 
grid, c) 4th order filter on a non- 
staggered grid, d) bl-haraonlc dif- 
fusion filter. 


Figure 8. Weighting coefficients In real space 
derived froa the daaplng function 
used at the latitude nearest to the 
poles for each filter described In 
Figure 7. 
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SHALLOW WATER EXPERIMENTS USING AN IMPROVED 
KREISS-OLIGER 2-4 SCHEME 

L. L. Takacs and R. C. Balgovind 


1. INTROnuCTION 

Kreiss and Ollqer (1973) Indicated the advantages of the 2-4 scheme for 
meteorological purposes. This scheme consists of the second-order leapfrog 
scheme for time derivatives, and fourth-order space differencing for spatial 
derivatives. Abarhanel and Gottlieb (1980) have shown that by modifying the 
Kreiss-Oliger scheme in a manner similar to that done by Wilson (1978) to the 
multi -dimensional leapfrog scheme, a substantial improvement of the maximum 
time-step could be obtained. As an experiment we applied the Aba rbanel -Gottlieb 
scheme to a global 1 -level shallow water model and compared the results to 
those using the standard Kreisf-Oliger sche»ne. 


2. REVIEW OF SCHEMES 


For a system of equations defined by 

la = iL + (U 

3t ax ly 

where g, F, and G are m-componet vectors, the Kreiss-Oliger (hereafter referred 
to as KO) scheme is given by 


g 


n+1 


i.i 
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^ ^i,j+2 ■ ®i,j-2 


( 2 ) 


4ax 


4ay 


The Abarhanel -Gottlieb (hereafter referred to as AG) scheme expands the 
domain of influence of the above scheme by including off-axis grid points. 
This scheme is given by 


q 



n-l 

Qii + 2At { 


^ “ ^i+2,.H2 ^ ^i+l.j+2 • ^1-l.j+2 ^1-2,j*2 
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^ ^ '"?+!. .i+1 “ ^ ’ ^?-2..i+i 

• Pf.2J - 2 Pi.lJ - 2 FU.S * F?.2J 

^ ^i+2J-1 * 2 Fi+i j.i - 2 - I'i-aj-l 

n n n n 

' ^i+2,.i-2 + Pi+1,.1-2 - fl-l,.i-2 + ^1-2..i-2 ) / 12 ax 

^ (3) 

< ‘ ^‘i+2,.i+2 * ®i+2J+l - ®i+2,.i-l ®1+2,j-2 

^ ^'1+1 ,.1+2 ^ 2 G!}+^ j+i - 2 Gi+ij-i - <'*1+1 J-2 

• «1,.i+2 * 2 . 2 + G?^^.2 

''’1-1 ,.1+2 + 2 - 2 g5_i j.i - g![.ij.2 

- <*1-2J+2 * <*1-2J+1 - <'1-2,j-l ■*■ <*1-2,j-2 )/^2 ay } 

Since the fluxes F and G are already calculated at all qrld points for the 
KO schene, the extra terms In the AG scheme do not constitute any additional 
flux calculations but only additional flux sun«!)at1ons. 


3. TIME STEP ANALYSIS 

Since our main Interest here Is to examine the effectiveness of the AG 
scheme when applied to the shallow water equations on a sphere. It Is useful to 
perform a time step analysis based on the linearized gravity wave equations. 

For clarity, the analysis will be done In both cartesian and spherical coordi- 
nates. 


The gravity wave equations In cartesian coordinates are given by 


du B 
at 

. g ah 
ax 

(M 

av B 
al 

• q ah 
ay 

(5) 

ah B 
aT 

- H c iu + av ] 
ax ay 

(6) 


where h Is the height of the atmosphere and H Is the corresponding scale 
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helqht, and u and v are the two-dlraenslonal wind components* 
Assufflinq wave solutions of the form 


1(e xi+eyi) 

. ’"..I . * 1 


( 7 ) 


where ex “ k^x » 2 hAx/Lx 

ey = lAy » ZnAy/Ly 


and using the appropriate difference schemes* the maximum allowable time-step 
qoverninq stability is qiven by 


At, 


max 


^ 1//^ T“ /Wxv^ ^ -r 

Ax' 'Ay J, 


- 1/2 


max 


where for the KO scheme we have: 


Wx = 4 sin ex" 1 sin2ex 


Wy * 4 sin ey" 1 sin2ey 

1 ? 


while for the AG scheme we have: 


Wx * 1 C 2( sinex“Sin2ex)cos2ey 


+ 2{3sinex+s1n20x) cosOy 
+ (3sinev-sin20v) ] 


Wy = 1 C 2{ siney-sin20y)cos20x 
1? 


+ 2(2sin0y+sin2ey)cos0x 
+ i2siney-sin26y) 


For slmplir.ity* let dx“dy=d. Then equation (8) becomes 


At, 


max 


d//^ 


j^(w| ♦ My/n^x ^ 


- 1/2 


( 8 ) 


( 9 ) 

( 10 ) 


( 11 ) 


( 12 ) 


(13) 


f 
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For thA KO schonp, and Wy are mutually Independent and thus each can be 
■^ximized separately. This is not true, however, for the AG scheme and therefore 
the maximization of the sum is required. Ooinq so we find that 


KO: 

*'*^x ■ ‘^ymax ~ *^x * 
when 6x = 9y = 1.804 

(14) 

while 



AG: 

(■< " ‘^y)max = ‘4x = 1*882 



when Qj( * 1.804 ; 6y = 0 

(15) 


OR Ox = 0 ; Oy = 1.804 



Note t'<at for the AG scheme, the non-dimensional wavenumber which qets zeroed 
out will he the one associated with the larqest grid length. Here, since 
dx=dv, the preference is unbiased. Hsinq (14) and (15) in (13), we obtain 

AtAf, = 4t|jo (16) 


or, an Increase in time step of /7 can he expected when using a uniform rectan- 
gular grid with equal spacing. 

We can now repeat the analysis by writing the governing equations in 
spherical coordinates. 


au 

= . JL_ 

ah 


on 

at 

acnS(^ 

ax 



av 

= . g 

1!L 


(18) 

at 

’ a 

a<^ 



ah 

= - H 

. T 

‘ au + a(vcos^) ~r 

(19) 

at 

acos^ [_ 

_ ax a<ji J_ 



where a denotes the radius of the earth. Using the appropriate difference 
schemes we find 

Atp,.„ = acos<^//W 1 + cos^(ji W»W<t l T 

AX WT'maxJ. (20) 


where Wj^ and Wx are defined analogous to equations (9)-(12). The term W, 
replaces with 9^+A^. Note again that for the KO scheme, and an 
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mutually Independent while for the AG scheme they are not. Also, the maximum 
time step has now become latitude dependent, reflecting the changing grid length 
In the zonal direction as one approaches the pole. 

Tables (1) and (2) show the maximum allowable time steps in seconds for the 
two schemes uslno a course grid (10** longitude x 8.18** latitude) as a function 
of latitude and scale height. Tables (3) and (4) show the same using a fine 
grid (8** longitude X 4** latitude). For the KO scheme a uniformly decreasing 
time step Is observed from the equator to the pole as expected. The AG scheme, 
however, has a constant time step up to a critical latitude, and then begins to 
decrease as one continues toward the pole. This critical latitude Is the latl* 
tilde where acos(a)AX«aA^ (~35®). Below this latitude the grid length In the 
zonal direction Is larger than the grid length In the meridional direction, 
while above the critical latitude the reverse Is true. This causes a latitudinal 
bias as to which non-dimensional wavenumber, 6x or 6^, gets zeroed out when the 
time step is maximized. The result, therefore. Is a constant time step associ- 
ated with Aa below the critical latitude, and a uniformly decreasing time step 
associated with AX above the critical latitude. As can be seen, the maximum 
gain the AG scheme can have over the KO scheme Is right at the critical latitude, 
given by a factor of ✓?. In practice, however, a much siraller time step than 
that associated with the critical latitude Is needed In order to Inhibit the 
growth of short waves near the poles. Thus, the resulting gain of the AG scheme 
over the KO scheme Is reduced. 
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TABLE 2 
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4. THF MOREL 

The moflel usefi tr> experimentally test the effectiveness of the AO scheme 
is based on the cotmilete shallow water equations on a sphere qiven by 


3(hu) 

at 


1 


acosi^ 


I 


a(hu) 

ax 


a(hu)vcos» 
a<) __ 


(21) 


_ ( f + utan» ) |,v t 
d 




a cos<^ 


ah 

ax 


a(hv) + 

at 


L. r 

acos<i> I 


a(hv) ♦ a(hv)vcosi) ”f 

ax ^ 


+ { f + utan» ) hv + qh 3h = 0 

d d 3^ 


( 22 ) 
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3h 4. 1 I a(hu) + a (hvcos») 

at acos<|) I ax a<|) 


0 


(23) 


All variables are HefIneH on a non-stan'iera*! grid with a 36x23 grid point 
domain (10® longitude X 8,10® latitude). 

The model uses Fourier filtering of u,v, and h at every time step for 
latitudes |♦|>6^ deg. to Inhibit the growth of the shorter waves near the 
poles. In addition, a 16th order Shapiro filter Is applied to the field vari- 
ables In both the zonal and meridional directions once every hour. With these 
filtering procedures, the total potential enstrophy remains quasi -constant for 
all stable runs In agreement with the results obtained by Kalnay and Hoitsma 
(ld7g). The model equations at the poles follow those of the GLAS 4th Order 
General Circulation Model (see Kalnay et al. , 1977) 


6. THK RF.SHLTS 

We performed three different model runs for our experiment. In each case 
a maxinum allowable time step was determined. A scheme was considered unstable 
when the height field decreased by one order of magnitude. The first run (Rl) 
consisted of updating the variables using the standard KO scheme. The second 
run (R?) replaced the KO scheme with the AG scheme in the continuity equation 
and In the pressure gradient terms in the momentum equation in order to Increase 
the time step defined by the fast moving Inertia gravity waves. The third run 
(R3) applied the AG scheme to the advective terms In the momentum equation as 
well. For all runs, the Initial condition was a Rossby-Haurwitz wave of zonal 
wavenumber 1 following Phillips (1959). The mean height of the wave was 9079 
m. Each run was Integrated for a forecast time of 120 hours. Table (5) shows 
the maxinum time step and the relative computing times for each run for which 
the scheme remained stable. 


TABLE 6 

Experimental Results (Course Grid) 


Model Run 

Maximum Time Step 

Rel. Comp, Time 

Rl 

1000 seconds 

1.000 

R2 

1350 seconds 

0.933 

R3 

1 

1400 seconds 

0.936 
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6. CONCLliniNG REMARKS 


It can he seen that a 35t to 40% increase In the maximun allowable time 
step was obtained throuqh application of the AG scheme. From the linear 
analysis In section (3), a qain of 20% could be expected using time steps cal- 
culated at ~60“. Neither scheme was able to give the theoretical maximum. 

The overall Improvement in the AG scheme found experimentally may be due to a 
better simulation of non-linear effects. 
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A MODEL TO DETERMINE OPEN OR CLOSED CELLULAR CONVECTION 
H. M. Helfand and E. Kalnay 


A siinple mechanism is proposed to explain the observed presence in the 
atmosphere of open or dosed cellular convection. If convection is produced by 
cooling concentrated near the top of the cloud layer, as In radiative cooling 
of stratus clouds. It develops strong descending currents which are compensated 
by weak ascent over most of the horizontal area, and closed cells result. Con- 
versely, heating concentrated near the bottom of a layer, as when an air mass 
is heated by warm water, results in strong ascending currents compensated by 
weak descent over most of the area, or open cells. This mechanism, unlike that 
of Hubert (1966), does not Invoke vertical variations of the eddy diffusion 
coefficients, and Is similar to the one suggested by Stommel (1962) to explain 
the smallness of the oceans' sinking regions. 

This mechanism Is studied numerically by means of a two-dimensional, non- 
linear Boussinesq model. An internal heat source-sink destabilizes a layer of 
fluid adding no net heating. A steady state Is attained. The horizontal con- 
figuration of the flow depends on the vertical distribution of the internal 
heating. Thus, when (a) Internal cooling Is concentrated near the top of the 
domain, (b) Internal heating is concentrated near the lower surface, or (c) the 
internal heating profile Is antisymmetric about the mid-plane of the domain 
(see Fig. 1), then the resulting flow pattern Is respectively (a) open, (b) 
closed, or (c) horizontally symmetric (see Fig. 2). The mean lapse rate is 
unstable in the upper half of the fluid and stable In Its lower half In the open 
case, stable In the upper half of the fluid and unstable in its lower half in 
the closed case, and unstable throughout in the symmetric case. 

The numerical results Indicate that the width of the plume produced by the 
cooling in the upper part of the layer or by the heating In the lower part of 
the layer Is largely Independent of the degree of vertical asymmetry of the 
internal heating profile (see Fig. 3). On the other hand, the compensating 
motion occupies a region which becomes broader as the heating profile becomes 
more asymmetric. In other words. If cooling Is very concentrated near the top 
of the layer with heating spread throughout the rest of the region or If heating 
Is very concentrated near the bottom with cooling spread throughout, the 
generated closed or open cells have an aspect ratio much smaller than 1, These 
results may help explain the small aspect ratio observed in atmospheric 
convection. 
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Fig. I 
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Fig. 2 


Fig. 3 


Fig. 1 Internal heating distribution (In non-dimensional units) as a function 
of non-dimensional height. The solid line represents the case of con- 
centrated cooling near the top. The dashed line represents the case 
of concentrated heating near the lower surface. The dashed-dotted 
line represents the case of heating antisymmetric about the mid-plane 
2 = 0.5. All of the profiles integrate to zero. The positive area 
under each curve (and the negative area, as well) Is the same. 


Fig. 2 The steady-state stream functions for the cases of (a) internal cool- 
ing concentrated near the top, (b) internal heating concentrated near 
the lower surface, and (c) antisymmetric internal heating. The cellular 
convective patterns are, respectively, (a) closed, (b) open, and (c) 
horizontally symmetric. 


Fig. 3 The horizontal profile of vertical velocity within the downdraft for 
the case of internal cooling concentrated at the top (solid line) and 
antisymmetric internal heating (dashed-dotted line). 


118 




25570 




LARGE AMPLITUDE STATIONARY ROSSBY WAVES IN THE SOUTHERN HEMISPHERE 

E. Kalnay and M* Halem 


1. INTRODUCTION 

The First GARP Global Experiment (FGGE) has provided the most extensive 
atmospheric data set ever compiled. An analysis of the atmospheric structure 
during the first Special Observing Period (SOP-1, January 5th, 1979 to March 
5th, 1979) was completed using the Goddard Laboratory for Atmospheric Sciences 
(GLAS) Analysis and Assimilation System. Preliminary results from this analysis 
were reported in Kalnay-Rivas et al. . 1980, where we pointed out the existence 
of large amplitude, short wavelength stationary waves in the lee of South 
America and Australia during the first 16 days of SOP-1. 

The purpose of this paper is to present further observational evidence of 
the characteristics of these waves during the complete SOP-1, and discuss their 
possible origin. The observed structure and evolution of the waves are presented 
in Section 2. Section 3 contains a brief discussion of the possible origin of 
the stationary waves. 


2. OBSERVATIONAL EVIDENCE OF STATIONARY ROSSBY WAVES 

One can expect significant differences in the structure and distribution 
of stationary waves in the Northern and Southern Hemispheres. In the Northern 
Hemisphere, the geographical distribution of orograhic and thermal forcing 
occurs on very large scales. Therefore it is not surprising that stationary 
waves have most of their energy in planetary wave numbers 1 to 3 (Charney and 
Eliassen, 1949; Smagorinsky, 1953). The Southern Hemisphere, on the other 
hand, is basically a homogeneous domain, with smaller isolated regions of land 
and sea contrast (South America, southern Africa and Australia) and one narrow 
major orographic chain (the Andes). No counterpart of the strong heating by 
the Gulf and Kuroshio western currents (Geller and Avery, 1978) exists in the 
Southern Hemisphere. 

Van Loon and Jerne (1972) and Van Loon et al. (1973) studied the standing 
(stationary) waves in the Southern Hemisphere and compared them with those of 
the Northern Hemisphere. They found that at 50®S, planetary wave numbers 1 and 
3 contain 99% of the variance of the annual mean, and contrary to their Northern 
Hemi‘jphere counterpart, are essentially barotropic. As in the Northern Hemi- 
sphere, they did not find significant stationary waves with synoptic scales. 

A very different picture emerges when one examines the FGGE data for the 
Special Observing Period. Figures la and lb present the 300 mb meridional 
velocity V averaged over two consecutive 30 day periods, Jan. 5- Feb. 3 and 
Feb. 4 -Mar. 5, 1979, respectively. In the Northern Hemisphere we observe for 
both months, as expected, the presence of planetary scale stationary waves. 
However, in the Southern Hemisphere, the January mean shows large amplitude 
synoptic scale waves immediately upstream and in the lee of the Andes, and 
weaker waves over Australia and in the South Pacific Ocean. In February, the 
synoptic stationary waves have disappeared and the circulation is more in 
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agreement with the results of van Loon et al. (1973), i.e. , very weak stationary 
waves of planetary acale. 

The amplitude of the January stationary waves In the South American region 
Is remarkable (15 m sec"' In the mean meridional wind), especially when one 
considers that this Is a summer, oceanic hemisphere. It Is also Interesting 
that the mean amplitude of transient waves (not shown here) Is about 50% stronger 
In February, when the stationary waves are not present. 

The January waves appear to have a virtually equivalent barotropic struc- 
ture. They are barely apparent at surface levels, with an anomalous trough In 
the subtropical South Atlantic center of high pressure. Their amplitude 
Increases with height, and they peak near the tropopause level. In the meri- 
dional direction, they occur In the 20® to 40® latitude band. 

The wavelength of the January stationary waves over South America Is about 
50® longitude, or zonal wavenumber 7 (Flo. 2a). Their position with respect 
to the Andes Is also anomalous. A very clrong ridge appears centered at a 
longitude of 70®W. As a result, the flow over the Andes Is, for the most part, 
from the North . On the other hand, the weak ridge upstream and trough downstream 
of the Andes that appears In February (Fig. 2b) Is typical of the winter circu- 
lation (Palmen and Newton, 1969). 

The orientation of the waves Is from the NW to SE in the Southern Hemisphere. 
This can be better observed In the map of absolute vortlcity (Figures 3a and 
3b) which also indicate the regions with significant convective precipitation 
according to the assimilation model. It should be noted In the January map 
that the the three major bands convective precipitation over the oceans (to the 
southeast of Indonesia and Brazil) coincide with regions of maximum advectlon 
of cyclonic vortlcity. The orientation of the waves implies that they transport 
westerly momentum Into the "roaring forties" latitudes. 


3. DISCUSSION 

We have shown the existence of synoptic scale stationary waves In the 
Southern Hemisphere for the month of January 1979. Confirming evidence comes 
from the examination of TIROS-N cloud pictures. Indicating the permanent pre- 
sence of a low level cyclonic vortex during the first two weeks of FGGE SOP-1 
at about 30®S, 40®W, in agreement with our analysis. The presence of the 
stationary waves over South America is also clear in the lower stratospheric 
temperature waves observed by Stanford and Short (1981). Other analyses, like 
the British Met. Office and ECMWF are also in agreement. That such waves are 
frequently present during periods of several weeks is very clear from the well 
known satellite-observed cloud and precipitation bands In the Southern Pacific 
and Atlantic Oceans (Rao et al. , 1976). 

The orientation and wavelength of these stationary Rossby waves appears 
similar to the numerical results obtained by Grose and Hoskins (1979) with a 
barotropic model forced by orography. The waves resemble In particular the 
solution forced by small scale topography at 30® latitude. In their results, 
higher latitude forcing generated more zonal ly oriented bands* 
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Several characteristics of the Rossby stationary waves should be explained: 
Why are they relati-.aly short? Why do they occur at subtropical latitudes? 

Why aren't they observed in longer time averages (seasonal to interannual). To 
answer these question: we may be guideu by several simple considerations: 

a) Tt.e scale of the forcing is very important. This may be seen in Figure 
4, which presents the solution of the 1-dimensional steady state barotropic 
vorticity equation on a 3-plane, forced by a "chapeau"-! ike mountain. When 

the mountain is of large scale (top panel) the response is primarily of the 
scale of the forcing, whereas when the forcing occurs in a small isolated 
region (bottom panel), the response occurs in the form of a train of stationary 
Rossby waves, in the lee o^ the mountain. This agrees well with the results 
obtained by Kalnay-Rivas and Merkine (1981). The first solution may be associ- 
ated with the Northern Hemisphere, where forcing, both thermal and orographic, 
occurs on planetary scales. The second solution is closer to the Southern 
Hemisphere, with forcing in relatively small isolated regions. 

b) Short stationary waves can occur in sub-tropical regions where the 
vertically-averaged zonal flow is positive but small as may be inferred from 
the frequency dispersion relationship of Rossby waves. The region where this 
condition is likely to hold is in the subtropical summer hemisphere, with weak 
easterlies at low levels and westerlies in upper levels. Because of their 
dependence on the zonal flow, the waves are not so apparent in longer term 
averages. 

c) The strong correlation between oceanic convective bands originating in 
Brazil and the Indonesian convergence zones, and the monthly averaged regions 
of cyclonic vorticity advection associated with the stationary waves, suggest 
that latent heat release may be the cause of their lar^ amplitude . It is 
possible that a CISK-type of instability occurs, in which the stationary Rossby 
waves determine the regions of low level convergence and ascending motion, 

and the resulting convective activity amplifies the Scutionary waves. Longer 
stationary waves occurring at higher latitudes where the vertically-averaged 
zonal wind is stronger, would not be as effective in producing convective 
precipitation. The fact that the waves appear to be equivalent barotropic also 
suggests that they are not a result of baroclinic instability. Sea surface 
temperature anomalies may also play a role if their phase is favorable for 
intensification of the waves through a CISK-mechanism. 

d) The large amplitude stationary flow over the Andes is poleward, rather 
than equatorward as required in order to generate perturbation (cinetic energy 
(Kalnay-ftivas and Merkine, 1981, and Figure 4). This also suggests that oro- 
graphy may not be generating the waves, but rather anchoring them at preferred 
regions. 

e) The extent of the tropical easterlies is smaller in January than in 
February 1979. In January, at 500 mb, the westerlies were present throughout 
most of the tropical Pacific, with speeds of up to 10 m sec"^ NE of Australia. 

In February there were easterlies throughout most of the 20“S to 15“N latitude 
band. The release of latent heat, in the Pacific convergence zone which accord- 
ing to the assimilation model is maximum at 500 mb, was stronger during January 
than during February. One may speculate that the reduced extent of easterlies 
at level of maximum release of latent heat was more favorable (Chamey, 1973, 

p. 293) for the propagation of waves originated In the Indonesian- South Pacific 
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convergence zone into the South American domain. However, it is not clear 
whether this had an effect on the generation of waves in the South American 
region. 


4. CONCLUSION 

We have printed out the existence of large amplitude, short wavelength 
stationary Rossby waves in the Southern Hemisphere during the month of January 
1979. They appear to be associated with observed stationary cloud bands in the 
South Atlantic and Pacific Oceans. Although we suggest that they attain their 
amplitude through i CISK mechanism, much further study is necessary to explain 
their origin and ?ir role in the general circulation of the Southern Hemi- 
sphere. 
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MULTIPLE EQUILIBRIA AND VACILLATION ON THE SPHERE. 
WITH APPLICATION TO BLOCKING* 

B. Legras^ and M. Ghil^ 


We study the barotropic vorticity equation on the sphere, including forcing 
and dissipation. A simplified orography, symmetric around the Equator, and a 
forcing which corresponds to a mid- latitude jet is used. The non-dimensional 
parameters of the problem are a Rossby number of the forcing Rq and a dissipa- 
tion coefficient o. 

The model is discretized in spherical harmonics, with a variable triangular 
truncation. Most of the study retains 25 real modes, while both lower and 
higher truncation is used in some experiments. 

Multiple equilibrium solutions obtain for the same value of (Rg»«); among 
these are low-index and split-jet, blocked flows. The dependence of equilibria 
on (Ro.a) is studied, along with their stability and bifurcation properties. 

Stable periodic, as well as aperiodic, vacillatory solutions are also 
obtained. The change in the character of solutions and in their stability is 
associated with foldings of the stationary solution sheets in parameter space. 

In the aperiodic case, solutions exhibit persistence periods longer than two 
weeks, as well as irregular, strongly oscillatory behavior. We also study the 
frequency spectrum of solutions as a function of (Rq,®). 

Both equilibria and vacillatory solutions seem to be structurally stable 
with respect to the number of modes used in the discretization. 


^ This work was presented at the 3rd Conference on Atmospheric and 
Oceanic Waves and Stability of the Amer. Meteor. Soc. , 19-23 January 
1981, San Diego, CA. It will also be presented at the lAMAP Symposium 
on the Dynamics of the Atmospheric General Circulation, 3-7 August 1981, 
Reading, England. 

2 laboratoire de Meteorologie Dynamique, Centre National de la Recherche 
Scientifique, 75231 Paris, Cedex 05, France. 

3 Courant Institute of Mathematical Sciences, New York University, New York, 
NY 10012, U.S.A. 
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THE ROLE OF ZONALLY ASYMMETRIC HEATING IN THE VERTICAL AND TEMPORAL STRUCTURE 
OF THE GLOBAL SCALE FLOW FIELDS DURING FGGE SOP-1^ 

J* Paegle^, E. Kalnay, and W. E. Baker 


1. INTRODUCTION 

The global scale structure of atmospheric flow is best documented on time 
scales longer than a few days. Theoretical and observational studies of ultra- 
long waves have emphasized forcing due to global scale variations of topography 
and surface heat flux, possibly interacting with baroclinically unstable or 
vertically refracting basic flows (Charney and Eliassen, 1949; Smagorinsky, 
1953; Bates, 1977; Grose and Hoskins, 1979, Holopainen, 1970; Lindzen and Tung, 
1979, and many others). Additionally there have been studies of the long term 
effects of tropical heating gradients (Ramage, 1968; Webster, 1972, Dickinson, 
1971; Rowntree, 1976; Julian and Chervin, 1978, Krishnamurti et al ., 1973, and 
others). 

Although these studies have produced a number of mechanistic explanattdns 
of period wave generation, it is probably fair to conclude that our understand- 
ing has remained insufficient to produce uniformly accurate deterministic pre- 
dictions beyond a few days for any space scales, including global scales. One 
reason for this may be that until the FGGE year, observational evidence about 
completely global scale short term variations was rather incomplete. Conse- 
quently, it has been difficult to verify theories of short term variations of 
global scale waves. Indeed, virtually all of the above mentioned theories that 
have quantitative form either emphasize local mid-latitude structure or long 
time (or time mean) interactions with the tropics. 

The present paper documents our analyses of SOP-1 data in terms of global 
scale spherical harmonics, with emphasis upon weekly transitions. 


2. APPROACH 

Level Il-b observed data have been analyzed (Baker, 1981) to generate a 
set of level Ill-b data on a 4® x 5® lat-lon grid. The 4th order GLAS GCM 
(Kalnay-Rivas et al. , 1977; Kalnay-Rivas et al. , 1979) served as the assimila- 
tion model. From this data set we have computed the vorticity (c) and the 
divergence (5) from 

c»ik.vxV (1) 

6=v.y (2) 

^ Published in Condensed Papers and Meeting Report: International Conference 

on Early Results of FGGE and Large-Scale Aspects of its Monsoon Experiments, 
Jan. 12-17, 1981; Tallahassee, Florida. Published by the World Meteorological 
Organization, Geneva, Switzerland. 

^On leave from the Dept, of Meteorology, University of Utah. 
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where V is the horizontal wind, Vis the finite difference gradient operator 
and IK is the unit vertical vector. The vorticity and divergence can be re- 
written in terms of a streamfunction (▼) and the velocity potential (X) 

(3) 

6= V^X (4) 

These equations have been solved for v and X employing a technique described 
by Paegle and Tomlinson (1975). The algorithm gives the exact solution of the 
algebraic system of equations arising from the finite difference approximation 
of the Poisson equation. 

Fields of v, X and geopotential height (Z) were subsequently spectrally 
decomposed as: 
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where X[J, Z[J are amplitudes of the spherical harmonics (P[J), of 
degree n, order m. 


3. INDIVIDUAL MODES 

Fig. 1 presents Z^ as a function of pressure level for a 2 uay period 
starting at OZ, January 28. ZS is the amplitude of a zonal harmonic possess- 
ing 3 nodal latitudes between the poles, but no longitudinal variation. In 
general, within th?* tropics, the low pressure belt at the surface reverses to 
a high pressure belt at high levels (because of the relative warmth of the 
tropical zone). In mid-latitudes, however, the zonally averaged pressure gener- 
ally decreases poleward at ar levels. The fact that Z^ maintains the same 
sign at all tropospheric levels suggests that the mid-latitude structure pro- 
jects more amplitude upon this zonal harmonic than does the vertically reversing 
tropical structure. 

A rather different conclusion holds for the Zj mode (Fig. 2). This 
model represents a field with three waves on each latitude circle and latitudinal 
nodes only at the poles. The complex amplitude of this mode reverses almost 
completely with height. This vertical variation is more symptomatic of monsoonal 
patterns (possessing low pressure at low levels and high pressure at high 
levels) than it is of typical mid-latitude waves (which generally display only 
mild vertical phase tilts or nearly equivalent barotropic vertical structure, 
et. Paegle et al. . 1979). 
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Figure 1. The amolitude of the zonal Figure 2. The complex ampli- 

height harmonic ZV. It is a real tude of the height harmonic 

quantity because this mode has no Z^. 

longitudinal variation. 



Figure_3. Complex amplitude of mode ¥3 of the strearafunction {left! apd of 
mode xj of the velocity potential (right). Each unit represents 10^ in‘/s. 
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Although an individual spherical harmonic of the height field has no direct 
physical relation to an individual spherical harmonic of the divergence agd 
vorticity fields, some of the characteristics of Zo are also evident in 
and Xj, (Fig. 3). Each of these also reverses with height within the ^ 
troposphere. Furthermore, both maximize in the upper troposphere where is 
about 50% larger than This means that the divergence is larger than the 
vorticity for this mode^at this level and time. 

These diagrams are representative of much of the SOP 1 data that we ha-'e 
examined. The planetary modes with increasing latitudinal structure tend to 
tilt only slightly with height in the troposphere, or to be nearly equivalent 
barotropic. Low order modes retaining east-west structure and latitudinal 
nodes only at the poles are strongly divergent and height revising. 


DIVERGENT U 


DIVERGENT V 



850 mb 



Figure 4. Divergent velocity components for the planetary scale modes, averaged 
over the 6 day period beginning OZ, January 28, and ending 18Z, February 2. The 
isoline interval is .5 m/s, and negative values are dashed. 
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4. HIGHLY TRUNCATED FIELDS 

Highly truncated sums of the rotational and divergent winds (obtained from 
truncated sums of f and X) exhibit an east-west structure that also tends to 
reverse with height. Figs 4 and 5 display 6 day average wind fields truncated 
by choosing N=4. Peak zonal winds to this truncation are about 30 m/s and are 
almost entirely rotational. Meridional winds, obtained by summing the rotational 
and divergent components, peak at about 6 m/s and the divergent component is 
about 50% of the rotational component. 

Thus, while the zonal flows are essentially rotational to this truncation, 
the distribution of meridional meanders is significantly impacted by the diver- 
gent flow field, especially in the upper troposphere. Additionally, the rota- 
tional wind field itself retains a markedly monsoonal character. In particular, 
the southerlies in the rotational meridional flow reverse to northerlies in 
height, and the maxima of the 200 mb zonal rotational flow are generally far 
removed from the 850 mb maxima. 


ROTATIONAL U 



ROTATIONAL V 




Figure 5. The same as Figure 4, for the rotational wind field. The isoline 
interval is 5 m/s for U and 1 m/s for V. 
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Figure 6. The same as Figure 4, for the second 6 day period, beginning OZ, 
February 3, ending 18Z, February 8. 


5. IMPLICATIONS 

Since we are describing the fields in terms of an arbitrary mathematical 
basis of expansion, possessing no a priori physical content, wo can only specu- 
late about the physical implications of these results. One implication of the 
presently described distribution is that tropical monsoons may be relevant to 
the longitudinal distribution of troughs and ridges, as well as jet streams 
and storm tracks (Krishnamurti , et al . , 1973). 

Similar hypotheses are also drawn in earlier studies by Djerknes (1969), 
j Ramage (1968) and Webster (1972), who have made comparisons of local Hadley 

I cell enhancement with subtropical jet streams on long (seasonal) time scales. 

The subject is currently receiving much attention in association with seasonal 
I abnormal ities. 

^ Much of this research suggests that an enhanced meridional overturning out 

of the tropical South Central Pacific, during the Northern Hemisphere winter, 
enhances the subtropical jet stream in similar longitudes. In the present 
study, for the sector from Central America, westward to the dateline, there is 
a rather enhanced poleward divergent flow in the first 6 day period relative to 
the second 6 day period (compare Figs 4, 5, with 6, 7). Our interpretation of 
this result is complicated by the fact that, to this truncation the zonal flow 
does not change substantially. 
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Figure 7. The same as Figure 6, for the second 6 day period, beginning OZ 
February 3, ending 18Z February 8. 



However, the potential relevance of tropical higher latitude interactions 
on weekly time scales cannot be dismissed, even here. The rotational flow 
couples the divergent flow to the extent that angular momentum is conserved, 

i.e. 


d(c ^ f) = -(f + (8) 

dt 

The appearance of the Coriolis parameter as well as nonlinearity, in the 
relation also tends to produce finer scale structures in the rotational flow 
response to changes of divergence, than are present in the divergence itself. 

In particular, the product of the Coriolis term with the divergence produces 
finer scales in che latitudinal structure of the evolving vorticity than are 
present in the divergence. 

At sufficiently high resolution, the discrepancy between the spectral 
truncation of the vorticity and of the divergence term of the vorticity equa- 
tion should not be important. Because of this, we have also decomposed the 
longitudinal variation of the stream- function and velocity potential in Fourier 
series truncated at wavenumber 4, and then computed the latitudinal amplitude 
structure of each of these modes to the full resolution of the 4® latitudinal 
grid spacing. 
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The 200 mb divergent meridional wind field that results at 18®N latitude 
is plotted as a function of longitude for each 6 day period in Fig. 8a. The 
strongest poleward flow (over 5 m/s) during the first period is found at about 
60®E over the Arabian Sea. This maximum diminishes markedly to about 3 m/s 
during the second 6 day period, while shifting eastward. A secondary maximum 
of about 3 m/s, centered at about 110®w (over the Eastern Pacific), during the 
first 6 days diminishes to a relative minimum of I m/s during the next 6 days 
as a new poleward miximum (emanating from the eastern Amazon Basin) strengthens 
over the Central Atlantic Ocean. 


The 200 mb rotational zonal flow structure at 30°N during the same two 
periods is displayed in Fig. 8b. At 90®W the 70 m/s zonal wind of the first 
period reduces to 43 m/s during the second 6 days, as a new maximum of 60 m/s 
establishes just off the East Coast of North America. A less dramatic weakening 
and eastward shift occurs in the jet stream around South-East Asia. These jet 
stream shifts appear to be consistent with the above-mentioned theories and 
regional shifts of the meridional overturnings. 




Figure 8a. Divergent meridional wind at 18°N, obtained as described in the 
text. The dashed curve represents the first 6 day average, the solid curve 
depicts the second six day average 

Figure 8b: The rotational zonal wind at 30°N, obtained as described in the 

text. The dashed curve represents the first 6 day average, the solid curve 
depicts the second 6 day average. 

The jet stream changes were accompanied by rather significant changes of 
weather over the United States. The cold wet weather over the South-Western 
states of the first 6 day period diminished during the second 6 days, as the 
storm track focused more strongly over the South-Eastern part of the country. 

In the instance of seasonal time scales, it appears that enhanced regional 
Hadley cell circulations are sometimes associated with warmer than normal tropi- 
cal sea surface temperatures and precipitation. It is possible that stronger 
local latent heating provides a regional energy source for the subtropical jet 
stream, thereby modifying subtropical jetstream meanders induced by other 
effects such as topography, baroclinity and wave interactions. 

Analagously, in the present case, it seems reasonable to suppose that 
regionally enhanced divergent meridional circulations over a period of a week 
are associated with locally enhance ITCZ or monsoonal activity. However, it is 
not possible to discern, without further analysis, whether mid-latitude effects 
produced meaners of the jet stream that induced stronger local tropical over- 
turnings or whether the latter may have had more independent origins. The 
question requires more detailed analysis, perhaps including numerical model 
experiments. One relevant comparison of such experiments to nature is to com- 
pare the model forecast global scale flows to their atmospheric counterparts. 
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INTEGRATED AND SPECTRAL ENERGY FLOWS OF THE GLAS GCM 

J* Tcnenbaum 


Analysis of the generation, transport, and dissipation of energy repre- 
sents one of the primary methods of studying geophysical fluid flows. Lorenz 
proposed subdividing atmospheric energy into kinetic and available potential 
forms and in turn subdividing these categories into zonal kinetic energy and 
eddy kinetic energy categories. Further subdivision into the spectral domain 
was suggested by Saltzman (1970) and followed by similar observational studies 
(Tenenbaum, 1976; Baker et al . , 1978). Our previous work in this area was 
summarized in last years annual report (Tenenbaum, 1980) and in a recent paper 
(Tenenbaum, 1981). 

During the past year these energetics analyses have been pursued in several 
directions. First, the longitudinal and time dependence on the energy flow to 
the stratosphere has been examined as a function of geographical sector. 

Second, the cases of strong and weak energy flows have been correlated with the 
skill of medium range forecasts. Third, the one-dimensional spectral results 
(Fourier services around latitude circles) have been extended to spherical har- 
monics over a global domain. Fourth, the validity of vertical velocities 
derived from mass convergence have been examined in terms of their effect on 
the conversion of eddy-available potential energy to eddy kinetic energy. 

One result in Tenenbaum (1981) is the existence of a weak vertical shear 
above the subtropical jet in all nominally 10-level general circulation models. 
When examined in more detail this phenomenon is correlated with an anomalous 
growth of kinetic energy in the GLAS model's stratospheric level (120 to 10 mb). 
This energy growth is closely related to the locations of strong initial jet 
streams and can be characterized geographically. For the North American and 
European sectors, the jet must be both quite strong and quite high for the 
anomalous energy rise to occur. Approximately one- third of the prominent jets 
in the initial states satisfy this criterion. In the East Asian sector, oro- 
graphy normally keeps the jets high. Thus 90% of the initial jets are anomalous. 

The forecast correlation studies dealt with whether the cases of signifi- 
cant energy rise were related to the downstream forecast skill. Two very sug- 
gestive results were found by examining the initial states of 14 forecasts 
made from FGGE data. The energy rise over North America was moderately corre- 
lated with the initial state divergence averaged over North America. The 
European 500 mb forecast skill at day 5 was even more correlated with this 
initial state divergence (Fig. 1). These results were independent of the size 
of the averaging region but possessed a very strong signal as one varied the 
location of the averaging region. 

For the two-dimensional spectra, preliminary analyses have been made at 
the end of the project year. The vertical velocity studies have been made 
using the geostrophic approximation to the diagnostic omega equation including 
both moist and dry terms. Initial results have shown unphysically smooth ver- 
tical velocity fields. We are currently examining whether the forcing terms 
in the omega equation (an elliptic partial differential equation) have the 
necessary scale of spatial variation to permit physically reasonable results. 
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Divergence (10~*per second) 






Fig. 2 


143 






9 


REFERENCES 

Baker, W. E., E. C. Kung, and R. C. J. Somerville, 1978: An energetics analysis 

of forecast experiments with the NCAR general circulation model. Mon. Wea. 
Rev. , 106 , 311-323. 

Saltzman, B. , 1970: Large scale ftmospheric energetics in the wavenumber 

domain. Rev. Geophys. Space Phys. . 8, 289-302. 

Tenenbaum, J., 1976: Spectral and spatial energetics of the GISS model atmo- 

sphere. Mon. Wea. Rev. , 104 , 15-30. 

Tenenbaum, J., 1980: Integrated and spatial energetics of the GLAS general 

circulation model. Atmospheric and Oceanographic Research Review - 1979 . 
Laboratory for Atmospheric Sciences, Nodding and^SimuVation facfTity. 

NASA TH 80650 . 

Tenenbaum, J. , 1981: Integrated and Spectral energetics of the GLAS general 

circulation model. (Submitted to Mon. Wea. Rev.). 


144 


55573 


A NUMERICAL MODEL STUDY OF LONG-TERM PLANETARY WAVE PREDICTABILITY 

R. L. Wobus 


A numerical model has been developed and used to study planetary wave 
behavior over periods up to four months under the influences of stochastic and 
diffusive forcinq by synoptic-scale motions and of sea surface temperature 
anomalies. The model is suitable for Monte Carlo experiments which provide 
ensemble mean and variance predictions and permit the response to heating ano- 
malies to be partly separated from random variability. 

The model incorporates the linear-balanced dynamical equations for wind at 
two tropospheric levels and temperature at one level. These variables are 
represented over the northern hemisphere by spherical harmonics truncated at 
global wavenumber 8. The model atmosphere is forced by semi -empirical January 
heating linearly dependent on local air temperature. Topography and the lower 
boundary layer are represented by vertical motion and drag functions, respec- 
tively, of the lower level wind. The linear phase speeds of Rossby waves have 
been modified to reflect the effect of the stratosphere. 

The truncation intentionally excludes synoptic disturbances and most baro- 
clinically active waves. The unpredictable effects of these motions are simu- 
lated by uniform random forcing; the predictable effects are simulated by 
vertically varying quasi -horizontal diffusion of heat and vorticity. The 
levels of random forcinq and horizontal diffusion are chosen by experiment to 
conform with the observed January mean circulation and the results of previous 
predictability studies* 

Random forcing is shown to be far more significant than variation of ini- 
tial conditions as a source of model unpredictability; in fact, under the 
influence of identical time-dependent random forcing, two model runs starting 
from unrelated initial states converge over a period of two months. Random 
forcinq of the largest-scale mode produces episodes of strong baroclinic wave 
growth, somewhat analogous to zonal index variations in the atmosphere. An 
ensemble of 8 parallel 120-day runs from a common initial state shows the 
model's climate, variability, and short-range predictability. 

In order to separate the relatively small responses to a North Pacific sea 
surface temperature anomaly (similar to that used in NCAR GCM studies) from the 
variability due to unpredictable synoptic-scale motions, the 120-day ensemble 
is rerun with the anomaly but with the same 8 sequences of random forcing. The 
ensemble of differences between corresponding runs with and without the anomaly 
provides a sensitive measure of the response. The ensemble mean response grows 
raoidly during the first two weeks and more slowly thereafter; the ensemble 
variance of the differences grows rapidly at first, then more slowly until day 
60, overwhelming all but the local mean response after day 30. This secondary 
variance growth occurs primarily through episodes of strong baroclinic wave 
growth. 

These experiments suggest, for the large-scale flow in winter, a time- 
spectral gap between the period over which detailed forecasts can be made from 
initial states (about 2 weeks) and that over which large-scale statistics 
follow seasonal and anomalous variations of thermal forcing (two months or 
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lonyor); In the range from two weeks to two months the uncertainty propagated 
into the planetary flow from synoptic- scale motions should significantly limit 
the skill of forecasts* even If the external forcing were completely known. 
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CLIMATE MODEL POSTPROCESSING 
J. Abeles, E. Pittarelli and D. A. Randall 


Enormous effort and expense is devoted at GLAS and at other centers through 
out the world to design, code, test, and perform simulations with global general 
circulation models, the effort and expense are justified, since these models 
represent our best hope to improve our understanding of and ability to predict 
the vagaries of weather and climate. But in order to learn from and utilize the 
model results, it is necessary to create elaborate "postprocessing" software to 
allow analysis of the large volume of data produced. The model results would 
be practically useless without postprocessing software. This report briefly 
summarizes the development of new postprocessing software for the GMSF Climate 
Modeling Group. 

To analyze the results of the potential -enstrophy-conserving general cir- 
culation model (PEGCM) and the seasonal cycle model (SCM), extensive new analy- 
sis codes have been developed. 

The models produce sigma history tapes. The sigma history records are 
interpolated to pressure history records, which are written on a pressure his- 
tory tape. This allows us to analyze the model results on pressure-surfaces, 
with snap-shots or time averages. 

The nominal pressure surfaces used for the pressure history tape are: 100, 

200, 300, 400, 500, 600, 700, 800, 900 and 1000 mb. For each of these pressure 
surfaces, a mask is created. If the pressure surface is below ground, the mask 
is zero, and if the pressure surface is above ground, the mask is one. The 
mask is utilized in the time averaging and plotting software. 

The following diagnostic fields are written to the pressure history tape: 


1 ) u -wi nd 

2) v-wind 

3) potential temperature 

4) temperature 

5) relative humidity 

6) specific humidity 

7) geopotential 


8) vorticity 

9) omega 

10) omega * alpha 

11) stream function 

12) velocity potential 

13) divergence 

14) sea level pressure 


The interpolation of the various fields is done by utilizing tailored 
methods: 

1) Potential temperature and relative humidity are interpolated linearly in 
(p/pop (Specific humidity is converted to relative humidity, and the relative 
humidity is interpolated. Specific humidity on the pressure surface is then 
calculated from the interpolated relative humidity.) The temperature on the 
pressure surface is calculated from the interpolated potential temperature. 

2) The geopotential and u and v winds are interpolated using methods sug- 
gested by Prof. M. Suarez of UCLA. 


149 


3) Omeqa and omeqa * alpha are Interpolated linearly in pressure. 

4) Vorticity is computed directly on the pressure surfaces, usinq the 
interpolated u and v wind components. 

5) Divergence on the pressure surface is calculated by differencing omegas 
in the straddling sigma layers. 

6) Stream function and velocity potential are computed from the vorticity 
and divergence, using codes provided by Prof. J. Paegle of the University of 
Utah. 


In addition, the sea level pressure is calculated. 

In order to analyze time means of these fields, we time average the pres- 
sure history records. The following fields are put to disk: 


1 ) 

2 ) 

3) 

4) 

5) 

6 ) 

7) 

8 ) 
9) 

10 ) 


^otal mask 

£ 

V 


T 

W 

I 


11) t2 

12) Sh2 

13) 52 

14) UV 

15) V5FT 

16) 7T 

17) U5H 

18) uT 

19) Tort 

20) Omega 


In addition, the following total, transient and stationary variances and 
covariances are calculated: 


1) KE 

2) Temp 

3) SH 

4) ^ 

5) UV 


6) VT 

7) UT 

8) VSH 

9) USH 


The following is a list of some of the plots currently being produced 
by the climate postprocessing package, utilizing the NCAR Graphics software: 

For each of the following, world map contour plots and zonal means for 
total, ocean and land are produced. 


(ieopotential Height* 
Vorticity* 

Temperature (K)* 
Relative Humidity* 

Sea Level Pressure 
Ground Temperature (C) 


♦Plotted at 200, 500, 700 and 800 mb 


Surface Temperature (C) 

Total Precipitation (imn/day) 

Convective Precipitation (mm/day) 

Large Scale Precipitation (mm/day 
Evaporation Flux at the Surface (mm/day) 

Latent Heat Flux at the Surface (watts/m^) 

Sensible Heat Flux at the Surface and 
at the Top of the Atmosphere (watts/m^) 

Long Wave Radiation at the Surface and 
at the Top of the Atmosphere (watts/m^) 

Solar Radiation Absorbed by the Atmosphere (watts/m^) 
Planetary Albedo 

Radiation Balance at the Surface and at the Top of the 
Atmosphere (watts/m^) 

Heat Storage in Water or Ice at the Surface (watts/m^) 
Energy Balance at the Surface (watts/m^) 


The following are latitude-height plots: 


Zonal average U-wind (m/s) 

Zonal average V-wind (m/s) 

Zonal average Temperature (K) 

Zonal average Temperature (K) and U-wind (m/s) 
Zonal average Specific Humidity (g/kg) 

Zonal average Relative Humidity (g/kg) 


All of these fields are plotted for total eddy, transient eddy and stationary 
eddy: 


Latitude Height Plots 


Kinetic Energy (m^/s^) 

Temperature Variance 

Specific Humidity Variance (g/kg) 

Geopotential Variance 

Meridional flux of Zonal Momentum (m^/s^) 

Meridional Temperature flux (m/s) 

Zonal Temperature flux (m/s) 

Meridional Mixing Ratio flux 
Zonal Mixing Ratio Flux 


World Map Streamline Plots 


Winds 

Moisture flux 
Transient Temperture flux 
Momentum flux 
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Postprocessing software development is an ongoing project. Code has also 
been developed to calculate and analyze moisture flux convergence, total dia- 
batic heating, total cloudiness, heating due to dry convective adjustment, and 
cumulus heating and drying. 
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SEASONAL CYCLE OF HEIGHTS, TEMPERATURES, AND SEA LEVEL PRESSURE 
OVER THE NORTHERN HEMISPHERE 

B. Doty and K. Mo 


1. INTRODUCTION 

Recently Hsu and Wallace (1976) and White and Wallace (1978) documented 
the annual and semiannual cycles of sea level pressure and surface temperature 
over the Northern Hemisphere, using data based on NMC analyses for the Northern 
Hemisphere, and data described by Taljard et al « (1969) for the Southern Hemi- 
sphere. Hsu and Wallace (1976) also examined the annual cycles in the geopo- 
tential height fields for several levels. They found that the continental 
winter maxima in the sea level pressure field extended upward to near the 700 mb 
level. However, at 500 mb and above, they found that the annual cycle exhibited 
summertime maxima over both sea and land. 

Van Loon and Jenne (1970) studied the annual cycle of 100 mb temperature 
using the same data sources. They found highs In the amplitude between 30 to 
45 degrees North, the maximum occurring in the late winter or early spring. 

Reed and Vicek (1969) carried out harmonic analysis of monthly mean tempera- 
tures over the tropics to obtain the amplitudes and phases of the annual and 
semiannual cycles for 17 levels between 300 mb and 100 mb. They found that the 
annual temperature variation in the lower tropical stratosphere and its maximum 
amplitude lust above the equatorial tropopause had nearly the same phase at 
that particular level throughout both hemispheres. They also discovered that 
the coldest stratospheric temperatures in the tropics occurred during the 
northern hemisphere winter. 

Most analyses of annual and semiannual cycles have been done using harmonic 
analysis. However, Heddinqhaus and Kung (1980) recently used eigenvector analy- 
sis to study annual and semiannual cycles of 700 mb and 500 mb temperatures, 700- 
500 mb thickness, etc. over the Northern Hemisphere using NCAR data for the 
period 1955 to 1974. 

The annual and semiannual cycles over the Southern Hemisphere have been 
studied in great detail by Van Loon et al (1972) for temperatures, heights, and 
winds at all levels from 100 mb to sea level. 


2. RESULTS 

Data used in the analysis was obtained from NCAR, consisting of 15 years 
of data for heights at all standard levels, 14 years of temperatures at all 
levels, and 14 years of sea level pressure. Data was given twice daily from 
Jan 1, 1963 for heights (Jan. 1, 1964 for temperatures and sea level pressure) 
to December 31, 1977 on the NMC octagonal analysis grid. The data was inter- 
polated to a 4x5 degree latlon grid, and missing and erroneous data was replaced 
by linear Interpolation using neighboring dates (Doty and Mo, 1981). 

Harmonic analysis was performed at each grid point for each variable at 
each level using the twice daily time series: 10958 (15*365*2+8) time steps 

for heights and 10228 (14*365*2+8) time steps for temperatures and sea level 
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pressure (leap years included). The annual cycle is defined as the 15th (or 
14th) Fourier coefficient giving a period of 365.3 days; the semiannual cycle 
is defined as the 30th (or 28th) Fourier coefficient giving a period of 182.6 
days. 


Calculations of the standard deviations of the annual cycle amplitudes for 
the 15 (or 14) years were also made. This was done by calculating the annual 
cycle for each year separately at each grid point. 

Results are presented in a vectorial format, with each vector representing 
the amplitude and phase, the phase being the date when the harmonic is at its 
maximum. The vector convection is for the length of the vector to indicate the 
amplitude (the scales are not uniform) and for the direction to indicate the 
phase. An arrow pointing south to north indicates a maximum on Jan, 1 for the 
annual cycle, and maxima on Jan. 1 and July 1 for the semiannual cycle. The 
dates become later as the arrow is rotated clockwise, so that an arrow pointing 
west to east indicates a maximum on Apr. 1 for the annual cycle and Feb. 15 and 
Aug. 16 for the semiannual cycle, etc. 

Figure 1 shows the vector diagrams of the annual cycle for 200 mb and 700 
mb temperatures. Figure 2 is the same as Figure 1 but for 300 mb height and 
sea level pressure. Results based on calculations at other levels are discussed 
below, but are not shown. 


3. DISCUSSION 

The phase of the sea level pressure annual cycle shows wintertime maxima 
over the continents with summertime maxima over the oceans. The continental 
wintertime maxima extends to the 850 mb level in the annual cycle of heights, 
especially over Asia. At 700 mb, except for a small area over Asia where weak 
amplitudes prevail, the maxima occur in summer. 

The semiannual cycle in the height fields is stronger in the lower levels 
when expressed as a percentage of the strength of the annual cycle. At 850 mb 
the amplitude of the semiannual cycle is in places more than half the strength 
of the annual cycle. At higher levels, 500 mb and above, the semiannual cycle 
rarely exceeds one third the strength of the annual cycle. 

The amplitudes of the annual cycle of temperatures from 350 mb have very 
similar patterns. Highs occur over North America and eastern Asia, with axes 
of low amplitude occurring over the easternmost portions of the oceans. The 
pattern changes markedly at 200 mb, where the amplitude is largest over the 
polar region. The standard deviations between years for the annual cycle are 
very noisy and are generally small (less than 1C). 

The maxima of the annual cycle temperatures generally occurs in summer. 

At 200 mb the phase of the maxima is about one month earlier over the con- 
tinents than over the oceans. The phase shifts sharply eastward from ocean 
to land over western North America; the shift is more gradual from ocean to 
land over Eurasia. The latest dates for the maxima occur generally over the 
eastern portions of the oceans. At higher levels the phase pattern is simi- 
lar, with the ocean/land difference less well defined. 
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The patterns of semiannual cycle amplitudes of temperatures are similar 
from 850 mb to 300 mb, with highs over the pole, over western North America, 
and over the western Mediterranean. These amplitudes all have maxima In 
January or February and June or July. There Is also a small relative high 
in the amplitudes near India, with the maximum occurring In April and October. 
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ANALYSIS OF THE SURFACE HEAT BALANCE 
OVER THE WORLD OCEAN 


S. K. Esbensen 


CoMoboratlon with Professor Y.-J. Han and graduate students S.-W. Lee and 
Y. Kushnir has resulted In new estimates of the net surface heat fluxes over 
the global ocean for all calendar months (Han et a1. » 1980). The method used 
In the present calculation Is patterned after Budyko (1963) and Bunker (1976). 

But further simplifications were made, following Haney (1971), to obtain a 
formula In the form Qj ■ 02 (T*a “ where Qs Is the net surface heat flux, 

Tj Is the sea surface temperature, Is the "apparent" atmospherlj^ equili- 
brium temperature, and Qo Is the proportionality constant. Here T^ and 03 , 
derived from the original heat flux formulas, are functions of the surface 
meteorological parameters (e.g., surface wind speed, air temperature, dew 
point, etc.) and the surface radiation parameters. This particular formulation 
of the net surface heat flux together with the necessary climatological atmo- 
spheric parameters, are Intended to provide a realistic and computationally 
efficient upper boundary condition for oceanic climate modeling. 

The basic climatological data set used to compute T^ and Q 2 was obtained 
from several sources. The surface atmospheric parameters were provided on a 
5* X 5® Marsden subsquare grid by the National Climatic Center. The surface 
Insolation under the clear sky condition was taken from the table In Kondratyev 
(1969), but reduced by 10 percent at all latitudes to be consistent with more 
recent estimates (Bunker, 1976). The total cloud cover needed for both long- 
and short-wave radiation calculations was obtained from Miller and Feddes (1971), 
who estimated the monthly mean cloud cover from the satellite measurements. 

The monthly mean values of T^ and Q 2 were computed over the global ocean 
for all calendar months, but only the annual mean values are shown in Fig., 1 
and Fig. 2 to show the significant geographical dependence of these parameters. 
Although values are calculated In polar regions, the density in both the Arctic 
and Antarctic Oceans is not^large enough to provide stable climatological means. 
As a result, the computed T^ and Q 2 fields near the climatological sea-ice 
boundaries are noisy even In the annual mean maps. 

The annual mean net surface heat fluxes (Fig. 3) were computed and comiiared 
with the earlier estimates by Budyko (1963) for the world ocean ar'i Bunker 
(1976) for the North Atlantic Ocean. The flux patterns were in good agreement 
at middle latitudes but some significant discrepancies were found at high and 
low latitudes. These discrepancies are also apparent In the zonal mean values 
shown in Fig. 4, where more recent estimates by Oort and Vonder Haar (1976) 
are also plotted for comparison. A more detailed comparison of heat flux 
component shows that the major discrepancies are due to the estimates of the 
net radiation In the tropics and of the sensible heat fluxes In the high lati- 
tudes. He suspect that the solar Insolation Is somewhat overestimated In the 
present calculation mainly due to the use of the satellite-inferred clot') data 
which In the tropics, show systematically less clouo cover than the conventional 
ground observations (Fig, 5); the latter were used by both Budyko and Bunker. 

The relatively large heat fluxes In the high latitudes are probably due to the 
use of the heat transfer coefficient which are made to vary «<th the thermal 
stability and wind speed (Bunker, 1976). 


In order to investigate the seasonal variability of the net flux, the 
nonthly zonal mean values were computed and are shown in Fig. 6. The largest 
seasonal variabilities occur in the middle latitude oceans in both hemispheres. 
The relative magnitude of the variability, however, is greater in the northern 
oceans than in the southern oceans. This is consistent with the relatively 
stronger surface heat exchange processes in the northern ocean, especially in 
the vicinity of the major western boundary currents. For the purpose of com- 
parison, the present result was checked with the independent calculations of 
Oort and Vonder Haar (1976). Their estimates of the net surface heat flux 
which was obtained as a residual of the vertically integrated atmospheric heat 
balance equation using satellite radiation measurements at the top of the 
atmosphere and conventional radiosonde data shows relatively less heating 
(~50 watts/m^) during the summer half year and very much less cooling at the 
high latitudes during the winter. The relative accuracy of these two indepen- 
dent estimates, however, is difficult to judge because of many uncertainties 
involved in the methods and data used in each estimate. 

The parameterization of the net surface heat flux by T^, T^, and Q 2 
emphasizes the important role of sea surface temperature in determining the 
net surface heating of the oceans. The comparison between the heat fluxes 
obtained from model -simulated and observed sea-surface temperatures for the 
same T« and should give needed guidance for the development and verifica- 
tion of coupled ocean-atmosphere general circulation models. 
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SPATIAL AND TEMPORAL SCALES OF TIME-AVERAGED 700 MB HEIGHT ANOMALIES 


D. Gutzler 


1. INTRODUCTION 

The monthly and seasonal forecasting technique currently in use by the 
National Meteorological Center (NMC) is based to a large extent on the extrapo- 
lation of trends in the positions of the centers of time-averaged geopotential 
height anomalies. The complete forecasted height pattern is subsequently 
drawn around the forecasted anomaly centers. 

To test the efficacy of this technique, we have examined time series of 
observed monthly mean and 5-day mean 700 mb geopotential heights. Autocorrela- 
tion statistics were generated to document the tendency for persistence of ano- 
malies. These statistics were compared to a “red noise" hypothesis to check for 
evidence of possible preferred time scales of persistence. Space-time spectral 
analyses at middle latitudes were checked for evidence of periodicities which 
could be associated with predictable month-to-month trends. A local measure of 
the average spatial scale of anomalies was devised for guidance in the con^le- 
tion of the anomaly pattern around the forecasted centers. 


e. SPATIAL SCALE 

The average spatial scale of 700 mb height anomalies was determined by con- 
sidering the area surrounding each gridpoint over which contemporaneous heights 
were highly correlated. To estimate the tendency for anomalies to cover dif- 
ferent areas over different geographical regions, the complete set of spatial 
cross-correlation maps for both monthly and 5-day means was computed, following 
the procedure used by Wallace and Gutzler (1981). Each map shows the contempor- 
aneous correlation coefficient between heights at a specified base gridpoint and 
all other qridooints. As discussed by Wallace and Gutzler, the size of the 
region of high correlation surrounding each base gridpoint is an indication of 
the average horizontal extent of height anomalies in the vicinity of that grid- 
point. To quantify this extent, the approximate area within the 0.7 contour on 
each cross-correlation map was computed by summing the areas associated with 
the 4® by 5® patch centered on each gridpoint surrounding the base gridpoint 
which was correlated with a coefficient greater than 0.7. 

Figure 1 shows the geographical distribution of the horizontal extent sta- 
tistic for wintertime monthly means plotted at each gridpoint. Southward of 
the thick dashed line, the calculated area of high correlation extends to the 
boundary of the data set, and therefore is probably an underestimate of the 
actual area. 

Anomalies tend to be largest over the central Pacific Ocean, central Asia, 
northern Siberia, and Greenland. Secondary maxima in Figure 1 exist off the 
southeastern coast of the United States, over the northern Atlantic Ocean, and 
south of Japan. Regions where anomalies tend to be small include the Gulf of 
Alaska, the central United States, eastern China, and the North Sea. 
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3. PERSISTENCE 


The autocorrelation of monthly mean data (not shown) was found to be quite 
low northward of 30° N, 

The geographical distribution of autocorrelation for the 5-day mean time 
series with one lag» based on all months, is shown in Figure 2. Maxima in the 
autocorrelation field appear around the edge of the data set near the tropics, 
over the north central oceans, and over the polar region north of the Soviet 
Union. 


Klein (1951) postulated that height variability was well simulated by a 
first-order autoregressive or "red noise" model, in which anomalies relax expo- 
nentially toward a climatological mean value. If there exists a preferred time 
scale for atmospheric persistence, as suggested by Shukla and Mo (1981), then 
this red noise hypothesis should not serve as a good predictor of the observed 
5-day mean autocorrelation over regions of preferred persistence. The hypo- 
thesis was tested by comparing the observed 5-day mean autocorrelation with 
hypothetical values extrapolated from the observed daily autocorrelation at 
each gridpoint. Observed and hypothetical autocorrelation values generally 
agreed to within 30%, but distinct geographical variations could be seen. It 
was found that anomalies tend to relax toward climatology somewhat faster than 
red noise over the eastern oceans, western continents, and polar regions (par- 
ticularly Siberia and northern Canada). Anomalies persist longer than a red 
noise model would predict over the subtropics and western oceans. The pattern 
of discrepancies closely resembles the band-pass variance field calculated by 
Blackmon (1976), which included fluctuations with periods between 2.5 and 6 
days. Thus, the observed 5-day mean autocorrelation is larger than expected 
over regions associated with strong baroclinic wave development. 

In summary, significant autocorrelation exists in the time series of 700 mb 
height 5-day means. This autocorrelation could be explained in terms of two com- 
ponents. First, long-period fluctuations account for much of the tendency for 
persistence over the entire hemisphere. These fluctuations comprise the vari- 
ance of monthly means, and show up in the daily autocorrelation field. The 
large amount of this low frequency variance explains the reasonably good agree- 
ment with Klein's red noise hypothesis. Second, baroclinic disturbances create 
a significant component of variance at time scales too short to be reflected in 
daily autocorrelation statistics, but which enhance the 5-day mean autocorrela- 
tions over the storm track regions off the east coasts of the continents. Auto- 
correlation did not explain a significant amount of the variance of monthly means. 


4. PERIODICITY 

Periodicity in the 700 mb height time series was studied using a space- 
time spectral analysis. In this technique, the total space-time variance of 
heights along a specified latitude circle is decomposed into separate frequency 
and zonal wavenumber components using a Fast Fourier Transform (FFT), The 
routine used here is described in more detail by Straus and Shukla (1981). 
Height anomalies were used for the analysis, so that the stationary component 
of each wave was removed, and only the transient component, minus the annual 
cycle, was left. 
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The space-time spectral analysis of monthly mean 700 mb height anomalies 
at latitude 60°N is plotted in Figure 3. The wavenumber spectra (the columns 
of Figure 5) were smoothed by averaging over six discrete spectral estimates, 
so that each entry in Figure 3 contains twelve degrees of freedom. Isolines of 
variance per unit frequency interval (units nr*months) were drawn on the basis 
of the smoothed wavenumber spectra assuming that the variance was a slowly 
varying function of both frequency and zonal wavenumber. 

The significance of maxima in the space-time analysis can be estimated 
from a null hypothesis that the variance is independent of both frequency and 
wavenumber. The total variance would then be evenly divided among the 286 
spectral extimates included in Figure 3, This amounts to 340 units per estimate. 
Any observed spectral estimates greater than 450 units are then significant at 
the 99% level. It should be kept in mind that, in the absence of an a priori 
expectation that peaks would appear at specific frequencies and wavenumbers, 

1% of the estimates are expected by random chance to exceed the 99% significance 
level , 


The spectra for wavenumbers 1-3 are almost entirely above this level, and 
nearly all the variance (95%) is explained by fluctuations of the zonal mean 
and wavenumbers 1-4, This re-emphasizes the dominance of large-scale anomalies, 

A slight tendency for redness is demonstrated by the larger variances 
explained by lower frequencies in all wavenumbers. This is representative of 
the small (but generally non-zero) magnitude of autocorrelation of monthly 
means shown previously. The spectral peaks observed in higher frequencies are 
probably not significant because of the aliasing problems inherent in a time 
series of monthly means, each of which is itself a finite sample of observations 
containing a statistical uncertainty (Trenberth, 1980), 


5. CONCLUSION 

Nap et al , (1981) verified several operational forecasting schemes (includ- 
ing NMC's), and found that they all possessed little skill. The slight persis- 
tence and lack of dominant periodicities in the monthly mean time series indicate 
that the potential for improvement in the skill of trend-based long-range fore- 
casting techniques is not great. However, recent modeling experiments by Shukla 
(1981) suggest that the atmosphere is dynamically predictable for up to one 
month. The results here suggest that progress in long-range forecasting will 
come through new approaches to the problem, rather than more detailed analyses 
of longer time series of the same set of atmospheric variables. 
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Figure 1. Geographical distribution of the area within the 0.7 contour on the 
spatial cross-correlation map for each gridpoint, based on a 15 winter 
sample of monthly mean height anomalies for the months of December, January, 
and February. Contour interval 100 x 10^ km^. Southward of ♦’he heavy 
dashed line, the 0.7 correlation contour intersects the southern boundary 
of the data set, so that the calculated area may be underestimated. 
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Figure 3. Variance of monthly mean 700 mb height at bO^N as a function of 
frequency and wavenumber for the period 1953-1978. Contours represent 
isolines of variance per unit frequency interval (units m2*months). The 
450 unit contour is the 99% significant level, testing for departure from 
white noise. Light stippling indicates spectral peaks with magnitudes 
greater than 2000 units. The tick marks on the vertical axis denote the 
midpoint of the frequency bands, and indicate the degree of frequency 
resolution used. 
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HID-LATITUOE WIND FORCED OCEAN CIRCULATION STUDIES 
D.E. Harrison 


A simple barotropic vorticity equation model has been developed to study 
some of the various modeling factors that affect the characteristics of strong 
western boundary currents like the Gulf Stream and Kuroshlo. Successful pre- 
diction of sea surface temperature, both In the climatological mean and over 
periods as short as one month requires that the heating tendency, due to 
horizontal advectlon of heat by these currerts, be accurately modeled. Con- 
ventional, coarse resolution ocean models do not satisfactorily reproduce 
the dominant features of these currents. It Is important to understand why 
they do not and what must be done In order that they will be able to do so 
in the future. 

The first set of calculations, making use of a simple linear bottom drag 
law as used by Veronis (1966), has been completed and submitted for publica- 
tion (Harrison and Stalos 1981). This work establishes the vital Importance 
of understanding the mechanisms of vortldty loss In the western boundary 
currents and the necessity of studying adjoining circulation systems together 
rather than separately. The classical wind driven ocean studies, which focus 
on a single subtropical gyre driven by a anticyclonic wind stress pattern, 
do not allow for the poss1b111':y of vortldty transport by the circulation, 
out of the gyre. When this process u permitted by the model boundary condi- 
tions, geometry and forcing, drvamatically different-circulation patterns 
result (Figures 1 and 2). 

Apart from these qualitative results. It Is also Important to understand 
quantitatively how the model "Gulf Streams" depend on the modeling parameters. 
A simple theoretical argument has been developed that predicts that the speed 
of the eastward flowing "separated Gulf Stream" will vary as: 

Uj - Lx curl^T 

IT — S — 

where Lx Is the east-west basin dimension, K Is the linear drag Inverse decay 
time, curljT Is the maximum v.-ind stress curl and D Is the mean depth of the 
fluid layer. The width pf this current Is predicted to vary like: 

Lj ~ Uj 

wonm 

where 6 Is the variation of Coriolis parameter with latitude, so long as 

Uj < U* = KLx 

When Uj » Ux, the classic Inertial scaling 

Lj - Uj 1/2 
6 
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should be observed. These results are confirmed by numerical experiments 
(Figures 3, 4). To our knowledge this is the first non-linear and viscons 
system whose quantitative behavior has been understood from first principles. 

Further work is underway to understand how these results generalize when 
other frictional mechanisms are used, and baroclinic studies will be initiated 
after the barotropic calculations are complete. 


REFERENCES 

Harrison, D.E., and S. Stalos, 1981: On the wind-driven ocean circulation, 

(submitted) 

Veronis, G., 1966: Wind-driven ocean circulation. Part II Deep-Sea Res. 

13, 31-55. 


3 


171 


Figure 1 


ORIG'NAL PAGE IS 
OF POOR QUALITY 


NORnRLIZECI BI 9.8S3tMc’ 


HOOEL B 
R = 2.000E-02 
E = 5.000E-0? 

















SEA SURFACE DATA STUDIES 


D.E* Harrison 


The HSST form of the TDF-11 data set from the National Climate Center Is 
being used to produce a high sp&tlal and temporal resolution climatology for 
sea surface temperature, air-sea temperature difference kinetic energy Input 
to the mixed layer, wind and wind stress over the world ocean. High spatial 
resolution, as close to 1° x 1° resolution as the data will usefully permit. 

Is needed for mar\y ocean modeling purposes because there are many special 
oceanic regions (primarily along or adjacent to coasts and the Equator) where 
there are strong gradients across a few hundred km and because It Is oHen 
necessary to know the curl of the wind stress In order to know the dominant 
ocean forcing. Estimates of wind stress curl from spatially smooth data may 
seriously underestimate the magnitude of the curl. Considerable attention Is 
also being given to making estimates of the error bars on the various climato- 
logical quantities, as error estimates are necessary to establish the climato- 
logical "noise" that will affect climate anomaly studies. 

The Atlantic Ocean climatology calculations using data between 1850 and 
1975 have been made, and detailed analysis of these results Is now underway. 
Figure 1 shows SST results for January, off the Africlan coast (10®W*30°W. 
20“N*30"N). Isotherms have been sketched In; the data In each 1® box are: 
number of observations, average value and standard deviation. Fig. 1 also 
Illustrates several features characteristic of the data set: the number of 

observations per 1® square per unit time varies greatly within a 10® square 
region (from 1 observation to 1331 observations In this case), obviously 
spurious values can be found In the 1® means when the number of observations is 
very small (here < 5obs.) but overall the data vary quite consistently on the 
1® scale without any smoothing, and there Is strong spatial variation on the 
1“ scale that would be seriously smoothed out If the data were averaged over, 
say. 5®, squares as has been done In previous empirical orthogonal function 
analyses of the Atlantic. 

After the climatological results are well understood and documented, 
anomaly studies will be carried out, to try to better document the variations 
from climatology that have occurred over the last 130 years in the Atlantic. 

The Indian and Pacific Oceans will then be studied In a similar fashion. 























GLOBAL FIELDS OF SOIL MOISTURE AND LAND-SURFACE EVAPOTRANSPIRATION 


Y. Mintz and Y. Seraflnl 


Zeroth-order estimates of the monthly normal soil 
surface evapotranspiratlon are obtained by Integration 


moisture and land- 
of the equation » 


dW * P - 
It 


E. 


•^max “ 


W*. 


E - B Ep. 

where U Is the evapotranspi ration-avail able soil moisture and E Is evapotrans- 
piratlon. P Is precipitation, taken from the observed monthly normals given 
by Jaeqer (1976). Ep Is potential evapotranspiratlon, calculated from the 
observed monthly normal surface air temperature using the empirical formula of 
Thornthwalte (1948). M* Is the maxlnum available moisture that the soil can 
hold, and Is taken as B, the evapotranspiratlon coefficient, is 

taken as S = 1 exp 1 when P < Ep, and 3 = 1 when P > Ep. 

Starting from an Initial state in which W = 0, the governing equation Is 
Integrated forward In time until a steadily seasonally varying state Is reached. 
The calculations are made at Intervals of of longitude and 4** of latitude 
over the globe. The output data, U and E, as well as the Input data, P and 
Ep(T), are stored at half-monthly Intervals. 

Figures 1 and 2 are examples of the calculated soil moisture and evapo- 
transpiratlon (Mintz and Seraflnl, 1981). 
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Soil nroisture (n*n) in mid-January (top); and mid-July (bottom). 
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SEASONAL VARIATION OF BLOCKING 
J. Shukla and K. C. No 


We have calculated the seasonal variation of frequency of blocking and Its 
geographical location by examining the grid point values of dally 500 mb geopo- 
tential height over the Northern Hemisphere for 15 consecutive years (1963-1977). 
Blocking events are Identified objectively by requiring that a large positive 
anomaly of a specified magnitude persist for 7 days or ir. 'e. The magnitude of 
the threshold anomaly Is assumed to be 200 gpm for winter. 100 gpm for summer, 
and 150 gpm for fall and spring. 

It Is found that the geographical locations of the maximum frequency, char- 
acterized by three distinctly different maxima, remains nearly the same In all 
the four seasons. These maxima coincide with the maxima of the low frequency 
and total variance. If the persistence criteria Is changed to 1-3 days, the 
geographical distribution of frequency for winter Is very similar to the band- 
pass variance, showing maximum values In the areas of storm tracks. Seasonal 
Invariance of the locations of maximum blocking events suggests the Important 
role of seasonally Invariant forcing due to mountains. 

Our principal conclusions are as follows: 

1) Large scale persistent anomalies of positive sign are more frequent 
than those of negative sign. As also noted by Charney et al. (1981). amplified 
blocking ridges persist longer than either deeper troughs or weaker ridges. 

2) There are three distinct centers of maximum blocking activity. They 
occur In the Pacific to the west of the Rockies. In the Atlantic to the west of 
the Alps and Scandanavlan Mountain ranges, and over land to the west of Ural 
Mountains of USSR. These relatively preferred locations of persistent anomalies 
do not change with season . 

3) A persistence criteria of 1-3 days represents storm tracks, and the 
resulting features are distinctly different from the persistence criteria of 
seven days or more. 

4) Occurrence of anomalous seasonal mean anomalies accounts for less than 
25X of the persistent anomalies identified as blocking events. 


REFERENCES 

Charney. J. 6., J. Shukla. and K. C. Mo. 1981: Comparison of barotropic blocking 
theory with observation. J. Atmos. Sd.. 38. 762-779. 


181 


B. 


CLIMATE MODELING 


» wtCED;r<vi i-'Aat Nor ^xiwcD 


preceding page blank not filme6 


‘ 25579 



THE SEASONAL-CYCLE CLIMATE MODEL 
L. Marx and D. Randall 


The GLAS climate model has undergone several modifications in preparation 
for the 400 day "seasonal cycle" run. The seasonal cycle run will become the 
control run for the comparison with runs utilizing codes and parameterizations 
developed by outside investigators. The GLAS climate model currently exists in 
two parallel versions: one running on the Amdahl at GNSF and the other running 
on the CYBER 203 at Langley Research Center. These two versions are as nearly 
identical as machine capability and the requirement for high speed performance 
will allow. Developemental changes are made on the Amdahl /CMS version for ease 
of testing and rapidity of turnaround. The changes are subsequently incorporated 
into the CYBER 203 version using vectorization techniques where speed improvement 
can be realized. 

The changes leading to the seasonal cycle version of the model can be sum- 
marised as follows: 


1. Conversion of model from SSS to CMS 

2. Removal of split-grid code 

3. Use of Fourier filtering to maintain computational stability near the 
poles 

4. Change to new tape format 

5. Faster long wave radiation code 

6. Use of potential temperature as the basic thermodynamic variable 

7. Revised negative specific humidity adjustment 

8. Revised dry convective adjustment 

9. New planetary boundary layer parameterization 

10. New moisture availability parameterization for land points 


The conversion of the code from SSS to CMS was designed to ensure that the 
model will be running on a maintained system with minimum overhead. This con- 
version required some code rewriting and creation of a completely new CMS EXEC 
procedure for operational runs. Many features available in SSS (such as sense 
switch •responses) have been included and in some cases expanded under CMS. The 
CMS version of the model allows for direct interactive debugging facilities for 
testing new codes and parameterizations. 

The split-grid in the original climate model has been found to generate 
spurious noise near the boundaries of the grid mesh changes. The new model uses 
the standard B-Grid mesh. The impact of this change on the CYBER 203 version 
of the model was substantial, since that code was designed exclusively for the 
split-grid mesh. 

The routine for maintaining computational stability near the poles was also 
found to be Inadequate. This led to the use of a Fast Fourier Transform (FFT) 
scheme to remove shorter waves in a more rational manner. With the use of the 
FFT filter, the same time step as was used in the split-grid code can also be 
used in the no split-grid code, thus minimizing the speed impact of the removal 
of the split-grid code. For land points, moisture availability is now computed 
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with the aid of a new prognostic “lyslmeter temperature," as suggested by Mintz 
and Serafim (1981). Preliminary results Indicate that the simulated distribu- 
tion of precipitation Is greatly Improved by this change. 

The new tape format was Introduced for the following reasons: 


1. Replacement of TAU (absolute hour) by a time/date system 

Z, Remove need for second tape for diagnostics and accumulate all useful 
diagnostics for model evaluation on one tape. 

3. Remove duplicate fields and remove unneeded variables from the C array. 

4. Segregate C array variables by type so that they can be treated together 
when transferring data from one machine to another. 

5. Reorder arrays and group them together at each latitude to accommodate 
a paging environment. 

6. Save code and boundary conditions for each run on the tape so that all 
code and boundary condition modifications can be seen. 

7. Maintain detailed records for a limited number of points to study the 
Impact of experiments over certain key regions of the world. 


Separate changes to the longwave radiation code have been made on the 
Amdahl and CYBER 203 versions of the model, yielding a 40% speed Improvement on 
each machine. This yields a 20% speed Improvement for the model overall. The 
results from both versions agree with the original longwave radiation code to 
within .1 degrees per day. 

The basic variable of temperature was replaced by potential temperature so 
that the thermodynamic advectlon equation could be simplified. 

The treatment of negative specific humidity has been changed so mat the 
adjustment Is made only after all updates have made to the specific humidity 
field, thus allowing other sources to compensate for the negative humidity 
where possible. Negative specific humidities are removed by a routine which 
ensures conservation of total moisture, thus eliminating spurious sources of 
moisture In the model. 

To ensure better control of the static stability, dry convective adjust- 
ment Is performed at every hydrodynamics and physics call instead of only at 
every physics call. 

A new planetary boundary layer parameterization has been incorporated into 
the model. This parameterization incorporates the surface flux parameterlza- 
atlon of Oeardorff (1972) as modified by Randall (1976). Realistic stability 
dependence Is Included, and the surface roughness has distinct values for land, 
ocean, and sea ice points. 

The changes to the model have yielded a single, up-to-date, and easily 
modified code that will serve as the reference run for Investigations of new 
numerical schemes and physical parameterlzations. The 400 day seasonal cycle 
run will serve as a control run for both medium and long range climate forecasts 
and sensitivity studies. By keeping the code simple a better understanding of 
model response to various Influences can be achieved. This In turn will lead 
to a better understanding of the atmosphere. 
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QUASI-LAGRANGIAN MODELS OF NASCENT THERMALS 
S. Rambaldl and D. A. Randall 


1. INTRODUCTION 

Scorer and Ludlam (1953) suggested that atmospheric buoyant convection 
occurs as discrete warm bubbles or "thermals" grow near the earth's surface 
and break away to rise under the influence of the buoyancy force. The motions 
in and around an isolated thermal were studied in the laboratory by Woodward 
(1959), who found rising motion in the core, and sinking motion on the outside; 
the circulation resembled that of a vortex ring. 

In an "entity cloud model," cloudy thermal is tracked, in a Lagrangian 
fashion, as a discrete entity; the field of motion in and around the thermal 
is not explicitly simulated. Such models cannot simulate the mutual interac- 
tions of neighboring clouds, since these interactions occur partly through 
mutual advection. 

In recent years, much effort has been devoted to the development of "field- 
of-motion cloud models," in which the equations of motion are numerically inte- 
grated on an Eulerian grid. 

There are several drawbacks to field-of-motion cloud models. They are 
expensive, and analysis and interpretation of the enormous volume of numerical 
results produced is a difficult task. Moreover, as Turner (1969, p. 35) has 
noted, "No way has been found so far to deal numerically with the sharp bound- 
ary between the turbulent buoyant fluid and its surroundings, which is so 
characteristic of the experiments." A very fine Eulerian mesh is required 
to resolve the sharp gradients which arise near the thermal's boundary; in 
other regions a much coarser mesh would suffice. The boundary must be resolved 
somehow, because that is where vorticity is generated, and where entrainment 
and detraininent occur^ 

These considerations show the great potential of a "hybrid" cloud model 
which could combine the simplicity of the entity models with the generality and 
flexibility of the field-of-motion models. A key problem to be overcome in the 
development of a hybrid model is the formulation of a mathematical framework 
within which the cloud dynamics can be represented. Details of the mathematical 
analysis are given by Rambaldi and Randall (1981). 


2. THE INITIAL ACCELERATION 

We have shown analytically that the initial acceleration of a cylindrical 
thermal of buoyancy B, which starts from rest in an unbounded neutral envir- 
onment, is as depicted on the right side of Fig. 1. Buoyancy generates positive 
vorticity on the right side of each thermal, and negative vorticity on the left 
side. The generation is most vigorous where the thermal wall is vertical. The 
induced environmental circulation is outward from the cap, downward along the 
sides, and inward and upward over the base. The rising motion is strongest on 
the center of the base. 
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Outside the thermal, there is a purely upward acceleration on the vertical 
plane of synmetry, and a purely downward acceleration on the horizontal plane 
of symmetry. On the boundary, the acceleration is purely normal, and the verti- 
cal component is nowhere downward. Across the boundary, the normal acceleration 
is continuous, while the tangential acceleration is discontinuous. Inside the 
thermal, there is a uniform upward acceleration. 

Similarly, we have shown that the initial acceleration of a spherical ther- 
mal in an unbounded neutral environment is as depicted on the left side of Fig. 1. 
There is clearly a great similarity to the slab- symmetric result, but several 
interesting differences are seen. One is that, for a given buoyancy B, the 
interior acceleration is more rapid for an axisymmetric thermal. The physical 
interpretation is that in the axisymmetric case the environmental air which must 
be pushed aside to allow the thermal to pass can spread out in two dimensions, 
while in the slab- symmetric case it can only spread out in one. As a result, an 
axisymmetric thermal can push through its environment more easily than a slab- 
symmetric thermal. For a similar reason, the environmental motion decays more 
rapidly with distance from an axisymmetric thermal. This is reminiscent of the 
numerical result, obtained by Murray (1970) and others, that compensating envi- 
ronmental subsidence is more vigorous around a slab- symmetric cloud than around 
an axisymmetric cloud. Finally, the initial acceleration on the boundary of the 
axisymmetric thermal is not purely normal, but has a tangential component. 


3. DISCRETE MODELS 

’n general, we cannot solve our equations exactly, so we adopt an approx- 
imate method suitable for use with a computer. The method closely parallels 
that used by Fink and Soh (1978) to study the roll-up of the trailing vortex 
sheet shed by a wing. 

We introduce finite-difference approximations to (2.11) and (2.12-13), 
taking care to avoid the singularities. These finite difference equations 
allow us to determine the velocity field associated with the distribution of 
vorticity along the boundary. The velocity field is used to predict the future 
positions of a finite number of equidistant points along the boundary. A "redis- 
tribution" algorithm is used to ensure that the points remain equidistant as 
the themwl evolves. Finally, the future values of the vorticity are predicted 
for each point. Gy repetition of this procedure we march forward in time. 


4. NUMERICAL EXPERIMENTS 

To study the organization of thermal circulations, and to demonstrate our 
models, we have simulated the motions of slab- symmetric and axisymmetric thermals, 
which are initially at rest in unbounded, neutral environments. We have used 
40 segments, and a finite ( nondimens ional) time increment of 2xl0“2, with the 
second-order modified Euler time-differencing scheme (Carnahan et al. , 1969). 

Fig. 2 shows the evolution of the thermal boundary up to a nondi mens ional 
time t* » 3.2, for both slab- symmetric and axisymmetric thermals. The more 
rapid rise of the axisymmetric thermal is clearly evident. Both thermals evolve 
from initially circular cross-sections into flattened domes with concave bases, 
in qualitative agreement with the laboratory results of Richards (1963) and 
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Walters and Davidso.. (1962, 1963). Dimple formation Is to be expected In view 
of the distribution of Initial acceleration. Although both thermals have 
broadened, the broadening Is much less pronounced for the axisymmetric thermal, 
apparently because It has spread out Into two dimensions. 

The growth of the "dimple" on the base Is accompanied by a shrinking of 
the cap thickness, for both thermals. The cap of the axisymmetric thermal Is 
close to breaking at the last time shown. The thermal would then become tor- 
oidal. The evolution of Initially spherical bubbles Into toroids has been 
observed In the laboratory by Walters and Davidson (1963). A cap-breaking 
process Is also occurring (more slowly) In the slab- symmetric thermal. 

Advectlon removes vorticity from the source regions on the sides, carrying 
It towards the base and Into the dimple. Inside the dimple buoyancy acts to 
destroy vorticity as It Is advected In. As a result, a local maximum of vorti- 
city occurs near the trailing edge of each thermal. A second very strong maxi- 
mum develops within the cusp. We did not anticipate the double vorticity 
maximum; Instead, we expected to find a single maximum, so that the thermal 
structure would approximate a vortex ring (or pair of vortex lines). If our 
model were generalized to Include realistic dissipation, the very strong and 
sharp vorticity maximum In the cusp would probably be weakened more severely 
than the broader, less Intense maximum on the trailing edge. The two maxima 
might then be more comparable In magnitude. 

Slightly beyond the last times shown, both thermal walls develop small- 
scale features near their cusps. This leads the walls to fold across themselves, 
so that the simulations cannot be continued. 


5. SUMMARY AND CONCLUSIONS 

We have presented new models of slab- symmetric and axisymmetric thermals, 
based on the Idea that the sharp, moving boundary of a thermal can be considered 
as a self-advecting vortex sheet. The analytically-obtained initial accelera- 
tions show that both thermals tend to form dimpled bases, and that the axisym- 
metric thermal accelerates upward more rapidly. Numerical simulations with dis- 
crete models confirm these results, and reveal the formation of double vorticity 
maxima. 

Our results encourage us to consider further applications. We hope that 
the models can eventually be used to simulate cumulus clouds. 
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AXISYMMETRY SLAB-SYMMETRY 


Figure 1. The analytically determined initial accelerations outside, on, and 

inside the thermal boundaries. The lengths of the arrows are measured 
relative to the buoyancy AB. The soVid arrows at the boundaries show 
the accelerations just outside the thermals, while the dashed arrows 
show the accelerations of the boundaries themselves. 
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Figure 2. The time-evolution of the thermal boundary. 
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SECOND-ORDER CLOSURE STUDIES OF ENTRAINMENT INTO A STRATOCUMULUS 
LAYER WITH DISTRIBUTED RADIATIVE COOLING 

D. A« Randall and C.-H. Moeng 


Recently the role of radiative cooling in detennining the structure of 
stratocumulus-topped mixed layers has become controversial {Deardorff 1976, 

Kahn and Businger 1979, Lilly and Schubert 1980, Randall 1980, Shaller and 
Kraus 1981, and Stage and Businger 1981). Mixed-layer modelers have recognized 
the importance of radiative cooling and tried to clarify this aspect. However, 
because mixed-layer models are highly sensitive to the choice of entrainment 
assumption, the conclusions of these studies are questionable. To avoid expli- 
cit dependence on an entrainment assumption, a second-order closure turbulence 
model is used here to show the sensitivity of a cloud- topped turbulent layer to 
the depth of radiative cooling. 

The one-dimensional second-order closure model is based on the Level 4 
model of Mellor and Yamada (1974), with some modificiations described in Moeng 
and Arakawa (1980). In order to be able to control the depth of the radiative 
cooling, we have omitted an explicit radiation calculation, and Instead simply 
prescribed the distribution of the radiative cooling. A net upward radiative 
flux R = 70 watts/m2 given at the cloud top. The flux is assumed to decrease 
linearly to zero at the level 6zr below the cloud top. Three 60-hour simulations 
with 6zr=20 m, 100m, and 200 m have been performed, all with the same initial 
and boundary conditions. The results show that 6zr increases, both cloud top 
and cloud base descend toward the sea surface, and the depth of the layer of 
negative buoyancy flux increases. The surface evaporation rates at the ends of 
the runs are 1.33, 1.21, and 0.86 mm/day for 6zr=20, 100, and 200 m, respectively 
the evaporation rate decreases as 6zr increases. 

These results are in qualitative agreement with those obtained by Randall 
(1980) using a mixed-layer model. The sensitivity of cloud- topped mixed layer 
structure to the depth of the radiative cooling suggests that in order to simu- 
late a realistic stratocumulus-topped boundary layer by a mixed-layer model, an 
accurate parameterization of the distributed radiative cooling is required. 
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SOME ASPECTS OF SOUTHWEST MONSOON AS SEEN IN A GLAS GCM SIMULATION 


V. Satyan 


1. INTRODUCTION 

In the past decade« several studies have been made of GCM simulations to 
see If they simulate the southwest summer monsoon or aspects of It. A review 
by Gilchrist Q977) discusses these attempts and their relative successes. In 
particular^ t^ie GLAS GCM simulations have also been studied with monsoon In 
view. This paper describes an attempt to Investigate In some detail two Impor- 
tant aspects of the summer monsoon, viz., disturbances In the monsoon flow and 
periodicities In the monsoon fluctuations as revealed In a GLAS GCM simulation. 

Results of GLAS GCM climate simulation for 50 summer days from June 15, 
1979 through August 3, 1979 (0162 C11 Run) were used In the present study. 

Dally mean fields were computed from these data sets for the monsoon region. 

The monsoon region here Is defined as the region bounded by latitudes 30**S-42*’N 
and longitudes 45®E-120®E. 


2. MONTHLY MEAN FLOW PATTERN 

First, It Is necessary to examine the model climatology for July In the 
monsoon region. Examination of July mean horizontal wind field at 850 mb and 
200 mb levels shows that the model simulates well the observed July mean south- 
westerly wind In the lower troposphere over the Arabian Sea and peninsular 
India. The 200 mb mean flow shows the anticyclonic flow over Tibet and the 
easterly jet over 6®N-14®N. Next, It Is also seen from the July mean sea level 
pressure that the ridge over the Arabian Sea, the monsoon trough over northern 
India, the Mescarene High east of Madagaskar and the low over Pakistan are 
simulated by the model and compare well with the observed July climatology. 
Thus, the overall monthly mean picture of wind flow and sea level pressure 
pattern are satisfactorily reproduced by the model. 


3. DISTURBANCES IN THE MONSOON FLOW FIELD 

A convenient way of studying propagating or stationary disturbances In a 
flow field Is by constructing Hovmoller diagrams. These are diagrams showing 
longitude-time cross sections of the ^low at fixed latitudes. Examination of 
the variation of sea level pressure 80®E and 22®N during the 50 day period 
shows that a major low formation occurs about the 13th day (1st day Is June 15, 
1979). Furthermore, the plots of precipitation and the wind speed at 1000 mb 
level also show peaks around this time. We present the Hovmoller diagram at 
22®N In Figures la and lb. It Is seen that a low forms around 80®E on the 
13th day. Furthermore, this disturbance propagates westward with a phase speed 
of about 2.7 degrees longitude per day (slope of the trough-trough line). 
Observational statistics show that disturbances frequently originate In the 
Bay of Bengal region during July and they usually move In northwesterly direc- 
tion along the monsoon trough axis. It Is also Interesting to note the forma- 
tion of another disturbance around the 33rd day and this Is accompanied by more 
rainfall than the major low described above. Hovmoller diagrams at latitudes 
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1A*’N and 26**N (not shown here) are similar to Figure 1 and reveal the lati- 
tudinal extent of the disturbance. He conclude that the model has produced a 
genuine physical monsoon disturbance and not a numerical Illusion. 


4. PERIODICITIES OF A SOUTHWEST SUMMER MONSOON SYSTEM 

The question of periodicity In the fluctuations of a monsoon has received 
some attention In the past. The studies made so far have been confined to the 
analysis of observed data. The subject has been reviewed by M. Murakami (1977). 
The conclusion thot emerges from these studies Is that the monsoon system exhi- 
bits two kinds of quasi -periodic variations, one around the 5-day period and 
the other around the 15-day period. The question arises, therefore, whether a 
ACM simulation also shows these periodicities In the monsoon region and whether 
they resemble the observed periodicities. We Investigate this question In the 
foliowlnq, using the results of a GLAS-GCM simulation. 

The time series analyzed consists of 50 days of dally mean fields computed 
from the 0162 summer simulation referred to earlier. The fields analyzed are 
the pressure In the monsoon trough area, pressure In the Mescarane High, tropi- 
cal easterly Jet and the monsoon rainfall. These are some of the major large- 
scale observed features which define, along with others, the monsoon system. 

In the present study, these fields are averaged over almost the same region as 
In the observational study of Krishnamurty et a1. (1976) so as to afford an 
easy comparison with their findings. Before carrying out spectral analysis, we 
detrend the time series (departure from the 50-day mean) by subtracting out the 
linear and parabolic trends so that the red noise and the long-term variations 
(like the seasonal) are removed. The spectrum Is studied by plotting power x 
frequency against logarithm of the frequency. 

Figures 2a and 2b show the day-to-<idy variation of the pressure In the mon- 
soon trough and Its spectral analysis., respectively. The pressure here is the 
dally mean sea level pressure averaged over the area bounded by latitudes 18®N- 
26®N and longitudes 60®E-100®E. The spectral diagram shows several peaks; the 
prominent ones havl.ig periods around 16.6-day, 10-day, 6-7 day and 5-day. Of 
these, the 16.6-day and 5-day periods may be compared with the 15-day and 5-day 
periods found in the analysis of observed data (M, Murakami, 1977). Ue also 
note that the time series In Figure 2a does show oscillations around the 16th 
day. 


We present In Figures 3a and 3b the results of our analysis for the 
Mescarane High, This Is the high pressure area to the east of Madagascar. 
Figure 3a shows the dally march of the pressure (averaged over the area bounded 
by latitudes 22®S-30®S and longitudes 45®E-70®E), The spectrum Is shown In 
Figure 3b. He find three distinct periodicities at 16.6-day, 7.1-day and 5,1- 
day. This result agrees well with the observational studies (Krishnamurty et 
a1 . . 1976) except for the 7.1-day oscillation. 

The dally march of precipitation and Its spectrum are shown In Figures 4a 
and 4b. Figure 4a shows the dally mean monsoon rainfall averaged over the area 
covered by 10®N-26®N and 70®E-90®E (cjintral India). Considerable fluctuations 
seen In the precipitation vs time plot Indicates that Its spectrum is ccwnposed 
of several periodicities. In fact, we see In Figure 4b that there are six dis- 
tinct periods. Besides the 16.6-day and 5-day oscillations as In the oscll la- 
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tions of other elements of the monsoon system one sees here prominent ones at 
7,1-day, 10-day and 2.6-days. It is Interesting to note that the 2.6-day 
period agrees very well with the 2-day period oscillation found by Krishnamurty 
et al. (1976) in their analysis of observed data. 

We also analyzed the tropical easterly jet for periodicities. The para- 
meter studied is the daily mean wind speed at 200 mb level averaged over area 
covered by latitudes 2®N-10®N and longitudes 55®E-75®E, We find that the 
16.6-day period is very prominent in the spectrum. Again we see periodicities 
of 7.1, 5.5 days and some additional ones. The high frequency end of the 
spectrum shows several peaks, which may be the consequence of the local baro- 
tropic instability of the jet (Krishnamurty et al,, 1976). W6 also find that 
the wind speed fluctuates considerably with time. 

Thus, the GLAS GCM simulates the periodic variations of the summer monsoon 
system fairly well. The 5-day and 15-day periods found in the analysis of 
observed data are reproduced closely by the model. In addition, the model 
results show a 7.1-day oscillation in all the elements of the monsoon considered 
above. It is interesting to note that the existence of a mode around 7-day 
period has been discussed by Keshavamurty (1972, 1973) in his analysis of 
observed data* 
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AN ANALYSIS OF THE CLOUD FIELDS SIMULATED 8Y THE GLAS CLIMATE MODEL 

C.-H. Moeng and D. A. Randall 


In the present GLAS Climate model, It Is assumed that clouds of any type 
fill an entire grid area, which Is typically about 400 km x 400 km In horizon- 
tal and 1-2 km in vertical. This cloud cover Is then used to calculate the 
radiative heating/cool Inq of the atiiosphere-earth system. Subgrid-scale cumulus 
convection Is observed to have a very small horizontal extent and, because no 
parameterization of subgrid-scale fractional cloudiness Is currently available. 

It seems reasonable to Ignore the existence of cumulus clouds for the radiation 
balance. 

The purpose of this study Is to analyze the cloud field produced by the 
model, and try to answer the following questions: (a) does the model produce 

reasonable cloudiness? (b) Is the amount of model -generated cumulus cloudiness 
comparable to supersaturation cloudiness? (c) Is the effect of cumulus clouds 
on radiation balance large? 

To compare the model's cloudiness with observations, we have chosen February 
(D150) and July (0162) simulations. In order to make a clean comparison with 
satellite observations (or ground observations), we define cloudiness such that 
It represents the horizontal distribution of cloud-cover as one looks straight 
down from the top of the atmosphere (or straight up from the ground); the cloud- 
iness Is 100% If any cloud exists in the vertical grid column in any model 
layer, and zero If there Is no cloud in any layer. These analyses have been 
made for cumulus clouds, supersaturation clouds, and both types clouds. Individ- 
ually. The method Is as follows: (a) unpack cloud Information at each grid 

point from log-9 model tapes; (b) compute the horizontal distribution of cloud- 
iness from cumulus clouds, supersaturation clouds, and total clouds; (c) take 
the zonal mean for each type of cloud; (d) repeat the above steps for one-month 
of data, and take a monthly average. 

The February and July zonal -mean cloudinesses generated by the model are 
shown In Fig. 1. Comparing with those of observations (Fig. 2), we find good 
agreement. There Is about a 10% overestimation of the global cloudiness. The 
overall shape of the zonal mean cloudiness Is very well -simulated. 

The zonal mean cloudiness due to cumulus clouds between ^ 60** latitudes Is 
only 1/3 of that due to supersaturation clouds. The total cloudiness of the 
model Is mainly due to supersaturation clouds. In other words. If we completely 
Ignore the cumulus clouds, the model -produced zonal-mean cloudiness Is still 
comparable to that of observations. 

Although the cloudiness due to cumulus clouds is small, this does not mean 
that the cumulus clouds occur less frequently In any particular model layer. 

We have analyzed the July simulation to see the distribution of the monthly 
cloud frequency In each model layer. The monthly cloud frequency Is defined as 
the fraction of the time for which cloud occurs within the one month period. The 
cloud Information Is retrieved from log-9 model tapes for every 5-hour output, 
then the time average of the appearance of cloud at each grid point Is taken. 

The monthly cloud frequency Is obtained for cumulus clouds and supersaturation 
clouds, separately. For cumulus clouds In each layer, we sum over different 
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types of convective clouds. For example, to obtain the cloud frequency of 
cumulus clouds in layer 7 we sum that of low level convection and deep con- 
vection in layer 7. Examination of these results (not shown) reveals three 
features. First, in the mid-troposphere (layers 6 and 7) cumulus clouds occur 
more frequently in the low and middle latitudes and supersaturation clouds are 
confined to the polar regions. Second, most supersaturation clouds appear in 
the lowest layer. The supersaturation cloud frequency in the lowest layer is 
unrealistically large. The sudden decrease of supersaturation cloud frequency 
from layer 9 to layer 8 suggests that too much moisture is deposited inside the 
lowest layer; the model is unable to transport low-level moisture upward. 

Third, the model produces unrealistic supersaturation clouds in the stratosphere 
(layers 1, 2, and 3). 

Many GLAS climate simulation experiments have shown that the model is sen- 
sitive to the parameterizations of radiation and to the treatment of the inter- 
actions between cloud and radiation (Wu et al., 1978; Herman et al., 1980; 

Shukla and Sud, 1980). This report shows that the scattered, small-scale cumulus 
clouds occur more frequently than the horizontally wide-spread supersaturation 
clouds in the low latitudes and mid-troposphere. Hence, a systematic bias 
exists in the model because the model treats the horizontal extent of cumulus 
clouds and supersaturation clouds the same. We suggest that all the subgrid- 
scale cumulus clouds should be ignored for the radiation budget. 
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Fig. 1. Latitudinal distribution of zonal ly averaged cloudiness for 
(a) February and (b) July, from model output. (Dashed for 
cumulus clouds, dotted for supersaturation clouds, and solid 
for total clouds). 



Fig. 2. Latitudinal distribution of zonal ly averaged cloudiness for 
(a) January and (b) July. The Landsberg and Vowinckel and 
Van Loon results are from ground observation. Other curves 
from satellite observations as marked. (After Clapp, 1964). 
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CLIMATOLOGY OF BLOCKING IN THE GLAS CLIMATE MODEL 
J. Shukla, K. C. Mo, and M. Eaton 


1. INTRODUCTION 

Durinq the past several years, the GLAS climate model has been Integrated 
with a variety of Initial and boundary conditions. Results of winter and summer 
simulation by the GLAS climate model have been discussed by Halem et a1 . , (1979) 
and Straus and Shukla, (1981). The model has shown remarkable succesTTn simu- 
lating the storm tracks. The purpose of this note Is to examine the model's 
ability to simulate the occurrence of persistent anomalies to be referred to as 
blocking. 


2. DATA 

We have used 17 winter and 7 summer simulations with the GLAS climate model 
(Halem et al. . 1979; Shukla and Bangaru, 1980). The time series at each grid 
point for each run was fitted to a parabola to define the seasonal cycle, ^s, 

«^ 5 (t) = a + bt + ct^ 

The mean values of coefficients 7 , F, and c were determined by averaging the 
coefficients over all runs In a set and seasonal cycle was defined as: 

^j(t) * 7 + Ft + 7t^ 

at each grid point. Anomaly at each grid point was defined as: 

f (t) = m -?s(t) 

if, at any grid point, an anomaly of +100 gpm or more persists for 7 days 
or more, we count it as a blocking event. We have calculated the frequency of 
blocking events at each grid point for winter and summer simulations separately 
and presented the results In the next section. 


3. RESULTS 

Figures la, lb, Ic, and Id show the standard deviations of dally values of 
qeopotentlal height for winter and sMmmer, respectively. Model standard devia- 
tions look more like observed band-pass variances rather than observed total 
variances. This suggests that the model does not simulate the low frequency 
components realistically (Straus and Shukla, 1981). 

Figures 2a and 2b show the map of number of events for which anomaly of 
+100 gpm or more persisted for 7 days or more. Figure 2a shows only the Northern 
Hemisphere, and Figure 2b also shows the Southern Hemisphere. In the Northern 
Hemisphere, the centers of maxima are located near 50**N, IBO^U (Pacific), 75"N, 
30“W (Greenland), 60®N, 40«*E (USSR), and 70®N, 135«W. An analysis of 15 years 
of observations (Shukla and Mo, 1981) showed only three centers of maxima; more- 
over, the locations of the maxima are not correctly simulated. The Pacific and 
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the Russian maxima are shi<^ted southwestwards, and the Atlantic maxima Is shifted 
northwestwards with respect tu the observations. It Is not appropriate to com- 
pare the magnitude of the frcr^uency because rvdel nins have very similar Initial 
conditions. Shukla and Bangaru (1979) found that a large SST anomaly In the 
Pacific produced a persistent anomaly In the atmospheric circulation over North 
America. This situation Is being further analysed by Chen and Shukla and It Is 
found that SST anomaly run produced persistent low wave numbers. 

Figure 2b shows that during northern winter the Southern Hemisphere has 
relatively lower frequency of blocking. The maxima In blocking occur between 
70**-80®S. We do not have corresponding observations to compare the Southern 
Hemisphere results. 

Figure 2c shows the frequency of blocking for summer simulations. There 
Is almost no blocking event In the Northern Hemisphere (the maxima over Green- 
land Is not a realistic one), and surprisingly, there Is no blocking activity 
In the Southern Hemisphere either. 

We have also examined the number of events for winter simulations for ano- 
malies of -100 gpm or less persisting for 7 days or more. In agreement with 
the observations, maximum frequency maps for simulated positive and negative 
anomalies occur In the same region. 

The model's day-to-day fluctuations are dominated by synoptic scale distur- 
bances which give maximum variance In storm track regions, whereas. In the 
observations, the maxima of dally standard deviations are related with large 
persistent positive and negative anomalies associated with blocking. 

We have calculated the characteristic time Tq for winter model runs and 
15 observed winters, respectively. 


T 

To = 2 / U- T J R(^) dT 

where R(t) is lagged autocorrelation at lag t. and T * 30 days. We find that 
the model underestimates the value of Tq. Indicating the model's inability to 
simulate persistent patterns and low frequency components. 


4. CONCLUSIONS 

Although the GLAS climate model simulates the storm tracks (band-pass varl 
ance) quite realist'ically, the model Is not successful in simulating the occur- 
rence of persistent anomalies. Examination of model simulated synoptic maps 
shows that the blocking type of configurations appCi^r frequently, but thqy do 
not persist. 
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Fig. Id- Global standard deviations for summer simulations. 
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Fig. 2a. Number of events at each grid point N.H. for which positive 
anomaly of 100 gpm or more persists for 7 days or more for 
winter simulations. 
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Fig. 2b. Number of events at each grid point (global) for which 
positive anomaly of 100 gpm or more persisted for 7 days 
or more for winter. 
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STRUCTURE AND .DYNAMICS OF MONSOON DEPRESSIONS; THE MONEX 
' DEPRESSION (JULY 1979)1 

J. Shukla 


1. INTRODUCTION 

It is qenerally believed that althouqh the Initiation of the onset of the 
Asiatic monsoon Is related to the planetary scale heat sources and sinks (due 
to asymmetric continentallty and different thermal properties of land and ocean), 
the maintenance and the short term variability of the monsoon Is primarily 
determined by the release of the latent heat of condensation, a major portion 
of which occurs In association with synoptic and large scale disturbances. The 
monsoon disturbances therefore play a very Important role in the dynamics of 
the planetary scale monsoon circulation. Preliminary studies based on the FGGE/ 
MONEX data tend to support the conjecture that even If the large scale thermal 
and dynamical circulations are conducive for onset, the onset process requires 
a finite amplitude perturbation for northward propagation and establishment of 
the monsoon over India. 

There are several classes of monsoon disturbances (depending upon their 
space-time scales), the most Important of these being the east-west oriented 
monsoon troughs which move north and south, and the monsoon depressions which 
'appear' over the Bay of Bengal and move northwest over India. In this paper 
we shall confine our discussion to the mechanisms responsible for the formation 
of monsoon depressions. 

There have been many suggestions about the mechanisms for the formation of 
monsoon depressions. During the first quarter of the century, there were attempts 
to explain the dynamics of monsoon depressions based on extra-tropical frontal 
dynamics models which were subsequently abandoned when upper air soundings failed 
to show sufficient ba rod Ini city. Subsequent attempts to explain the dynamics 
of monsoon depressions were Influenced by the quasi -geostrophic development 
theory. In recent years, empirical evidence has been presented to show the 
relationship between the formation of monsoon depressions, the presence of warm 
pools at 200 mb (Moula, 1968), and decreased vertical shear (Raman et al. , 1978). 

The possibility of monsoon depression formation being forced by dynamical 
Instabilities of the mean monsoon flow was examined by Shukla (1977, 1978). 
Barotropic instability analysis was carried out for each level separately and 
It was found that both the lower and the upper levels were barotropically 
unstable. However, since the upper levels, due to the presence of the westerly 
and the easterly jet streams, had the largest available kinetic energy, the 
growth rates and the amplitudes were largest for the upper level jet. The most 
unstable mode for the joint barotropic-barocllnic Instability was also found 
to be dominated by the upper level Instability whose phase speed was larger 


1 Published In Condensed Papers and Meeting Report: International Conference 
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217 


than that of a typical rwnsoon depression. A joint CISK-barotropic-baroclinic 
instability analysis of the monsoon flow suggested that the latent heat of 
condensation is the most crucial source of energy for the growth and maintenance 
of the monsoon depressions. It is interesting to note that during the 1880's 
Eliott had suggested that the latent heat of condensation may be the primary 
energy source for initiation and maintenance of monsoon depressions; however, 
a quantitative treatment of interaction between latent heating and dynamics 
was not possible until recently when Charney and Eliassen (1964) and Ooyama 
(1964) developed the CISK theory. Based on CISK*barotropic-baroclinic insta- 
bility studies it was concluded that the barotropically unstable lower layers 
are the primary triggers for the onset of the instability which is further 
amplified by the latent heat of condensation. In addition, it was conjectured 
that the role of the terrain is to force a large scale flow which is barotropi- 
cally unstable at lower levels. 

Recently, Mishra and Salvekar (1980) have shown that baroclinic instability 
alone can account for the growth of the monsoon depressions. This erroneous 
conclusion is the consequence of the choice of a highly unrealistic meridional 
temperature gradient at the ground and arbitrarily chosen vertical wind shear 
at the lower levels. In the classical baroclinic instability problem (Eady, 
1949), the primary forcing comes from the meridional temperature gradient at 
the surface, and the choice of an unrealistic vertical shear at the lower levels 
is equivalent to specifying an unrealistic meridional temp'* ’'ure gradient at 
the ground. In the real atmosphere the observed vertical s .. near the surface 
is largely determined by the frictional forces and the *:nund' j layer dynamics. 
The internal let instability (Charney and Stern, 1962) of the upper level flow 
calculated by Mishra and Salvekar is in agreement with the earlier calculations 
(Shukla, 1977; Goswami et al. , (1980)); however, it does not provide the expla- 
nation for the monsoon depressions. Satyan et al. (1981) have calculated the 
barotropic-baroclinic instability of the monsoon flow which includes vertical 
shear of the meridional wind. As it was implicit in the work of Eady (1949), 
vertical shear of the meridional wind is always unstable because there Is no 
stabilizing beta effect. The results of Satyan et al. therefore appear to be 
a mere artifact of the choice of large vertical shear in the meridional wind 
for which there is no observational evidence. 

In most of the instability studies of the monsoon flow there has been too 
much emphasis on the co.icept of a 'preferred scale.' It is fairly clear by now 
that the latent heat of condensation is the most important energy source for 
the monsoon depressions. It is important to note that the concept of maximum 
growth rate as a criterion for the preferred scale, which Is quite reasonable 
for shear instabilities, is not necessarily suitable for the condensation-driven 
instabilities. As was pointed out by Shukla (1978), a large growth rate does 
not necessarily imply that the wave will attain the maximum amplitude and domi- 
nate over the other waves. The fastest growing wave will equilibrate faster 
and therefore the most preferred wave will be determined by its ability to max- 
imize the utilization of the available moisture. The potential for dominance 
is therefore calculated as the ratio of the imaginary and the real parts of the 
complex phase speed. These considerations indeed gave the most preferred scale 
to be about 3000 km, which is contrary to the statements by Keshavamurty et al. 
(1978) and Mishra and Salvekar (1980) that the earlier studies failed to get a 
preferred scale* 
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2. WHERE no THE MONSOON DEPRESSIONS FORM? 


It has tM*en nenerally assumed that the monsoon depressions originate over 
the Bay of Bengal and, therefore, most of the Instability studies have examined 
the flow over this region. A recent study by Saha et a1. (1981) has shown that 
during the 10-year period, 1969-1978, more than 80% o^ the monsoon depressions 
appearing over thg Bay of Bengal were associated with predecessor disturbances 
coming from the east. This suggests that the cyclogenesis over the Bay mainly 
Involves the amplification of a pre-existing disturbance and one need not have 
to Investigate the conditions for the growth of Infinitesimal disturbances. 
Since the sea surface over the Bay Is the warmest at that latitude and the low 
level flow Is convergent due to the presence of the monsoon trough, the CISK 
mechanism appears to be the most appropriate to explain the growth of the mon- 
soon depressions. It remains to be clarified as to why for some weak distur- 
bances CISK mechanism leads to rapid growth and for others It does not. 

It Is Interesting to note that the westward propagating disturbances In 
the tropics seem to be ubigultous as they can be traced back to the Pacific, 
the Atlantic and the African land masses. The formation of the monsoon depres- 
sions, the Intense easterly waves, and the cyclonic storms, etc., depends upon 
the existence of suitable large scale environment where these pre-existing weak 
disturbances can grow. 


3. FORMATION OF THE MONEX DEPRESSION (JULY 1979) 

The cyclonic circulation over the Bay of Bengal was first observed on 
July S at about 500 mb (there was no data above 500 mb) and It later descended 
to the lower tropospheric levels. This behavior was also noticed In the study 
of the past cases by Saha et al. (1981) and Raman et al. (1978), Me know that 
the upper level easterly jet always satisfies the necessary conditions for baro- 
troplc Instability and Internal jet Instability. Is It possible that the 
growth of the amplitude In the lower troposphere Is caused, at least In part, 
by the downward radiation of wave energy? This possibility Is being examined 
by Held and Desmukh (personal communication). However, since the upper level 
waves propagate westward with a large phase speed. It needs to be examined 
whether their stay over the Bay Is long enough to build sufficient low level 
amplitude. 

It has been shown by Saha and Shukla (1980) that there was some evidence 
of a westward propagating disturbance which may have amplified Into a monsoon 
depression over the Bay of Bengal. It has also been shown by Nitta and Murakami 
(1980) that the lower level flow over the Bay was barotropically unstable dur- 
ing the growth of this depression. 

It was noticed during the field phase that the extent and the organization 
of the clouds associated with the depression was not well defined and In spite 
of a well defined dynamical circulation, the cloudiness was poorly organized. 

It was also observed that the amplitude of the disturbance was prominent only 
in the lower and middle troposphere; the amplitude at 900 mb and below was rather 
weak. This leads to a speculation that rapid growth of the monsoon depression 
over the Bay of Bengal depends upon Its ability to utilize the latent heat of 
condensation by drawing on the boundary layer moisture convergence. If the 
Initial disturbance does not have sufficient amplitude In the lowest layers. 
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where the mixing ratio Is the highest, moisture convergence and associated 
latent heating may not be adeguate to Intensify the disturbance. We do not 
understand why some of the disturbances develop large amplitudes In the lowest 
layers and others do not. Ue also do not understand what determines the down- 
ward propagation of amplitude intensity. We hope that further theoretical 
studies may clarify the mechanism responsible for some of these observed fea- 
tures. It should be pointed out, however, that the phenomena of downward build 
up of the disturbance amplitude Is also noticed for the growth stage of the 
tropical cyclones and, therefore. In some respects, the growth of a monsoon 
depression Is not very different from the growth of a tropical cyclone. The 
monsoon depressions, however, move over land only after a few day's stay over the 
Bay and therefore do not Intensify Into tropical cyclones. 


4. MOVEMENT OF MONSOON DEPRESSIONS 

It Is generally observed that the phase velocity of the monsoon depressions 
over the Bay Is smaller than that over the adjoining land. P re -depress Ion dis- 
turbances coming from further east have different phase velocities depending 
upon the level at which their amplitude Is most pronounced. The low level flow 
over the Bay Is generally westerly and If the depression has large amplitude at 
the lower levels, the steering effect will produce eastward movement. The beta 
effect and the lower boundary slope effect will produce westward movement. The 
role of the quasi -geostrophic dynamics, by the combined effects of differential 
vorticity advectlon and thickness advectlon In determining the phase velocity, 
will depend upon the structure of the large scale and the embedded disturbance. 
Preliminary computations by Sanders (1981) Indicated that the quasi -geostrophic 
dynamics was not very dominant In determining either the growth or the movement 
of the depression during the developing stage. 

It Is suggested that the rapid westward propagation Is related to the ver- 
tical structure of the monsoon depressions. If the depression attains suffi- 
cient amplitude at the upper levels, the steering effect of the upper level 
easterlies contributes to the westward movement. The slow phase velocity of the 
depressions over the Bay permits a faster growth and development of a vertically 
coupled disturbance which can then be Influenced by the upper level easterlies. 

If the disturbance does not have sufficient vertical coupling, which may 
be either due to weak dynamical forcing and/or due to lack of organized moist 
convection, the vorticity maxima at different levels move with different phase 
velocity In different directions and the disturbance develops strong vertical 
tilt which Is followed by the weakening and decay of the disturbance. This 
was the case for the MONEX depression on 7-8 July, 1979, which developed strong 
vertical tilt and the magnitude of the tilt was accounted for by the vorticity 
advectlon at different levels (see Sanders, 1981). 

The foregoing survey of the mechanisms of formation and movement of monsoon 
depressions can be summarized as follows: 

(1) Monsoon depressions over the Bay of Bengal are caused either by ampli- 
fication of westward propagating weak disturbances or by downward propagation 
of the Internal Jet Instability of the easterly Jet. 



(2) The barotropic Instability of the low level flow over the Bay of Bengal 
Is conducive to the growth of weak disturbances. The presence of the surrounding 
terrain contributes to the establlshtnent of such a large scale flow (nx>nsoon 
trough over the head bay) which Is barotroplcally unstable. 

(3) CISK Is the primary driving mechanism for the rapid growth of a pre- 
existing weak perturbation. However, If the amplitude of the disturbance at 

the lower levels Is not sufficient to draw on the low level moisture convergence, 
the disturbance does not grow to a deep depression. 

(4) Once the disturbance has attained adequate amplitude In the lower levels 
(either by downward propogatlon of wave energy or by CISK), the low level mois- 
ture convergence and latent heat of condensation is utilized more efficiently 

for the development of a vertically coupled deep disturbance which along with 
upper level easterlies contribute to the westward propagation of the disturbance. 
Absence of strong vertical coupling leads to vertical tilt and decay of the 
disturbance. 

(5) Our ability to forecast the formation of monsoon depressions over 

the Bay of Bengal should be greatly Improved with the establishment of upper air 
stations over the Burmese coast and further east. 
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NUMERICAL PREDICTION OF THE MONSOON DEPRESSION OF 5-7 JULY. 1979^ 
J. Shukla. R. Atlas, and W. E. Baker 


1. INTRODUCTION 

One of the primary objectives of the FGGE/MONEX effort is to determine the 
three dimensional structure of the atmosphere more accurately, thereby improving 
the quality of deterministic short range prediction. The detailed instantaneous 
structure of the monsoon circulation, which is perhaps the largest asymmetric 
perturbation of the general circulation of the atmosphere, had never been defined 
earlier because of the paucity of data over the Indian Ocean and the adjoining 
countries. Numerical predictions of the monsoon fluctuations, therefore, pro- 
vide a good test to determine the impact of the FGGE/MONEX observing system on 
short range prediction. Although the large scale monsoon circulation is forced 
by planetary scale heat sources due to asymmetric continental ity and differen- 
tial heat capacity of land and ocean, the short term fluctuations are mainly 
caused by the synoptic scale disturbances referred to as the monsoon depressions. 
One such depression appeared over the Bay of Bengal during the first week of 
July 1979. This, fortunately, happened at a time when a large group of MONEX 
research scientists had gathered at Calcutta, India to observe and study the 
summer monsoon circulation. There were two MONEX research aircrafts (Electra 
from NCAR and P3 from NOAA) which could probe the formative and the growth 
stages of the monsoon depression. A well-defined monsoon depression was observed 
on July 7, 1979 which then moved over India and dissipated. 

We have chosen this particular synoptic situation for two assimilation 
and forecast experiments. Starting from the initial conditions (provided by the 
NMC global analysis) of 0000 GMT 1 July 1979, one experiment utilized only the 
conventional surface and upper air data, and the other utilized the same conven- 
tional data plus all the available FGGE/MONEX data. We refer to these two 
assimilations as the control assimilation and the FGGE/MONEX assimilation. From 
the two initial conditions valid for 1200 GMT 7 July 1979, arrived at by assi- 
milating the two different data sets, we have made numerical predictions with 
the GLAS general circulation model. In the following sections, we describe the 
model, the analysis and assimilation procedure, the differences in the analyses 
due to different data inputs, and the differences in the numerical predictions. 


2. ANALYSIS AND ASSIMILATION OF FGGE/MONEX DATA 

In this section we describe the objective analysis and assimilation proce- 
dure used with the FGGE/MONEX data. 


1 Published in Condensed Papers and Meeting Report: International Conference 
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2a. The Model 


The model is a global, primitive-equation model, discretized in finite- 
difference form. There are nine vertical layers equal in sigma with a horizontal 
resolution of 2.5" in latitude and 3“ in longitude. The horizontal differences 
are second-order accurate on a staggered grid (Arakawa, 1966; 1972). The model 
differs from that described in Somerville et al . (1974) in its horizontal 
resolution and by the introduction of a split -grid, which at high latitudes 
modifies the longitudinal mesh size. 

Time differencing is performed following the space-centered Euler-backward 
(Matsuno) scheme (Arakawa and Lamb, 1977). Due to a number of features of the 
difference scheme, such as the averaging required by the staggering of variables 
and the additional smoothing at high latitudes, the scheme is strongly dissipa- 
tive. As a result, spurious high-frequency oscillations generated at the begin- 
ning of a forecast decay rapidly. 

2h. The Data 


For this study we utilized FG6E data collected from 1-7 July 1979 and the 
special MONEX data (dropwindsondes, aircraft, and enhanced TIROS-N retrievals). 
The data nominally available to the objective analysis program over the entire 
globe is listed by type in Table 1. Table 2 summarizes the MONEX data utilized 
over the Bay of Bengal. 


Table 1. Typical number of "Soundings" for each 24 h period 
during the assimilation cycle (July 1-7, 1979). 


Rawinsondes: 

3000 

Satellite winds: 

5500 

Satellite temp.: 

8500 

Aircraft: 

3500 

Constant level balloons: 

250 

Enhanced satellite temp. 

(Bay of Bengal): 175 


The rms wind vector fit (final analysis vs. observations) for dropwindsondes 
and Wisconsin Indian Ocean clouds tracked winds was compared with that for 
rawinsondes. The rms fit for the dropwindsondes was quite similar to that for 
the rawinsondes. The Wisconsin winds exhibited a poorer fit due to less weight 
given in the analysis of that data and the higher data density. 

An enhancement of the TIROS-N sounding data for the Bay of Bengal was per- 
formed in order to supplement the operational sounding data set with higher 
resolution soundings in meteorologically active areas, and with new soundings 
where data voids or soundings of questionable quality existed. Man-computer 
Interactive Data Access System (McIDAS) terminals, developed by the Space 
Science Engineering Center of the University of Wisconsin, were utilized for 
the display and enhancement of the TIROS-N soundings. 
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The alqorithms for retrieving temperature profiles from the TIROS-N radi- 
ances were essentially the same as those used by NESS for the objective genera- 
tion of operational temperature profiles. Three types of temperature retrievals 
are possible; (1) clear column, (2) partly cloudy retrievals utilizing Infrared 
observations from the High Resolution Infrared Sounder (HIRS) Instrument aboard 
TIROS-N, and (3) cloudy retrievals utilizing only microwave observations from 
the Microwave Scanning Unit (MSU) on TIROS-N. The operational temperature 
profiles have a horizontal resolution of 250 km, whereas enhanced temperature 
profiles can be retrieved at the resolution of the measurements. (30 km for 
HIRS and 150 km for MSU). 

Subjective comparisons of enhanced and operational soundings and rawinsonde 
reports ^r the FGGE "Special Effort" (Atlas, 1980) have shown that the enhanced 
and operational retrievals tend to be similar In cloud free areas. However, 
large Improvements In thickness and mandatory level temperature, and Intensifi- 
cation of atmospheric thermal gradients occasionally result from the enhancement 
of cloudy and partially cloudy areas. 

2c. Objective Analysis Procedure 


In the GLAS objective analysis scheme (Baker et al. , 1981) eastward and 
northward wind components, geopotential height and relative humidity are 
analyzed on mandatory pressure surfaces. The 6 h model forecast provides a 
first guess for these fields at 300 mb and sea level, where sea level pressure 
and sea level temperature are also analyzed. The first guess for the other 
levels Is obtained from the model first guess modified by a vertical interpola- 
tion between the two closest completed analyses. Vertical consistency Is main- 
tained through static stability constraints. The analysis at each level Is 
performed with a successive correction method (Cressman, 1959) modified to 
account for differences In the data density and the statistical estimates of 
the error structure of the observations. The average distance d between data 
points Is found In a circle with a radius of 800 km centered at each grid 
point. Three scans are performed with a radius of Influence R^ = c^d, where 
the coefficients c^ (1.6, 1.4, 1.2) were chosen to minimize the analysis error 
(Stephens and Stitt, 1970), However, the radius of Influence is not allowed 
to become smaller than 300 km. During this process, all data are checked for 
horizontal consistency. The completed analyses are smoothed and then Interpo- 
lated to the model sigma levels. 

The assimilation procedure provides for the intermittent analysis of batches 
of data grouped In a 3 h window about each synoptic time. In these experi- 
ments, the wind and height fields were analyzed independently with no explicit 
coupling or balancing. 

2d. Differences between the control and the FGGE/MONEX assimilations 


Fig, 1 shows the zonal ly averaged root mean square (rms) vector wind error 
between the control assimilation and the FGGE/MONEX assimilation for 1200 GMT 7 
July 1979. It can be seen that the maximum differences occur In the southern 
hemisphere, especially near the poles. Since there Is virtually no data In the 
control assimilation over that part of the globe, the control analysis on July 
7 could not be very different from the model prediction except for the Influences 
that might have propagated from the data-rich regions. These differences are 
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therefore, an indication of the systematic model forecast errors. The differ- 
ence is the smallest over the northern hemispheric mid-latitudes which contains 
the highest density of the conventional upper air network. Large differences 
in the northern hemispheric tropics are mainly due to the large number of fairly 
accurate cloud-tracked winds. Riven such large differences in the initial con- 
ditions, it is natural to expect differences in the forecasts. 


3. RESULTS 

Due to the limited space available here, we present the results of verifi- 
cation over India only. Table 3 gives the SI skill scores and rms error for 
sea level pressure and 500 mb geopotential height for forecasts starting from 
control and F6GE/M0NEX assimilations on 1200 GMT 7 July 1979. Both forecasts 
were verified against the NMC global analysis. In general, there Is a positive 
impact of FGGE/MONEX data; however, the differences after 24 h do not appear to 
be significant. The rms error for 500 mb geopotential height shows a negative 
impact for the first 36 h. It should also be pointed out that the location of 
the center of the monsoon depression in the FGGE/MONEX assimilation at 1200 GMT 
7 July 1979 does not agree with a careful hand analysis. 

In performing the numerical experiments, a number of problems were encoun- 
tered with the tropical analyses. The horizontal resolution (2.5® x 3®) of the 
prediction model was too coarse to accurately define the monsoon depression. 

The assimilation of geopotentia'i height data derived from satellite soundings 
generated gravity waves whose amplitudes were comparable to the meteorologically 
significant features we were attempting to predict. Gravity waves were parti- 
cularly troublesome in the 6 h first guess during the assimilation cycle. The 
difficulty in utilizing satellite temperature soundings in the tropics was also 
compounded by the precipitahle water contamination of the TIROS-N microwave 
retrievals (Phillips, 1980). The results of this study should be considered 
preliminary for these reasons. We plan to carry out further experiments using 
the official FGGE/MONEX database with improved assimilation and initialization 
techniques. It is only then that a more definitive statement can be made about 
the impact of FGGE/MONEX data on the predictability of the monsoon disturbances. 
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Table 2 . MONEX data over Bay of Bengal (10N-25N, 80E-100E) 


July 1 


July 2 


July 3 


July 4 


July 5 


July 6 


July 7 


July 8 


July 9 


July 


Dropwindsonde Aircraft Wisconsin 

Winds 


OOZ 

06Z 

12Z 

18Z 

OOZ 

06Z 

12Z 

18Z 

OOZ 

06Z 

12Z 

18Z 

OOZ 

06Z 

12Z 

18Z 


9 

3 


18 

34 

1 

1 

0 

1 23 

2 


6 11 

0 25 

1 

1 

1 9 

38 

1 


OOZ 

06Z 15 

12Z 7 

18Z 


2 7 

2 20 

2 


OOZ 

06Z 6 2 15 

12Z 3 3 9 

18Z 


OOZ 

06Z 8 

12Z 8 

18Z 


2 

7 18 

3 19 


OOZ 

06Z 17 

12Z 1 

18Z 


1 

4 27 

1 38 

1 


OOZ 

06Z 11 3 24 

12Z 3 1 49 

18Z 


OOZ 


1 


Enhanced 

Soundings 


no 


290 


188 


172 


135 
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Table 3. Skill score (SI) and root mean square error (RMS) over India (6N-26N, 
70E-100E) for control and FGGE/MONEX initial conditions of 120U GMT 
7 July 1979 for sea level pressure (SLP) and 500 mb geopotential 
height (Z500). 


Impact for India (6-26N I.C. 7/7/79, 12Z 



70-100E), 



SLP 

-SI 


Hours 

Control 

FGGE 

Impact 

12 

66.8 

62.6 

+4.2 

24 

60.6 

48.8 

+ 11.8 

36 

65.7 

63.1 

+2.6 

48 

52.8 

52.5 

+ .3 


SLP 

RMS 


12 

4.1 

2.2 

+ 1.9 

24 

4.0 

2.4 

+ 1.6 

36 

2.5 

3.3 

-.8 

48 

2.7 

1.9 

+ .8 


Z500 

SI 


12 

92.1 

85.1 

+7.0 

24 

99.4 

94.0 

+5.4 

36 

100.9 

94.1 

+6.8 

48 

104.1 

106.2 

-2.1 


Z500 

RMS 


12 

32.4 

43.9 

-9.5 

24 

53.6 

62.4 

-8.8 

36 

49.3 

51.2 

-1.9 

48 

48.4 

46.9 

+ 1.5 
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GLOBAL AND LOCAL FLUCTUr^i'IONS OF WINTER AND SUNNER 
SIMULATIONS WITH THE GLAS CLIMATE MODEL 

0. M. Straus and J. Shukla 


U INTRODUCTION 

Winter and suimner simulations have been carried out with an Improved ver- 
sion of the GLAS general circulation model. The model used for the simulations 
Is basically that presented by Halem et al. (1979). The changes that were made 
to the model code Involved the use an Improved method of computing the boundary 
layer fluxes, and a more realistic specification of the albedo of ^now and Ice 
covered surfaces. (See Shukla et a]. , 1980). The Winter Integration was Ini- 
tialized with data valid for OZ o1^ 1 January, 1975, and the summer Integration 
was Initialized with data valid for OZ of 15 June, 1979. Both Integrations 
were. carried out for 90 days, during which all the boundary conditions were 
prescribed to vary smoothly. The observed data consists of the NMC twice-daily 
analyses for the years 1963-1977. Each particular diagnostic quantity was com- 
puted from the model data and each of the 15 years of observations In precisely 
the same way, wherever possible. The reported observational results are aver- 
aged over the 15 winters or summers, as appropriate. 


2. STATIONARY WAVES IN THE GEOPOTENTIAL HEIGHT FIELD 

Results are presented for the stationary (time averaged) component of the 
geopotential height. Figure la shows the observed winter variance around a 
latitude circle of the stationary geopotential height, summed over wavenumbers 
1-4. The term "stationary" refers to a 90-day average commencing on January 1. 
The corresponding stationary ultra-long planetary wave variance for the winter 
model' simulation Is presented In Figure lb. 

The model results show good overall agreement with the observations, with 
both variances reaching a maximum In the mid-latitude upper troposphere. The 
model variance Is slightly -too weak up to the 250 mb level, above which It fails 
to show the decrease with height apparent In the observations. It Is Important 
to note that the stationary variance of this version of the GLAS climate model 
shown In Figure lb shows remarkable Improvement over the same variance calcul- 
ated for the previous version of the model, as presented by Straus and Shukla 
(1981, their Figure 6a). 

The observed summer variance of the stationary component of the gepoten- 
tlal height, summed over wavenumbers 1-4, is presented In Figure Ic. The cor- 
responding model variance Is shown In Figure Id. These quantities were computed 
In the same manner as for the winter season, except that the averaging period 
was defined to be 90 days starting on June 15. The simulated ultra-^ong wave 
variance agrees well with the observations with regard to the locat’ion and 
strength of the three maxima: one In mid-latitudes at about 300 nb, and two In 
the subtropics, at 850 mb and 200 mb. Discrepancies Include the fact that the 


^ The model results presented will be limited to the domain for which observa- 
tions were available, namely 20°N to 90**N. 
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simulated mid-latitude variance does not decrease above 300 mb* the slight 
northward displacement In the model of the subtropical maximum at 850 mb, and 
the weakness of the simulated subtropical maximum at 200 mo. 


3. LOW FREQUENCY PLANETARY WAVES 

Host of the zonal and time variance of the transient geopotentlal height 
field (exclusive of the annual cycle) Is contained In low frequency, planetary 
scale motions. More specifically, we are referring to fluctuations with periods 
of 7.5 to 90 days and consisting of zonal wavenumbers 1 to 4. 

The vertical and meridional structure of the low frequency planetary wave 
variance (LFPW) observed In winter Is shown In Figure Za^, This variance 
reaches a maximum slightly above the 300 mb level at about 64N, and shows a 
decrease above this level. The model variance (Figure 2b) has a similar struc- 
ture, except that the variance Increases up to at least 200 mb. In addition, 
the model variance Is somewhat too small throughout the northern nemispheric 
extratropical domain. As was the case for the stationary planetary waves, this 
result represents a great Improvement over the simulation of the previous 
version of the model. (See Straus and Shukla (1981), Figure 2a.) 

The observed summer LFPW variance Is shown In Figure 2c The structure of 
the summer variance Is quite similar to that In winter, except that the maximum 
has moved northward. The magnitude of the variance is much less In the summer. 
The model summer LFPW variance (Figure 2d) has the same overall dependence on 
latitude and height as observed, but again It does not decrease above 300 mb (as 
do the observations), and again the model variances are. In general, too small. 


4. LOCAL ANALYSIS OF VARIANCE 

A useful description of the fluctuations of any basic field Is afforded by 
the analysis of local variability discussed by Blackmon (1976). This method 
portrays the spatially local fluctuations of a particular field by the construc- 
tion of maps of Its RMS deviation In time. The variability due to barocllnic 
activity Is separated by filtering the data with a 'band-pass' filter which 
retains fluctuations with periods of 2.5 to 6 days. In the application of this 
method, an estimate of the annual cycle Is removed from the data before either 
the filtering or the computation of the RMS deviation. The band-pass RMS of 
the geopotentlal height Is closely associated with regions having a high fre- 
quency of cyclonic activity, and so Indicates the location of the 'storm tracks' 
(Blackmon, et al. , 1977). 

A comparison of the observed winter band-pass RMS of geopotentlal height 
(Figure 3a) with that of the model (Figure 3b) shows that there is excellent 
agreement In terms of both the location and strength of the major areas of 
cyclonic activity In the north central Pacific and western Atlantic. The 
observed summer band-pass RMS height field (Figure 3c) shows the two oceanic 
maxima that were present In the winter map. The Atlantic maximum has been 


^ For details regarding the calculation of observed and model LFPW variances, 
see Straus and Shukla (1981). 
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shifted to the east, and the maqnitude of the sunaner RMS is generally much 
smaller. The corresponding map ^or the model (Figure 3d) Indicates that the 
Atlantic maximum has been realistically simulated, both with respect to position 
and maqnitude. The Pacific maxlnum, however. Is too weak and Is located too 
far to the west. 
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Figure 1. Stationary planetary wave variance of the geopotential height field in the Northern 
Hemisphere, a) Observations, winter (units of 10^ m^); b) Model, winter (units of 
10^ m2); c) Observations, summer (units of 10^ m2); d) Model, summer (units of 10^ m2) 
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Figure 3. RMS deviation of the 500 mb band-pass geopotential height field, in units of meters. 

a) Observations, winter (contour interval is 5 m) ; b) Model, winter (contour interval 
is 5 m); c) Observations, sunrier (contour interval is 5m); d) Model, summer (contour 
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INFLUENCE OF SURFACE-ALBEDO IN SUBTROPICAL REGIONS 
ON JULY CIRCULATION 

Y. Sud and M. Fennessy 


A simulation study to examine the influence of surface-albedo on July cir- 
culation in subtropical regions is presented. The results are based on two 47- 
day integrations. In the first integration, called the control run, surface 
albedos were normally prescribed, whereas in the second integration, called the 
anomaly run, the surface albedo was modified in four regions: the Sahel in 
Africa, the Great Plains in the United States, the Thar Desert border in the 
Indian subcontinent, and Brazil in South America. Each run was started from 
observed initial conditions for June 15, 1979 based on NMC analysis. The sur- 
face albedo in each of the regions detailed in Table 1 was arbitrarily made 30%. 

The key results for examining the Charney mechanism are summarized in 
Table 2. 

The short-wave radiation absorbed by the surface was reduced in the anomaly 
runs. The percentage reduction was less than the percentage change in albedo 
because the reduced cloudiness partially mitigated the effect of the increased 
surface albedo. Simultaneously, there was an increase in the outgoing long-wave 
radiation at the surface. Altogether the surface radiation balance was reduced. 
This was largely balanced by a reduction in the surface fluxes of sensible heat 
and evaporation. However, the precipitation reduced even more than the evapor- 
ation. It implies even larger cooling due to lack of release of latent heat. 
Thus, the diabatic heating in these regions reduced substantially. Consistently, 
subsidence was induced and adiabatic heating of the order of the diabatic cool- 
ing was obtained. The region which did not follow the above pattern precisely 
was the Great Plains. There, the total cloudiness and total precipitation were 
virtually unchanged. Although there was some reduction in the convective pre- 
cipitation, it was offset by an increase in the large-scale precipitation. This 
region is somewhat unigue due to the influence of the Rocky Mountains on large- 
scale flows. Thus, the lack of response of precipitation to increased surface 
albedo in this one realization is not surprising. 

Our results differ in many ways from those of Charney et al . (1977). We 
did not find an increase in the solar radiation absorbed by the surface in the 
anomaly simulation. Perhaps, the excessive evaporation in their simulations 
was responsible for this effect. The evaporation and precipitation fields in 
our studies are more realistic compared to the observations. The reduction in 
precipitation based on 5-day averaged plots is not consistently maintained in 
any region except the Sahel, whereas in Charney et al. (1977) the differences 
were large and well maintained. Excessive evaporation in their simulation may 
have resulted in excessive moisture divergence due to subsidence, and enhanced 
its effect on cloudiness and precipitation. Thus, excessive evaporation is 
likely to amplify the albedo effect on precipitation. 

However, the thermal processes induced by an increase in the surface 
albedo were consistent with Charney's (1975) mechanism in all the regions, 
including the Great Plains, where the precipitation did not reduce. 


1 An extension of Charney's albedo experiments. 
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The results show that Charney's hypothesis that an arid region's large 
surface albedo through radiative cooling provides a feedback mechanism for the 
production of a drought, is well founded. 

Of the several large differences in the SLP and 500 mb geopotential height 
fields, the one persistent difference was over Central Eurasia (see for example 
Fig. 1). Elsewhere there were large differences, but these were not persistent 
over the two 15-day averages. The occurrence of distant circulation anomalies 
due to albedo anomalies in the tropics is very intriging and needs further 
investigation. 

In any event, the changes in the surface albedo of the anomaly run used in 
this study were smaller than those in the Charney et al « (1977) study. Since 
the effects are non-linear, and the results were di Inferential ly affected by 
horizontal advection, the Impact of albedo changes varies in the different 
regions. The purpose of this exercise was to re-examine Charney's hypothesis 
with a better GCM. This has been accomplished because the hypothesis is vali- 
dated in the Sahel and the Rajputana region, which belong to the Sahara-Arabian- 
Indo-Pakistani complex of deserts, as well as in subtropical Brazil. These 
regions are relatively free of strong advective effects, a condition essential 
for the Charney mechanism to be effective in reducing the rainfall. 
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Table 1 


Areal Extent of Albedo Anomaly Regions 



Region 

Number of 
Grid Points 

Latitude 

Longitude 

1 

Sahel 

22 

12-16 

17.5"W-47.5“E 

2 

Thar Desert 
Border 

4 

24"N-32®N 

67.5°E-77.5“E 

3 

Brazil 

10 

4°S-8“S 

4“S-20“S 

32.5'’W-37.5“W 

37.5“W-42.5‘’W 

42.5‘*W-47.5"W 

4 

Western Great 
Plains 

8 

32“N-48®N 

97.5“W-107.5‘’W 
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Impact of Albedo Anomaly on Simulated JULY 
Regional Differences (Anomaly - Control) 
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June 26, 1980 


DOES EL NINO AFFECT OUR WINTER CLIMATE? 
J. M. Wallace/University of Washington 


Atmospheric phenomena associated with the interannual variability of sea 
surface temperature in the equatorial Pacific are examined, with emphasis on 
vertical structure and teleconnections to middle latitudes. This paper is speci- 
ficially concerned with the interannual variability of seasonal means for the 
Northern Hemisphere winter during the period 1951-1978. Among the variables con- 
sidered are sea surface temperature in the equatorial Pacific, 200 mb height at 
stations throughout the tropics, a "Southern Oscillation Index" of sea-level 
pressure, the mean temperature of the tropical troposphere, precipitation at 
selected equatorial Pacific stations, and the Northern Hemisphere 700 mb height 
field. Selected statistics derived from surface data are also examined for the 
period 1910-45. Results are presented in the form of time series and correlation 
statistics for the variables listed above. 

Results concerning the relationships between sea surface temperature, sea 
level pressure and rainfall are consistent with the major conclusions of previous 
studies by J. Bjerknes and others. Fluctuations in mean tropospheric temperature 
and 200 mb height are showr* to vary simultaneously with equatorial Pacific sea 
surface temperature fluctuations, not only in the Pacific sector, but at stations 
throughout the tropics. The zonal ly symmetric component of these 200 height 
fluctuations is considerably larger than the Southern Oscillation in 1000 mb 
height, and the corresponding fluctuations in the mean temperature of the tropi- 
cal troposphere arc on the order of nearly 1 K. 

The tropical time series exhibit a well defined pattern of correlation with 
Northern Hemisphere 700 mb height. Warm episodes in equatorial Pacific sea 
surface temperature tend to be accompanied by below normal 700 mb height in the 
North Pacific and and the southeastern United States and above normal 700 mb 
height over western Canada. These correlations with mid-latitude 700 mb height 
tend tc be samewhat stronger than the corresponding correlations with mid-latitude 
surface temperature or pressure. The mid-latitude 700 mb heights are more strongly 
correlated with equatorial Pacific sea surface temperature and tropical 200 mb 
heights than with the Southern Oscillation Index. The pattern of mid-latitude 
700 mb height correlations resembles steady- state solutions of the linearized 
barotropic vorticity equation on a sphere, with forcing from the tropics. 

Much of the material presented in this talk is included in a manuscript by 
John D. Horel and John M. Wallace entitled "Planetary Scale Atmospheric Pheno- 
mena Associated with Interannual Variability of Sea Surface Temperature in the 
Equatorial Pacific" which has been submitted to Monthly Weather Review . 
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July 2, 1980 


AN ATTEMPT TO SEPARATE CHANGES DUE TO STATISTICAL 
SAMPLING FROM REAL CLIMATE CHANGE 

R. A« Madden/NCAR 


The variance of time- averaged meteorological data is assumed to consist of 
two components. One is a statistical sampling variability resulting from fast- 
varying weather that is present even without any climate change. The other 
results from shifts in mean values that occur on time scales long relative to 
that of weather variation. The first component is presumed to be unpredictable 
while the second is considered potentially predictable. This paper discusses an 
attempt to separate these two components in time series uf daily temperature for 
stations in the United States. 
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July 10, 1980 


SIMULATION OF THE JANUARY 1977 BLOCKING EVENT 
K. Miyakoda/NOAA/GFDL 


January 1977 was a month noted for its extraordinary weather over North 
America. The winter was dominated by two 1ong*sustaincd, large amplitude, 
stationary ridges positioned over the west coast of North America and the Ice- 
landic region of the Atlantic Ocean. A very intense trough reached deep into 
the eastern United States and caused one of the coldest Januaries on record, 
while western North America including Alaska remained warm. Several one-month 
integrations of various GCM's were conducted in order to test their ability to 
simulate the January 1977 blocking event. The GCM's used were of the finite 
different and spectral model types of the highest spatial resolution available 
at GFDL and ANMRC (Australian Numerical Meteorological Research Center). Two 
different initial conditions and two different external forcing conditions, 
were utilized by these models. In total, about 15 one-month integrations were 
performed. One month is the range just beyond the limit of predictability for 
deterministic forecasting, so the individual forecast, if successful, should 
be sufficiently close to reality. 

In some cases, a fairly accurate forecast was obtained, and considerable 
skill was recognized in the simulation of the 30-day evolution in terms of the 
5 day or 10 day mean flow fields, including the coldest period over the eastern 
United States which occurred in the last five days of January. It is concluded 
that a proper treatment of the subgrid-scale processes is essential for the 
simulation of the blocking phenomenon, in addition to the importance of a 
proper space resolution of the GCM. 
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July 17. 1980 


ON THE GLOBAL DISTRIBUTION OF THE HEAT SOURCES AND SINKS OF THE ATMOSPHERE 
AND THE THERMALLY-FORCED PLANETARY SCALE RESPONSE 

D. R. Johnson/University of Wisconsin 


The time-averaged structure of the Asiatic summer and winter monsoon cir* 
culatlon and their link with planetary scale transport processes are determined 
from a FGGE operational data set generated by the U. S. National Meteorological 
Center. The rotational and irrotatlonal components of the horizontal transport 
of mass and energy are determined through Isentropic diagnostics. 

The relations between 1) the global distribution of boundary energy flux 
through the earth's surface and the top of the atmosphere. 2) heat sources and 
sinks Internal to the atmosphere and 3) planetary scale mass and energy trans- 
port processes within the Isentropic stratification of the atmosphere are broadly 
established from transport theory and results. For the months of January, April, 
July and October of 1979, the global distribution of diabatic heating Is deter- 
mined indirectly and linked to the time-averaged mass and energy transport within 
large scale Hadley-type and Walker-type circulations. The results verify the 
global nature of monsoonal circulations that are linked to the planetary scale 
of differential heating in the form of large scale Hadley-type and Walker-type 
circulations. The primary planetary source region of energy In the region of 
Indonesla-Phll Ippines-Southeast Asia is linked through quasi -horizontal mass 
transport with primary sink regions of energy In the two circumpolar vortices, 
the Sahara and subtropical anticyclonic circulations. This primary center of the 
source of energy In the Indonesia- Southeast Asia area moves seasonally from one 
hemisphere to the other. Net heating within a given region and net cooling In 
another region In association with differences in boundary flux through the 
earth's surface and the top of the atmosphere result In net energy transport from 
the source to the sink region, a result basic to thermally- forced circulations. 

The net energy transport Is realized from the condition that more energy Is trans- 
ported by a branch of a circulation In higher-valued isentropic layers from the 
heat source to the heat slnkethan Is returned by a branch of circulation within 
lower-valued Isentropic layers. At the same time no net mass change occurs. 

Thus a steady time-averaged circulation that Is enfcedded In the Isentropic struc- 
ture may be realized. 

The winter monsoon circulation In Asia links radiative cooling over Asia to 
sensible and latent heat release within mid-latitude barocllnic waves and also 
to latent heat release In the Indonesia, Philippine and New Guinea regions 
through horizontal transport of mass and energy. The role of mid-latitude baro- 
cllnic waves In providing for mass and energy exchange between polar and subtrop- 
ical regions within the time-averaged Isentropic structure is briefly discussed. 
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July 31. 1980 


COMPARISONS OF VISIBLE AND INFRARED DISTRIBUTIONS 
OF GLOBAL CLOUD COVERS 

M. T. Chahine/JPL 


Day and night mapping of the global distributions of the horizontal cloud 
cover and the corresponding cloud- top pressure levels can be obtained from the 
same set of Infrared radiance data used to retrieve dear-column temperature 
profiles. General formulation of the problem Is presented with Illustrations 
for the simple case of a single layer of non-reflecting clouds. Experimental 
verifications are obtained using 15 urn data measured by the NOAA-VTPR Infrared 
sounder. 

The upwelling radiance from a planetary atmosphere Is a function of the 
thermal state of the atmosphere, the concentration of radlatively active gases, 
and the extents, heights, and radiative transfer properties of clouds and 
aerosols. Thus, In principle. It should be possible to recover useful Informa- 
tion about the physical and chemical structure of an atmosphere from analysis 
of the upwelling radiance. However, the problem In analyzing such data lies 
In finding ways to uncouple the effects of these variables and retrieve the 
true values of each unknown parameter separately. By treating the cloud effects 
as short period oscillations over the clear column radiance, an analytical 
method was developed by Chahine to retrieve dear-column vertical temperature 
profiles from radiance measurements made In the presence of clouds. The 
uncoupling of the effects of clouds Is carried out analytically without any ^ 
priori Information about the amounts, heights and optical properties of the 
clouds In the fields of view. Once the dear-column temperature profiles are 
determined the same radiance data could then be used to determine the heights, 
amounts, and radiative transfer properties of clouds. 

The determination of the dear-column temperature profiles from the VTPR 
data requires a priori knowledge of the surface temperature Tj. We obtained 
Tj from the NOAA surface analysis. We investigated the effects of errors In 
the assumed surface temperature on the accuracy of the values of the effective 
cloud cover JT which Is the product of the geometrical cover N and the cloud 
emissivlty and the mean cloud top pressure, P^. We concluded that the 
effects of an error In Tj of ±2K on N^ and Pc are small, especially for 
Pc < 700 mb. Correction for the effects of water vapor on the atmospheric 
transmission functions were made before generating the dear-column radiances. 

It should be noted here that while the determination of the clear column 
radiance Is obtained without any assumptions about the properties of clouds, 
the determination of the amount and height of clouds requires the use of cloud 
models. In the case of the VTPR data we assumed that the difference between 
the reconstructed clear column radiance and the radiance measured In a given 
field of view Is due to the presence of a single layer of non-reflecting cloud. 
We applied the nonreflecting cloud model to analyze radiance data from the 
NOAA-VTPR sounder for a period of one week from January 1-7, 1975. The VTPR 
global cloud distributions were calculated In 1977. The results were then 
averaged for a grid size of 4° latitude by 5° longitude, and only the averaged 
results were stored on a magnetic tape for subsequent comparison with other 
cloud maps to be obtained from other sources. A typical comparison of the 
results Is shown In Figures 1, 2, and 3. 
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Figure 1 shows contours of cloud amounts for a region across the Pacific 
Ocean between 40“N-30®S snd 75®W-255“W. Fig. 2 shows contours of cloud amounts 
derived from computations made by J. Sadler of the University of Hawaii for the 
same period of time and the region from photographs obtained from the Vidicon 
cameras of NOAA's satellite. 

The results shown in Fios. 1 and 2, therefore, compare asynoptic infrared 
cloud maps and synoptic visual maps. Consequently only persisting cloudiness 
appears to be common between the two cloud maps. For this reason the zonal ly 
averaged values shown in Fig. 3 give a more realistic comparison between the 
amount of clouds observed in the infra»*ed and visible. The conclusion that the 
effective infrared cloud amount N~ = Ne is smaller than the cloud amount ob- 
served in the visible could be due to th#» facts that the cloud emissivity In 
the 15 ym is less than one, and the VTPR sounding channels are not sensitive 
to detecting low level clouds below 800 mb. The average cloudiness for the 
region shown in Figs. 1 and 2 is 0.39 for the infrared and 0.52 for the visible, 
and the ratio of the infrared to the visible cloud cover is ~ 0.75. 
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Fin. 1. Contours of cloud amounts in decimals derived from 15 um satellite 
data for the period of January 1-7, 1975. 
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Fiq. Contours of cloud amounts in decimals derived from visible data by 
Sadler for the period of January 1-7, 1975. 
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June 23. 1981 


DYNAMICS OF SUDDEN STRATOSPHERIC WARMINGS 
J. R, Holton/The University of Uashinqton 


The complex of events associated with sudden warmings in the Northern 
Hemisphere polar winter stratosphere provides perhaps the best natural example 
of planetary scale wave-mean flow interaction processes. Recent work by Andrews 
and McIntyre (1976) and others (cf., review by Holton, 1980) has provided new 
diagnostic methods for studying both observed and simulated stratospheric 
warmings. With the aid of these methods rapid pro'jress is being made in under- 
standing the origin and evolution of sudden warmings. 

In this lecture I first review the observed behavior of planetary waves 
and the mean flows during sudden warmings. Next I review wave-mean flow Inter- 
action theory in a quasi-qeostrophic framework and derive the so-called "non- 
acceleration" theorem. By an elementary transformation of the Eulerian mean 
equations I show that the net forcing of the mean flow by the waves may be 
expressed in terms of the divergence of the Eliassen-Palm flux, a vector in 
the meridional plane whose northward component is proportional to the merid- 
ional heat flux. This next forcing can alternatively be written as the poleward 
flux of quasi-qeostrophic potential vorticity. 

In the final portion of the lecture I discuss the similarities of the 
transformed Eulerian analysis to the more fundamental Lagrangian mean formula- 
tion, and show how the net forcing of the mean flow by the waves can be Inter- 
preted in terms of the divergence of the wave radiation stress. 


REFERENCES 

Andrews, D. G., and M. E. McIntyre, 1976: Planetary waves in horizontal and 

vertical shear: The generalized Eliassen-Palm relation and the zonal 

mean acceleration. J. Atmos. Sci. , 33 , 2031-2048. 

Holton, J. R., 1980: The dynamics of sudden stratospheric warmings. Ann. 

Rev. Earth Planet. Sci., 8, 169-190. 
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THE BEHAVIOR OF EDDY HEAT FLUXES ON SHORT TIME SCALES 
P. H. Stone^ /Massachusetts Institute of Technology 


July 7, 1981 


NMf. data from 30®N to 70®N for three Januaries are used to calculate time 
series of the tropospheric eddy flux of sensible heat and of the baroclinic 
stability, as measured by the difference between the zonal wind shear and the 
critical value of the shear in two-level models. On synoptic time scales the 
variance of the eddy flux is about 35% of the mean, and the variance of the 
stability parameter is about 2 m/s. Typically, the zonal wind shear is about 
2 m/s greater than the critical value of about 10 m/s. Time-lagged correla- 
tions between the two series show a highly significant negative correlation 
for short time lags, indicative of "free" variations in the eddy flux on short 
time scales. The negative correlation peaks at about -0.46 when the stability 
parameter lags about one half day behind the eddy flux. The autocorrelation 
of the eddy flux is compared with that expected for various low-order autore- 
gressive processes. 

Application of a Bayesian information criterion indicates that the short- 
term behavior can be represented best by a second-order model, with a damping 
time of about 0.8 day and periodicity of about 5 days. Because of the rela- 
tively long periodicity, the behavior is also represented guite well by a 
first-order or red-noise model, with an autocorrelation time of about 0.8 day. 
These results support theoretical models for the finite-amplitude eguilibra- 
tion of unstable baroclinic waves which yield a Stuart-Hatson type of eguili- 
bration equation. Such models simulate well the observed short-term fluctua- 
tions of the eddy flux of white-noise forcing is added to them, and they 
predict that the autocorrelation time is just one half the e-folding time of 
the linearly most unstable wave. The fluctuations of the stability parameter 
are analyzed in a similar way, and are found to be matched best by a red-noise 
model, with an autocorrelation time of about two days. However, the cross- 
correlation indicates that the forcing by the flux should be included in the 
model for the stability parameter. When this is done, a better fit is obtained 
for the autocorrelation function of the stability parameter, and an excellent 
fit for the cross-correlation, if the relaxation time for the stability para- 
meter is about 1.5 days. The red-noise models are used to calculate that the 
variance of the January mean eddy flux is 9% of the mean, and the variance of 
the January mean stability is 0.5 m/s. 


I Also authored with S. J. Ghan, D. Spiegel, and S. Rambaldi. 
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July 14, 1981 


NORMAL MODE INITIALIZATION AT ECMWF 
C. Temperton/European Centre for Medium Range Weather Forecasts 


Normal mode initialization is an essential component of the data assimila- 
tion/anal.vsis/forecast cycle at ECMWF; without it, unrealistic high-frequency 
oscillations soon degrade the quality of the analyses. 

In this seminar, we review the procedure for obtaining t'.i normal modes of 
a multi-level global qridpoint model (linearized about a state of rest), and 
show how they can he used in both linear and nonlinear initialization techniques. 
The linear approach does not give very good results, but the nonlinear procedure 
is highly successful in preventing contamination of the subsequent forecast by 
spurious oscillations. Moreover, it is computationally very efficient. 

However, nonlinear normal mode initialization as currently applied is not 
entirely free of defects. It is difficult to take into account the highly non- 
linear forcing by diabatic processes, and - partly as a consequence - the 
changes made to the original analysis (especially the surface pressure field) 
can he uncomfortably large. We conclude by outlining some possible approaches 
towards improving the procedure. 
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July 21, 1981 


LARGE-SCALE BOUNDARY LAYER WINDS 
J. A. Young/The University of Wisconsin 


Recent research conducted at Wisconsin on the behavior and dynamics of 
synoptic-scale wind systems in the lower atmosphere over the tropical oceans 
has been reviewed. Our knowledge about these winds has been substantially 
improved by recent observational studies such as the Monsoon Experiment and 
the Global Weather Experiment (1979), and by developing models of cumulus- 
filled planetary boundary layers. For example, one year's data for the entire 
low-level monsoon circulation was provided by a geostationary satellite, 
allowing calculation of kinematic and dynamic fields discussed below. In 
addition, gust-probe research aircraft data allowed quantitative estimation 
of scale-dependent turbulent fluxes as well as vertical profiles through the 
depths of the planetary boundary layer (PBL). These measurements support the 
need for a multi-layer approach to modeling these flows. The results presented 
here demonstrated 

a) the usefulness of the FGGE Il-b data for diagnostic dynamics studies; 

b) the way that boundary layer processes determine wind response to 
pressure forces; 

c) the expectation that these understandings can be applied to the 
problem of data assimilation and verification influencing the surface 
winds; 

d) the expectation that the resulting Ill-b data sets, with appropriate 
PBL assimilations, will allow a complete air-sea interaction history 
for the global oceans for all of 1979. 

The satellite wind fields for the summer monsoon show a dramatic intensi- 
fication between Hay and July, with an absolute vorticity distribution which 
indicates (a) non-conservation along streamlines, and (b) possible inertial 
instability between the equator and 8N. The wind fields appear to be of 
sufficient quality to calculate terms in the equations of motion: the mean 

pressure gradient force field was calculated as a residual for the entire 
western Indian Ocean, showing how it changed as an accelerating air parcel 
tra.iectory was followed. The constraint that this force must be irrotational 
allows quantitative estimation (by relaxation of Poisson equations) of: (a) 
the pressure field and (b) the dominant part of the friction field at 1 km 
altitude. 

The difficulties of PBL dynamics include the facts that the flow is 
highly baroclinic and accelerated in the mean and the turbulence processes 
depend upon scale and the role of cumulus clouds. These features were clearly 
seen in profiles of mean flow and fluxes taken from research aircraft data; 
it was noted that the stress component along the mean flow decreased from 
the surface, but could experience a reversal or secondary maximum when the 
cloud-scale eddy fluxes were included. The time scale of frictional decay 
of the jet flow was in the range of 0. 8-3.0 days. 
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A predictive PBL model with wind and mass adjustments has been developed 
which attempts to take into account all of these observed features: it repre> 

sents the PBL mean, flux and cloud activity profiles at three levels whose 
variation is predicted physically. The model allows the effective friction 
layer to be shallower than the thermal or moisture boundary layers, as ob- 
served. Results were shown for trade-wind simulations indicating dynamically 
determined wind structures, sensitivity to cloud flux processes, and depend- 
ence of clouds and wind profiles upon sea surface temperature. The results 
indicated good potential for application to Ill-b data assimilations at low 
levels* 
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July 28, 1981 


THE ROLE OF THE OCEAN IN THE TRANSIENT CLIMATIC RESPONSE TO 
INCREASING ATMOSPHERIC CC^ 

K, Bryan /NOAA/Geophysi cal Fluid Dynamics Laboratory 


Many studies have been made of the role of the oce is in absorbing excess 
atmospheric C 02 * Less attention has been paid to the role of the ocean in any 
transient response of climate to rising atmospheric CO 2 . Due to its enormous 
heat capacity the ocean can delay a warming due to radiative changes in the 
atmosphere associated with CO 2 . The size of the delay depends on the rate of 
penetration of heat from the ocean surface to lower levels. 1 he more effective 
downward pathways are, the slower will be the climatic response to rising CG 2 . 

At GFOL a coupled ocean>atmosphere model with an idealized geometry is 
used to study this problem quantitatively. First, climatic equilibria were 
determined for normal and 4X normal CO 2 concentrations in the atmosphere. This 
was done by a relaxation procedure described in previous studies. Next, the 
normal CO 2 equilibrium climate is perturbed by a step function, 4-fold increase 
of CO 2 . The response is studied by a detailed synchronous integration of the 
coupled ocean atmosphere system. Let the normalized response, R, be defined 
as. 


R = 



where T is the sea surface temperature, and T» and Tq correspond to the high 
and low CO 9 equilibrium climates. R rises much faster initially in the tropi- 
cal latitudes, but after a 25 year interval it is nearly uniform at all latitudes 
with a value between 0.6 and 0.7. The tendency to become uniform with respect 
to latitude is much greater if T is defined as the surface air temperature. 


Figure 1 shows the temperature anomaly after 10 years with respect to the 
normal CO 2 equilibrium climate for zonal average in the atmosphere and the 
ocean. Note that the penetration of heat is greatest in the subpolar region 
whore water mass formation normally takes place. Little penetration takes 
place in the tropical region which is dominated by upwelling. 


REFERENCES 

Bryan, K., F. 6 . Komro, S. Manabe, and M. J. Spelman, 1981: Transient response 
to increasing atmospheric CO 2 . (submitted to Science ). 
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August 4, 1981 


TRANSITION TO TURBULENCE IN SHEAR FLOWS 
S. A. Ors^aq/Massachusetts Institute of Technology 


A new three-dimensional linear instability, called catalytic instability, 
of two-dimensional (or axisymmetric) finite-amplitude flows is described. This 
instability seems responsible for transition to turbulence in the classical 
shear flows, including plane Poisenille flow, plane Couette flow, pipe Poisenille 
flow, boundary layers, jets, and mixing layers. Good agreement with experiment 
is achieved regarding critical Rynolds numbers and the structure of early 
transitional spots. The physics of the instability involves a subtle interplay 
of vortex stretching and tilting and is characterized by energy transfer from 
the mean (parallel) flow to the three-dimensional disturbance with the two- 
dimensional flow acting as a catalyst for this transfer. The subsequent evolu- 
tion of the evolving three-dimensional flow is nonlinear and chaotic. The 
resulting turbulence is in good agreement with experiment. For further details 
see Orszag and Patera, 1980; also Transition and Turbulence , ed. by R. E. Meyer, 
Academic, 1981, p. 127; also Phys. Rev. Letters , Sept. 21, 1981, and J. Huid 
Mech. , in press. 
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